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 4,603,785 A→G noncoding (71/87 nt) leuP ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,603,7850AG98.8% 489.2 / ‑4.5 176noncoding (71/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (97/74);  total (98/78)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 199 total aligned reads displayed.

CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAA  >  NC_000913/4603697‑4603874
                                                                                            |                                                                                          
ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTg                                                                                       >  1:3159587/1‑93 (MQ=255)
     aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATAGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAcc                                                                                   >  1:1741595/1‑93 (MQ=37)
      gCTAGCGCGTCTACCAATTCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAAcccc                                                                                 <  1:1662271/93‑1 (MQ=255)
      gCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAccc                                                                                  <  2:1490658/93‑1 (MQ=38)
       cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAcccc                                                                                 <  1:4618449/93‑1 (MQ=38)
       cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAcccc                                                                                 >  1:4103334/1‑93 (MQ=255)
        tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCa                                                                               <  2:6117242/93‑1 (MQ=255)
         aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGCGAGGGGGGGGGGCTTGAACCCCCa                                                                               <  2:1531837/93‑1 (MQ=38)
          gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCAc                                                                              >  2:2748960/1‑93 (MQ=255)
           cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTGGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACg                                                                             <  2:5012346/93‑1 (MQ=38)
            gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGt                                                                            >  1:3551001/1‑93 (MQ=255)
            ctgTCTACCAATTCCGCCACCTTGACACAGTCATCTTACTTTTTTTGATACCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTc                                                                           <  2:1238115/91‑1 (MQ=37)
            ccgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTc                                                                           <  2:4676407/92‑1 (MQ=255)
             cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGGCTTGAACCCCCACGTc                                                                           <  2:2932993/93‑1 (MQ=255)
              gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTcc                                                                          >  1:5885707/1‑93 (MQ=255)
              gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTcc                                                                          >  1:2700141/1‑93 (MQ=255)
              gTCTACCAATCCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCCGGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTcc                                                                          <  2:4797330/93‑1 (MQ=37)
                cTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGCTTGAACCCCCACGTCCGTaa                                                                      >  1:4125906/1‑93 (MQ=255)
                  aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTaa                                                                      >  2:3071124/1‑93 (MQ=255)
                   ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATGTCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAg                                                                     >  1:4434944/1‑93 (MQ=37)
                   ccAAGTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGGG‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg                                                                     <  2:5928417/93‑1 (MQ=37)
                    cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGGCTTGAACCCCCACGGCCCTAAg                                                                     >  1:2304863/1‑93 (MQ=38)
                    cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCCCGTCCGTAAGa                                                                    >  1:4292414/1‑93 (MQ=38)
                       ttCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAaca                                                                 <  1:1483157/93‑1 (MQ=37)
                       ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATACCGCCTCGTTTGGGGGGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAaca                                                                 <  2:5719590/93‑1 (MQ=37)
                          cGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTGGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTa                                                              <  1:5742140/93‑1 (MQ=25)
                           gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAac                                                            >  1:2903099/1‑93 (MQ=255)
                             cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAAGACTAacc                                                           >  1:1255363/1‑92 (MQ=37)
                             cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACCCTaaac                                                           >  2:3983699/1‑91 (MQ=37)
                             cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAacac                                                          >  2:618477/1‑93 (MQ=255)
                               ccTTCGCACATTCATCTTATTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                         <  2:4565708/93‑1 (MQ=25)
                               ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAAAACAATAACAcc                                                         >  2:1958508/1‑93 (MQ=37)
                                cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                        >  1:2817085/1‑93 (MQ=255)
                                 ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                       >  1:6255787/1‑93 (MQ=255)
                                 ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                       >  2:5781135/1‑93 (MQ=255)
                                  tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGGACCCCCACGTCCGTAAGAACAGTAACACCTGa                                                      >  1:5803056/1‑93 (MQ=25)
                                  tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGa                                                      <  1:4031510/93‑1 (MQ=255)
                                   cGCACAGTCTTCTTTTTTTTTTTGATATCGCTTCTTTGGGGGCGAGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGaa                                                     <  1:1203076/93‑1 (MQ=18)
                                   cGCACAGTCATCTTATTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                     <  2:4454077/93‑1 (MQ=25)
                                     cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGc                                                   >  2:1919608/1‑93 (MQ=255)
                                      acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGCCCGTAAAAACACTAAAACCTGAAGCt                                                  >  2:413668/1‑93 (MQ=25)
                                      acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                 <  2:5159062/93‑1 (MQ=255)
                                       caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                 <  2:3180994/93‑1 (MQ=255)
                                       caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                 >  1:3203491/1‑93 (MQ=255)
                                       caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                 >  2:3527407/1‑93 (MQ=255)
                                        aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                <  1:817040/93‑1 (MQ=255)
                                         tgtaTCTTACTTTTTTTAATATCGCCTGGTTGGGTGCGAGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                               <  1:3304109/90‑1 (MQ=18)
                                         gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                               <  1:1921583/93‑1 (MQ=255)
                                         gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                               <  2:4264013/93‑1 (MQ=21)
                                         gTCATCTTACTTTTTTCGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                               <  1:273636/93‑1 (MQ=21)
                                          tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAACCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                              >  1:2204169/1‑93 (MQ=21)
                                            aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                            >  2:5502572/1‑93 (MQ=255)
                                            aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGCCCGTAAGAAAACTAACACCTGAAGCTAGcgcg                                            >  2:5275043/1‑93 (MQ=21)
                                             tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGt                                           >  1:6177298/1‑93 (MQ=255)
                                              cTTACTTTTTTTTATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                         <  1:139000/93‑1 (MQ=255)
                                              cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGGG‑‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                          <  1:2609320/93‑1 (MQ=21)
                                               ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                         >  1:2284638/1‑93 (MQ=255)
                                               ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                         >  2:5951047/1‑93 (MQ=255)
                                               ttACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                        >  1:6213772/1‑93 (MQ=255)
                                               ttACTTTTTTATATATCGCCTCGTTTGGTGCGAGGGGGGGGGTCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                         <  1:4198125/93‑1 (MQ=17)
                                                 ttgttttttGATACCGCCTCGTTGGGGGGGAGGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                        <  2:432646/90‑1 (MQ=14)
                                                 aCTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                       >  2:1704337/1‑93 (MQ=255)
                                                   tttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                     >  1:3926931/1‑93 (MQ=255)
                                                     tttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                   >  1:473157/1‑93 (MQ=255)
                                                     tttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                   >  2:3149510/1‑93 (MQ=255)
                                                     tttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                   <  1:3774028/93‑1 (MQ=255)
                                                      ttttGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAAACCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                 >  2:2126853/1‑93 (MQ=255)
                                                        ttGATAGCTCCTCGTGTGGTGGGGC‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                <  1:2579740/93‑1 (MQ=14)
                                                         tGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTccc                               >  1:3075200/1‑92 (MQ=17)
                                                         tGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                               >  2:1208031/1‑93 (MQ=255)
                                                         tGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                               <  1:718356/93‑1 (MQ=255)
                                                         tGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                               >  2:3589717/1‑93 (MQ=255)
                                                          gATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                              <  2:3682494/93‑1 (MQ=255)
                                                          gATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                              >  1:512491/1‑93 (MQ=255)
                                                           atatCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACAGCTGAAGCTAGCGCGTCTACCAATTCCGcc                             >  1:3257421/1‑93 (MQ=17)
                                                           atatCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                             >  2:5462849/1‑93 (MQ=255)
                                                            aagcgccgcGTTTGGTGCGGG‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                            <  2:620456/90‑1 (MQ=14)
                                                             tgccccTCTTTGGGTGGGAGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                           <  2:3459229/89‑1 (MQ=11)
                                                             atCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                           >  1:3659702/1‑93 (MQ=255)
                                                              tCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCCCGTCTACCAATTCCGCCAcc                          >  2:821100/1‑93 (MQ=14)
                                                               cGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                        >  1:3737849/1‑93 (MQ=255)
                                                               cGCCTCGTTTGGGGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                         <  1:6093550/93‑1 (MQ=14)
                                                                gCCTCGTTTGGGGCGAGGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                         <  2:6283318/93‑1 (MQ=14)
                                                                 gcTCGTTTGGGGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                       <  1:1878041/92‑1 (MQ=14)
                                                                 ccTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                      <  2:5139746/93‑1 (MQ=255)
                                                                 ccTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                      >  1:4769310/1‑93 (MQ=255)
                                                                 ccTCGTTTGGCGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCCATTCCGCCACCTTc                       >  1:1074392/1‑93 (MQ=14)
                                                                  cTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                      >  1:5163771/1‑93 (MQ=255)
                                                                    cGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                    >  1:959614/1‑93 (MQ=255)
                                                                     gTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                   >  2:5559125/1‑93 (MQ=255)
                                                                     cttgggTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                  <  1:5410343/89‑1 (MQ=255)
                                                                      tttGGTGCGAGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                  <  2:4907933/93‑1 (MQ=11)
                                                                      tttGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTCCCAATTCCGCCACCTTCGCCCa                  >  2:2627145/1‑93 (MQ=11)
                                                                       ttGGTGCGGG‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                 <  2:1840166/93‑1 (MQ=11)
                                                                       ttGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                 <  2:3906980/93‑1 (MQ=255)
                                                                         ggTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCCCCCCGTg               >  1:4840528/1‑93 (MQ=11)
                                                                           tGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg             >  2:2962243/1‑93 (MQ=255)
                                                                           tGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg             <  2:748194/93‑1 (MQ=255)
                                                                                ggagggggggggCCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg           <  2:363390/89‑1 (MQ=14)
                                                                                         gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCaa  <  2:3701162/93‑1 (MQ=255)
                                                                                            |                                                                                          
CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAA  >  NC_000913/4603697‑4603874

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTG  >  NC_000913/4603692‑4603878
                                                                                                |                                                                                             
TAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGG                                                                                               <  SRR4302533.3123987/94‑1 (MQ=60)
TAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGG                                                                                               >  SRR4302533.3804432/1‑94 (MQ=60)
    ACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTG                                                                                           >  SRR4302533.3179031/1‑94 (MQ=60)
     CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGA                                                                                          >  SRR4302533.1565587/1‑94 (MQ=60)
        GAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAC                                                                                        >  SRR4302533.4180036/1‑94 (MQ=60)
        GAAGCGAGCGCGTCTACCAATCCCGCCCCCTTCGCACATTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAC                                                                                        <  SRR4302533.3496514/94‑1 (MQ=60)
         AAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACC                                                                                       >  SRR4302533.226871/1‑94 (MQ=60)
         AAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCC                                                                                      <  SRR4302533.3121078/94‑1 (MQ=60)
         AAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATAGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACC                                                                                       >  SRR4302533.1752157/1‑94 (MQ=60)
         AAGCTAGCGCGTCTACCAATTCAGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACC                                                                                       <  SRR4302533.5891245/94‑1 (MQ=60)
           GCTAGCGCGTCTACCAATTCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCC                                                                                    <  SRR4302533.1672320/94‑1 (MQ=60)
           GCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCC                                                                                     <  SRR4302533.2605961/94‑1 (MQ=60)
           GCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCC                                                                                     >  SRR4302533.4128296/1‑94 (MQ=60)
           GCTAGCGCGTCTACCAAATCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCC                                                                                     >  SRR4302533.1851731/1‑94 (MQ=60)
            CTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCC                                                                                    >  SRR4302533.1477028/1‑94 (MQ=60)
            CTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCC                                                                                    <  SRR4302533.4646437/94‑1 (MQ=60)
            CTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTCCGAGGGGGGGGGGCTTGAACCCCC                                                                                    <  SRR4302533.2728603/94‑1 (MQ=60)
               GCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACG                                                                                 >  SRR4302533.3832724/1‑94 (MQ=60)
               GCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACG                                                                                 <  SRR4302533.2314062/94‑1 (MQ=60)
                CGCGTTTACCAATTCCGCCCCCTTCGCCCAGTCATCTTCCTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTC                                                                               <  SRR4302533.3236131/94‑1 (MQ=60)
                CGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGT                                                                                >  SRR4302533.3572876/1‑94 (MQ=60)
                  CGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCC                                                                              >  SRR4302533.5921271/1‑94 (MQ=60)
                  CGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCC                                                                              >  SRR4302533.2716707/1‑94 (MQ=60)
                   GTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCG                                                                             >  SRR4302533.5923560/1‑94 (MQ=60)
                    TCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGCTTGAACCCCCACGTCCGTAA                                                                          >  SRR4302533.4150951/1‑94 (MQ=60)
                       ACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATGTCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAG                                                                         >  SRR4302533.4461753/1‑94 (MQ=60)
                        CCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGGCTTGAACCCCCACGGCCCTAAG                                                                         >  SRR4302533.2319140/1‑94 (MQ=60)
                        CCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCCCGTCCGTAAGA                                                                        >  SRR4302533.4318229/1‑94 (MQ=60)
                         CAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAA                                                                       <  SRR4302533.426033/94‑1 (MQ=60)
                         CAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAA                                                                       >  SRR4302533.6238274/1‑94 (MQ=60)
                          AATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAAC                                                                      >  SRR4302533.2486604/1‑94 (MQ=60)
                            TTCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACAC                                                                    <  SRR4302533.1492111/94‑1 (MQ=60)
                              CCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAA                                                                 >  SRR4302533.4900549/1‑94 (MQ=60)
                               CGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTGGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAA                                                                 <  SRR4302533.5776790/94‑1 (MQ=60)
                               CGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAAC                                                                >  SRR4302533.2920929/1‑94 (MQ=60)
                                ccctcctcCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAAC                                                                <  SRR4302533.4075527/86‑1 (MQ=60)
                                GCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAAC                                                                >  SRR4302533.4007973/1‑94 (MQ=60)
                                GCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAACCCCACGTCCGTAAGCACACTAAC                                                                >  SRR4302533.440599/1‑94 (MQ=60)
                                 CCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAAGACTAACC                                                               >  SRR4302533.1262868/1‑94 (MQ=60)
                                 CCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACA                                                               >  SRR4302533.5179003/1‑94 (MQ=60)
                                  CACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACAC                                                              <  SRR4302533.4143238/94‑1 (MQ=60)
                                    CCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCT                                                            >  SRR4302533.2834426/1‑94 (MQ=60)
                                     cctccctgtgcgtttttttttttttttggtatccgctttgtttggtgcgggGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACC                                                             <  SRR4302533.5561452/43‑1 (MQ=21)
                                     CTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTG                                                           >  SRR4302533.6293197/1‑94 (MQ=60)
                                      TTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGGACCCCCACGTCCGTAAGAACAGTAACACCTGA                                                          >  SRR4302533.5838107/1‑94 (MQ=60)
                                       TCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAA                                                         <  SRR4302533.4056050/94‑1 (MQ=60)
                                       TCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAG                                                        >  SRR4302533.741015/1‑94 (MQ=60)
                                        CGCACAGTCTTCTTTTTTTTTTTGATATCGCTTCTTTGGGGGCGAGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAG                                                        <  SRR4302533.1210307/94‑1 (MQ=36)
                                          caaactgtgatgttattttttttaatagcgccgtgttgggtgcgagGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAG                                                        <  SRR4302533.57941/48‑1 (MQ=22)
                                          CACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCT                                                      >  SRR4302533.6317795/1‑94 (MQ=60)
                                           ACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTA                                                     >  SRR4302533.3223207/1‑94 (MQ=60)
                                           ACAGTCATCTTACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTA                                                     <  SRR4302533.6252184/94‑1 (MQ=60)
                                            CAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAATCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAG                                                    >  SRR4302533.2869098/1‑94 (MQ=60)
                                             AGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGGCTTGAACCCCCACGCCCGTAAAAACACTAAAACCTGAAGCTAC                                                    >  SRR4302533.1751349/1‑94 (MQ=60)
                                             AGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGC                                                   <  SRR4302533.821827/94‑1 (MQ=60)
                                              tgtATCTTACTTTTTTTAATATCGCCTGGTTGGGTGCGAGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCG                                                  <  SRR4302533.3324513/91‑1 (MQ=37)
                                              GTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCG                                                  <  SRR4302533.1933312/94‑1 (MQ=60)
                                              GTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAAACCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCG                                                  >  SRR4302533.2217754/1‑94 (MQ=60)
                                              GTCATCTTACTTTTTTCGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCG                                                  <  SRR4302533.275300/94‑1 (MQ=60)
                                                 ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGT                                               >  SRR4302533.4769495/1‑94 (MQ=60)
                                                 ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGT                                               >  SRR4302533.6214409/1‑94 (MQ=60)
                                                 ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGT                                               >  SRR4302533.3742242/1‑94 (MQ=60)
                                                  tttttttttttttttattcccccttttttgggggggggGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCG                                                <  SRR4302533.20491/56‑1 (MQ=8)
                                                   CTTACTTTTTTTTATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTA                                            <  SRR4302533.139884/94‑1 (MQ=60)
                                                   CTTACTTTTTTTGATATCGCCTCGTTTGGTGCGGG‑‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCT                                             <  SRR4302533.2625340/94‑1 (MQ=58)
                                                   CTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCT                                             >  SRR4302533.2298777/1‑94 (MQ=60)
                                                   CTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTA                                            >  SRR4302533.6251031/1‑94 (MQ=60)
                                                    TTACTTTTTTATATATCGCCTCGTTTGGTGCGAGGGGGGGGGTCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTA                                            <  SRR4302533.4223535/94‑1 (MQ=43)
                                                     TACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAC                                           <  SRR4302533.1365887/94‑1 (MQ=60)
                                                       CTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCA                                         >  SRR4302533.3950904/1‑94 (MQ=59)
                                                         TTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAT                                       <  SRR4302533.218350/94‑1 (MQ=59)
                                                         TTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAT                                       >  SRR4302533.475986/1‑94 (MQ=54)
                                                          tttttttaaattcgcctgttttggggggggGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAA                                        <  SRR4302533.6323083/64‑1 (MQ=9)
                                                          TTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATT                                      <  SRR4302533.3797162/94‑1 (MQ=56)
                                                           TTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTC                                     >  SRR4302533.2943553/1‑94 (MQ=52)
                                                             TTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCC                                   >  SRR4302533.3094073/1‑94 (MQ=47)
                                                              TGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGC                                  >  SRR4302533.515536/1‑94 (MQ=45)
                                                              TGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGC                                  >  SRR4302533.5839285/1‑94 (MQ=45)
                                                              TGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGC                                  <  SRR4302533.722596/94‑1 (MQ=42)
                                                               ttgatagctcctcgtgtggtggggcGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCG                                   <  SRR4302533.2595615/69‑1 (MQ=8)
                                                               GATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACAGCTGAAGCTAGCGCGTCTACCAATTCCGCC                                 >  SRR4302533.3277510/1‑94 (MQ=38)
                                                               GATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCC                                 >  SRR4302533.5257842/1‑94 (MQ=42)
                                                                ATATCGCCTCGTTTGGTGCGAGGGGGGGGGGGGTTGAACCCCCCCGCCCGAAAGAACACTAAACCCTGAAGCTAGCGCGTCTACCAATTtcccc                                 >  SRR4302533.3592524/1‑89 (MQ=28)
                                                                ATATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCA                                >  SRR4302533.3279795/1‑94 (MQ=40)
                                                                 TATCGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAC                               >  SRR4302533.3682206/1‑94 (MQ=38)
                                                                   TCGCCTCGTTTGGTGCGAGGGGGGGGGGGTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCT                             >  SRR4302533.1254066/1‑94 (MQ=31)
                                                                   TCGCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTT                            >  SRR4302533.3760775/1‑94 (MQ=35)
                                                                    CTCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTT                            >  SRR4302533.99102/1‑94 (MQ=25)
                                                                    CGCCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTT                            >  SRR4302533.3180724/1‑94 (MQ=30)
                                                                    CGCCTCGTTTGGGGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTT                            <  SRR4302533.6130370/94‑1 (MQ=22)
                                                                     GCCTCGTTTGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCG                          >  SRR4302533.4798217/1‑94 (MQ=30)
                                                                     GCCTCGTTTGGCGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCCATTCCGCCACCTTC                           >  SRR4302533.1080725/1‑94 (MQ=19)
                                                                      GCTCGTTTGGGGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCG                          <  SRR4302533.1889487/94‑1 (MQ=10)
                                                                      CCTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCG                          >  SRR4302533.5194950/1‑94 (MQ=23)
                                                                       CTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCCCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCCCCACCTTCGC                         >  SRR4302533.578388/1‑94 (MQ=19)
                                                                       CTCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGC                         >  SRR4302533.4618036/1‑94 (MQ=21)
                                                                        gctgtttgtgtgggagGGGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCG                          <  SRR4302533.228447/78‑1 (MQ=7)
                                                                        TCGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCC                        >  SRR4302533.965259/1‑94 (MQ=21)
                                                                         CGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCC                       >  SRR4302533.3768326/1‑94 (MQ=21)
                                                                         CGTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGACCGTAAAGACACTAACACCAGAAGCTAGGGCGTCTACCAATTGCGCCACCTTCGCCC                       >  SRR4302533.5306146/1‑94 (MQ=57)
                                                                          cttgGGTGCGA‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAG                     <  SRR4302533.5442871/90‑1 (MQ=23)
                                                                          GTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCA                      >  SRR4302533.6152809/1‑94 (MQ=21)
                                                                           TTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAG                     <  SRR4302533.1932027/94‑1 (MQ=21)
                                                                           TTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAG                     >  SRR4302533.5197325/1‑94 (MQ=21)
                                                                             TGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTG                   <  SRR4302533.6127574/94‑1 (MQ=26)
                                                                             TGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCCCCCCGTG                   >  SRR4302533.4869867/1‑94 (MQ=5)
                                                                              ttggggggagGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGT                    <  SRR4302533.1268247/84‑1 (MQ=26)
                                                                               GTGCGG‑‑GTGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGA                <  SRR4302533.486799/94‑1 (MQ=31)
                                                                                TGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTTGCCCAGTGCGA                >  SRR4302533.2122093/1‑94 (MQ=27)
                                                                                 gggcgagGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCG                 <  SRR4302533.3640568/87‑1 (MQ=34)
                                                                                     G‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAAT           <  SRR4302533.3586366/94‑1 (MQ=47)
                                                                                     G‑‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACTGCGAGCAAT           >  SRR4302533.2812324/1‑94 (MQ=33)
                                                                                         GGGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATAT         <  SRR4302533.4726853/94‑1 (MQ=51)
                                                                                         GGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATAT         <  SRR4302533.873281/94‑1 (MQ=53)
                                                                                         GGGGGGGGCTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATAT         <  SRR4302533.342901/94‑1 (MQ=51)
                                                                                           GGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCA       <  SRR4302533.2668180/94‑1 (MQ=58)
                                                                                            GGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAA      >  SRR4302533.5709151/1‑94 (MQ=60)
                                                                                             GGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAC     <  SRR4302533.2416894/94‑1 (MQ=60)
                                                                                             GGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAC     >  SRR4302533.38712/1‑94 (MQ=60)
                                                                                              GGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACG    <  SRR4302533.2599919/94‑1 (MQ=60)
                                                                                              GGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACG    >  SRR4302533.5253444/1‑94 (MQ=60)
                                                                                              GGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACG    >  SRR4302533.3525723/1‑94 (MQ=60)
                                                                                               GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGT   >  SRR4302533.5713334/1‑94 (MQ=60)
                                                                                               GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGT   >  SRR4302533.2621594/1‑94 (MQ=60)
                                                                                                GCTTGACCCCCCACGTCCGTAAGAACACTACCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTG  <  SRR4302533.3409893/94‑1 (MQ=60)
                                                                                                GCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTG  <  SRR4302533.3417182/94‑1 (MQ=60)
                                                                                                |                                                                                             
TAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTG  >  NC_000913/4603692‑4603878

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