Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,421,120 Δ1 bp noncoding (60/121 nt) rrfF ← 5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,421,1190G.94.8% 276.6 / 14.3 115noncoding (61/121 nt)rrfF5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base . (49/60);  total (52/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-02

Only 100 of 129 total aligned reads displayed.

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207
                                                                                            |                                                                                        
cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          >  1:703768/1‑93 (MQ=255)
cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          >  1:3039237/1‑93 (MQ=255)
 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         <  1:4694263/93‑2 (MQ=39)
 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  2:5774927/1‑92 (MQ=255)
 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  1:2611368/1‑92 (MQ=39)
  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:1563974/1‑91 (MQ=255)
  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:1402407/1‑91 (MQ=255)
  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:1789481/1‑91 (MQ=38)
  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGAGCTAcggcgtt                                                                                        >  1:1022242/1‑91 (MQ=255)
   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       >  1:2113350/1‑93 (MQ=255)
   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       <  1:5447988/93‑1 (MQ=37)
    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      <  2:5339668/93‑2 (MQ=25)
    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      <  1:690018/93‑2 (MQ=25)
     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                                                    >  1:255358/1‑93 (MQ=255)
     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                                                    >  1:1284828/1‑93 (MQ=255)
      tCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAc                                                                                   >  2:3208144/1‑93 (MQ=255)
       ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTGCTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                                                  <  1:5935607/93‑1 (MQ=255)
       ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                                                  <  1:6076548/93‑1 (MQ=12)
        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                                                                 <  1:3583228/93‑1 (MQ=9)
        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                                                                 >  1:3098987/1‑93 (MQ=12)
        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                                                                 >  1:2676474/1‑93 (MQ=16)
         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                <  2:5317095/93‑1 (MQ=9)
         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                <  1:2155447/93‑1 (MQ=9)
         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                >  2:1700805/1‑93 (MQ=12)
          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                               >  2:1106583/1‑93 (MQ=16)
          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                               >  1:4754346/1‑93 (MQ=12)
          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                               >  2:744479/1‑93 (MQ=16)
          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                               <  2:835130/93‑1 (MQ=9)
           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTCCTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                                              <  2:5455358/93‑1 (MQ=9)
           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCGTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                                              >  2:4926531/1‑93 (MQ=255)
           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                                              >  2:3736350/1‑93 (MQ=16)
           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                                              >  1:5151176/1‑93 (MQ=9)
           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                                              >  1:4887814/1‑93 (MQ=16)
           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                                              >  2:4266416/1‑93 (MQ=16)
            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTgg                                                                             >  1:125350/1‑92 (MQ=16)
            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                                             >  1:771839/1‑93 (MQ=12)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCGCACTACCATCGGCGCTACGGCGTTTCACTTTTGAg                                                                            >  2:5153433/1‑93 (MQ=11)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                                                            >  1:1817180/1‑93 (MQ=9)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                                                            <  1:2524135/93‑1 (MQ=9)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                                                            >  1:291461/1‑93 (MQ=16)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                                                            >  2:3921995/1‑93 (MQ=16)
             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAg                                                                            <  2:4273798/93‑1 (MQ=9)
              aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGt                                                                           >  2:4045441/1‑93 (MQ=16)
              aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGt                                                                           <  2:1547098/93‑1 (MQ=9)
               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                                                          >  2:1092349/1‑93 (MQ=16)
               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                                                          >  1:2882613/1‑93 (MQ=16)
               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                                                          >  1:5055150/1‑93 (MQ=16)
               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                                                          >  1:3294038/1‑93 (MQ=16)
                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                                                         >  1:5375275/1‑93 (MQ=16)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTGCGGCGTTTCACTTCTGAGTTCg                                                                        <  2:2869136/93‑1 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        >  2:5074220/1‑93 (MQ=12)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        >  2:3833883/1‑93 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        >  1:5699108/1‑93 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        <  2:2366149/93‑1 (MQ=9)
                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACGCTCGCATGGGGAGCCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                                        <  2:5685310/93‑1 (MQ=11)
                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                                                       >  1:4858032/1‑93 (MQ=9)
                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                                                       >  1:3055950/1‑93 (MQ=9)
                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                                                       >  1:1970988/1‑93 (MQ=12)
                   tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                                                      <  2:2885359/93‑1 (MQ=9)
                   tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                                                      >  1:1725799/1‑93 (MQ=12)
                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                                                                     <  1:6276541/93‑1 (MQ=0)
                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                                                                     <  2:5600796/93‑1 (MQ=0)
                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                                                                     >  1:1154298/1‑93 (MQ=0)
                    ttGTTTTATTTGATGTCTGCCAGTCCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTCCGGCGTTTCACTTCTGAGTTCGGCa                                                                     <  2:3058236/93‑1 (MQ=1)
                                                                   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  2:4964186/93‑1 (MQ=9)
                                                                   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  2:1794561/93‑1 (MQ=9)
                                                                   cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      >  1:2702378/1‑93 (MQ=0)
                                                                    ccACACTACCATCGGCGCTATGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                     <  1:4203400/93‑1 (MQ=255)
                                                                     cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:1059927/93‑1 (MQ=12)
                                                                     cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:962756/93‑1 (MQ=12)
                                                                     cacaCTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:3359044/93‑1 (MQ=12)
                                                                      acacTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                   <  2:1909311/93‑1 (MQ=19)
                                                                       cacTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  <  2:119737/93‑1 (MQ=19)
                                                                       cacTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  >  1:2383779/1‑93 (MQ=19)
                                                                        acTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  2:1554006/93‑1 (MQ=255)
                                                                        acTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  2:3160725/93‑1 (MQ=255)
                                                                        acTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  2:2942558/93‑1 (MQ=255)
                                                                         cTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                <  1:5435902/93‑1 (MQ=255)
                                                                         cTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                >  1:4373139/1‑93 (MQ=255)
                                                                          tACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               <  1:5944513/93‑1 (MQ=255)
                                                                          tACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               <  1:3571448/93‑1 (MQ=255)
                                                                           aCCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCt              <  2:2352940/93‑1 (MQ=255)
                                                                            ccATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTg             <  1:1608026/93‑1 (MQ=255)
                                                                            ccATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTg             >  2:2244692/1‑93 (MQ=255)
                                                                             cATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc            <  2:1213051/93‑1 (MQ=255)
                                                                             cATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc            <  2:414512/93‑1 (MQ=255)
                                                                             cATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc            >  1:5480120/1‑93 (MQ=255)
                                                                              aTCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           >  2:470324/1‑93 (MQ=255)
                                                                              aTCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           <  1:64425/93‑1 (MQ=255)
                                                                                 ggCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGaa        <  1:280420/93‑1 (MQ=255)
                                                                                  gcgcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  1:812163/93‑1 (MQ=255)
                                                                                    gcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     <  2:6069915/93‑1 (MQ=255)
                                                                                    gcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     >  1:274350/1‑93 (MQ=255)
                                                                                    gcTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     <  1:2162413/93‑1 (MQ=255)
                                                                                      tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACttt  <  2:6049238/93‑2 (MQ=255)
                                                                                      tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   <  1:4104696/93‑1 (MQ=255)
                                                                                      tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   <  2:598632/93‑1 (MQ=255)
                                                                                      tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   <  2:911544/93‑1 (MQ=255)
                                                                                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  >  1:82623/1‑93 (MQ=255)
                                                                                       aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  >  1:492508/1‑93 (MQ=255)
                                                                                            |                                                                                        
CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAA  >  NC_000913/3421027‑3421212
                                                                                             |                                                                                            
CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGT                                                                                              >  SRR4302533.2627394/1‑94 (MQ=36)
 AGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTT                                                                                             <  SRR4302533.4722708/94‑1 (MQ=34)
 AGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGAGCTACGGCGTT                                                                                             >  SRR4302533.1028252/1‑94 (MQ=32)
  GTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTT                                                                                            >  SRR4302533.4177852/1‑94 (MQ=31)
  GTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTT                                                                                            >  SRR4302533.2126384/1‑94 (MQ=31)
   TCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTC                                                                                          <  SRR4302533.5480808/94‑1 (MQ=28)
    CTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCA                                                                                         >  SRR4302533.1292557/1‑94 (MQ=22)
    CTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCA                                                                                         >  SRR4302533.256920/1‑94 (MQ=22)
    CTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCA                                                                                         <  SRR4302533.694098/94‑1 (MQ=22)
     TTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAC                                                                                        <  SRR4302533.2276509/94‑1 (MQ=14)
       CCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTGCTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTT                                                                                      <  SRR4302533.5971493/94‑1 (MQ=6)
       CCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTT                                                                                      >  SRR4302533.2692886/1‑94 (MQ=7)
       CCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTT                                                                                      >  SRR4302533.3118039/1‑94 (MQ=7)
       CCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTT                                                                                      <  SRR4302533.5971492/94‑1 (MQ=7)
       CCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTT                                                                                      <  SRR4302533.6113300/94‑1 (MQ=7)
        CGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTC                                                                                     <  SRR4302533.3605296/94‑1 (MQ=4)
         GACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCT                                                                                    <  SRR4302533.2168719/94‑1 (MQ=4)
         GACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCT                                                                                    >  SRR4302533.4783164/1‑94 (MQ=4)
         GACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCT                                                                                    <  SRR4302533.761172/94‑1 (MQ=4)
          ACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTG                                                                                   >  SRR4302533.4917429/1‑94 (MQ=4)
          ACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTG                                                                                   >  SRR4302533.5182293/1‑94 (MQ=4)
          ACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTG                                                                                   <  SRR4302533.2602555/94‑1 (MQ=4)
          ACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTG                                                                                   <  SRR4302533.2799566/94‑1 (MQ=4)
           CTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGG                                                                                  >  SRR4302533.126135/1‑94 (MQ=4)
           CTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGA                                                                                  >  SRR4302533.776379/1‑94 (MQ=4)
            TGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAG                                                                                 >  SRR4302533.1828273/1‑94 (MQ=4)
            TGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAG                                                                                 >  SRR4302533.293236/1‑94 (MQ=4)
             GAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGT                                                                                <  SRR4302533.2539683/94‑1 (MQ=4)
              AGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTT                                                                               >  SRR4302533.5085709/1‑94 (MQ=4)
              AGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTT                                                                               >  SRR4302533.3314369/1‑94 (MQ=4)
              AGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTT                                                                               >  SRR4302533.2900345/1‑94 (MQ=4)
               GCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTC                                                                              >  SRR4302533.867885/1‑94 (MQ=4)
               GCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTC                                                                              >  SRR4302533.5407578/1‑94 (MQ=4)
                CCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCG                                                                             >  SRR4302533.5733526/1‑94 (MQ=4)
                CCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCG                                                                             <  SRR4302533.1485379/94‑1 (MQ=4)
                 CTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGG                                                                            >  SRR4302533.3010988/1‑94 (MQ=4)
                 CTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGG                                                                            >  SRR4302533.3074687/1‑94 (MQ=4)
                 CTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGG                                                                            >  SRR4302533.4887474/1‑94 (MQ=4)
                 CTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGG                                                                            >  SRR4302533.1983048/1‑94 (MQ=4)
                  TTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGC                                                                           >  SRR4302533.1736272/1‑94 (MQ=0)
                   TTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCA                                                                          >  SRR4302533.1161216/1‑94 (MQ=0)
                    TTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAT                                                                         <  SRR4302533.6314020/94‑1 (MQ=0)
                                             CCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGGGGTTTtcctttttctttttgggggggggggggggggggggggcgcgc                                                 >  SRR4302533.4444392/1‑52 (MQ=0)
                                                                    CCACACTACCATCGGCGCTATGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTC                         <  SRR4302533.4228829/94‑1 (MQ=0)
                                                                    CCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTC                         <  SRR4302533.1945644/94‑1 (MQ=0)
                                                                    CCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTC                         >  SRR4302533.190925/1‑94 (MQ=0)
                                                                     CACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCT                        <  SRR4302533.968414/94‑1 (MQ=1)
                                                                      ACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTA                       >  SRR4302533.3558221/1‑94 (MQ=4)
                                                                      ACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTA                       >  SRR4302533.2398515/1‑94 (MQ=4)
                                                                        ACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAT                     >  SRR4302533.4399538/1‑94 (MQ=9)
                                                                        ACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAT                     <  SRR4302533.4876469/94‑1 (MQ=9)
                                                                         CTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATC                    <  SRR4302533.5468636/94‑1 (MQ=11)
                                                                         CTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATC                    <  SRR4302533.5468637/94‑1 (MQ=11)
                                                                          TACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCT                   <  SRR4302533.3593458/94‑1 (MQ=20)
                                                                          TACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCT                   <  SRR4302533.5980463/94‑1 (MQ=20)
                                                                           ACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTG                  >  SRR4302533.5613448/1‑94 (MQ=22)
                                                                            CCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGC                 <  SRR4302533.1617782/94‑1 (MQ=24)
                                                                            CCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGC                 >  SRR4302533.6015360/1‑94 (MQ=24)
                                                                            CCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGC                 >  SRR4302533.5513160/1‑94 (MQ=24)
                                                                              ATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCG               >  SRR4302533.1843606/1‑94 (MQ=29)
                                                                              ATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCG               <  SRR4302533.64842/94‑1 (MQ=29)
                                                                              ATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCG               >  SRR4302533.2263040/1‑94 (MQ=29)
                                                                                 GGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGCAC            >  SRR4302533.122697/1‑94 (MQ=29)
                                                                                 GGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAC            <  SRR4302533.282131/94‑1 (MQ=37)
                                                                                  GCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACT           <  SRR4302533.3591968/94‑1 (MQ=39)
                                                                                  GCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACT           <  SRR4302533.816922/94‑1 (MQ=39)
                                                                                   CGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTT          >  SRR4302533.276019/1‑94 (MQ=41)
                                                                                   CGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTT          >  SRR4302533.2115103/1‑94 (MQ=41)
                                                                                    GCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTT         <  SRR4302533.5066013/94‑1 (MQ=44)
                                                                                    GCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTT         <  SRR4302533.2175719/94‑1 (MQ=44)
                                                                                     CTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTA        <  SRR4302533.6245538/94‑1 (MQ=46)
                                                                                      TACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA       >  SRR4302533.495452/1‑94 (MQ=52)
                                                                                      TACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA       >  SRR4302533.83147/1‑94 (MQ=52)
                                                                                      TACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA       <  SRR4302533.4129662/94‑1 (MQ=52)
                                                                                         GGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCT    <  SRR4302533.1451819/94‑1 (MQ=60)
                                                                                           gcGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTA   <  SRR4302533.5452355/92‑1 (MQ=60)
                                                                                           gcGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTA   >  SRR4302533.1578757/3‑94 (MQ=60)
                                                                                           gcGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTA   <  SRR4302533.741784/92‑1 (MQ=60)
                                                                                           gcGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTA   <  SRR4302533.3970114/92‑1 (MQ=60)
                                                                                           gcGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTA   <  SRR4302533.2671500/92‑1 (MQ=60)
                                                                                            CGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAA  >  SRR4302533.3962153/1‑94 (MQ=60)
                                                                                            CGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAA  >  SRR4302533.4551674/1‑94 (MQ=60)
                                                                                             |                                                                                            
CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAA  >  NC_000913/3421027‑3421212

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: