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 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,011,688 G→A P146P (CCG→CCA fliF → FliF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,011,6880GA96.4% 371.2 / ‑0.1 112P146P (CCG→CCAfliFFliF
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (47/61);  total (49/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 = 5.08e-01

Only 100 of 114 total aligned reads displayed.

AGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAA  >  NC_000913/2011600‑2011777
                                                                                        |                                                                                         
aGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTaa                                                                                       <  1:3363915/93‑1 (MQ=255)
 gTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTaaa                                                                                      >  2:4425374/1‑93 (MQ=255)
 gTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTaaa                                                                                      >  2:1167936/1‑93 (MQ=255)
  tttGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTaaaa                                                                                     <  1:2814986/93‑1 (MQ=255)
      gTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAgggg                                                                                 <  1:1193276/93‑1 (MQ=255)
      gTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAgggg                                                                                 <  1:1218937/93‑1 (MQ=255)
       tATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgc                                                                                >  1:290231/1‑93 (MQ=255)
       tATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgc                                                                                <  2:3386469/93‑1 (MQ=255)
        aTCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTAACATCGAAACTATCGGCCCAGTAAAAggggtg                                                                               >  2:2513079/1‑91 (MQ=255)
         tCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgcgc                                                                              >  1:1923051/1‑93 (MQ=255)
           aGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgcgcgc                                                                            >  2:1365053/1‑93 (MQ=255)
           aGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgcgcgc                                                                            >  1:1075623/1‑93 (MQ=255)
            gCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGcgcgcg                                                                           >  2:3652406/1‑93 (MQ=255)
             ccAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGt                                                                          <  2:3164153/93‑1 (MQ=255)
              cAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACGATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTa                                                                         <  2:1540279/93‑1 (MQ=255)
               aGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTAc                                                                        <  2:1688080/93‑1 (MQ=255)
                gTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACa                                                                       <  1:3218479/93‑1 (MQ=255)
                gTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACa                                                                       >  1:2015604/1‑93 (MQ=255)
                   cAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCt                                                                    >  2:1785478/1‑93 (MQ=255)
                   cAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCt                                                                    <  1:4598482/93‑1 (MQ=255)
                    aGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTg                                                                   <  2:3267637/93‑1 (MQ=255)
                     gCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTgg                                                                  <  2:1211991/93‑1 (MQ=255)
                      cGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGc                                                                 <  2:4183920/93‑1 (MQ=255)
                       gAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCa                                                                >  2:2271163/1‑93 (MQ=255)
                       gAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCa                                                                >  1:4015593/1‑93 (MQ=255)
                        aaCAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCaa                                                               >  2:10589/1‑93 (MQ=255)
                         aCAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAAt                                                              <  1:3101659/93‑1 (MQ=255)
                         aCAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAAt                                                              <  2:3888521/93‑1 (MQ=255)
                           aGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGc                                                            <  1:4421187/93‑1 (MQ=255)
                            ggTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGcc                                                           >  2:3541036/1‑93 (MQ=255)
                            ggTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGcc                                                           >  1:2475466/1‑93 (MQ=255)
                             gTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTGCCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCg                                                          <  1:294124/93‑1 (MQ=255)
                             gTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCg                                                          <  1:3719440/93‑1 (MQ=255)
                             gTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCg                                                          >  2:550492/1‑93 (MQ=255)
                                aaTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGaaa                                                       >  1:3422303/1‑93 (MQ=255)
                                aaTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGaaa                                                       >  1:1605872/1‑93 (MQ=255)
                                 aTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAAc                                                      <  1:882794/93‑1 (MQ=255)
                                 aTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAAc                                                      <  1:3547984/93‑1 (MQ=255)
                                  ttATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAAcc                                                     >  2:4398406/1‑93 (MQ=255)
                                   tATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCg                                                    <  2:61492/93‑1 (MQ=255)
                                    aTCAGGGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGt                                                   <  1:1681569/93‑1 (MQ=255)
                                    aTCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTATCATCGAAACTATCGGGCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGGAACCGt                                                   >  1:1361221/1‑93 (MQ=255)
                                    aTCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGt                                                   >  1:165418/1‑93 (MQ=255)
                                    aTCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGt                                                   >  2:741738/1‑93 (MQ=255)
                                     tCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTc                                                  <  2:530125/93‑1 (MQ=255)
                                      cAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCt                                                 <  1:2894775/93‑1 (MQ=255)
                                       aGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCtt                                                >  1:3244050/1‑93 (MQ=255)
                                          gggCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGTGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAt                                             <  2:3281153/93‑1 (MQ=255)
                                          gggCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAt                                             <  2:938887/93‑1 (MQ=255)
                                           ggCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAtt                                            >  1:259724/1‑93 (MQ=255)
                                           ggCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAtt                                            <  1:477521/93‑1 (MQ=255)
                                                tGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTcc                                       <  2:4476889/93‑1 (MQ=255)
                                                tGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTcc                                       <  2:2305975/93‑1 (MQ=255)
                                                tGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTcc                                       >  2:3577790/1‑93 (MQ=255)
                                                tGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTcc                                       <  1:180638/93‑1 (MQ=255)
                                                 ggAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCg                                      >  1:746646/1‑93 (MQ=255)
                                                 ggAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCg                                      >  2:4221644/1‑93 (MQ=255)
                                                    aGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGa                                   <  2:2617685/93‑1 (MQ=255)
                                                     ggCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGaa                                  <  2:2785281/93‑1 (MQ=255)
                                                     ggCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGaa                                  <  2:1426310/93‑1 (MQ=255)
                                                     ggCGAGCTTTCTCGTACCATCGAAACTATCGGACCAGTTAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGaa                                  >  1:391926/1‑93 (MQ=255)
                                                      gCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTTTGTGAAc                                 >  1:3726094/1‑93 (MQ=255)
                                                      gCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAAc                                 <  2:183956/93‑1 (MQ=255)
                                                        gAGCTTTCTCGTACCATCGAACCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCTTCTTTATTCGTCCGTGAACaa                               <  2:697823/93‑1 (MQ=255)
                                                        gAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaa                               <  1:3933348/93‑1 (MQ=255)
                                                        gAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaa                               <  1:2951265/93‑1 (MQ=255)
                                                         aGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaa                              >  1:64195/1‑93 (MQ=255)
                                                         aGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaa                              >  2:3619181/1‑93 (MQ=255)
                                                          gCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaa                             <  2:2894445/93‑1 (MQ=255)
                                                          gCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaa                             >  1:344586/1‑93 (MQ=255)
                                                           cTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaaa                            >  1:4412650/1‑93 (MQ=255)
                                                           cTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaaa                            <  1:1156457/93‑1 (MQ=255)
                                                           cTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaaa                            <  1:2155962/93‑1 (MQ=255)
                                                           cTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaaa                            >  1:2688709/1‑93 (MQ=255)
                                                              tcncGTACCATCGAAACTATCNNCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAAccc                         <  1:2302048/93‑1 (MQ=255)
                                                               ctcGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAAcccc                        <  1:1369124/93‑1 (MQ=255)
                                                                 cGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCt                      <  2:3424677/93‑1 (MQ=255)
                                                                 cGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCt                      <  1:3170664/93‑1 (MQ=255)
                                                                   tACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTc                    >  1:1106437/1‑93 (MQ=255)
                                                                     ccATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCtc                  >  1:1828666/1‑92 (MQ=255)
                                                                      cATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGc                 <  2:4466205/93‑1 (MQ=255)
                                                                       aTCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCa                >  1:818780/1‑93 (MQ=255)
                                                                        tCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCAt               <  2:4310177/93‑1 (MQ=255)
                                                                        tCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCAt               >  1:831938/1‑93 (MQ=255)
                                                                         cGAATCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATc              <  2:347246/93‑1 (MQ=255)
                                                                          gAAACTATCGGCCCGGTAAAAGGGGCGGGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCg             <  1:2244686/93‑1 (MQ=255)
                                                                          gAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCg             >  1:4300811/1‑93 (MQ=255)
                                                                           aaaCTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCgg            >  1:4207220/1‑93 (MQ=255)
                                                                            aaCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGt           >  2:3230393/1‑93 (MQ=255)
                                                                            aaCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAAGGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAGCCCCCTTCTGCATCGGt           <  1:2202901/93‑1 (MQ=255)
                                                                             aCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTg          <  1:744145/93‑1 (MQ=255)
                                                                              cTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGa         >  1:4043376/1‑93 (MQ=255)
                                                                              cTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGa         <  1:2455307/93‑1 (MQ=255)
                                                                               tATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGAc        >  2:2661414/1‑93 (MQ=255)
                                                                               tATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGAc        <  2:977233/93‑1 (MQ=255)
                                                                                aTCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACg       <  1:135750/93‑1 (MQ=255)
                                                                                aTCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACACCCTTCTGCATCGGTGACg       >  2:2871216/1‑93 (MQ=255)
                                                                                 tCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACgg      >  2:2447307/1‑93 (MQ=255)
                                                                                    gCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTaa   >  1:1592289/1‑93 (MQ=255)
                                                                                     cccAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTaaa  <  2:1955411/93‑1 (MQ=255)
                                                                                        |                                                                                         
AGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAA  >  NC_000913/2011600‑2011777

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATC  >  NC_000913/2011595‑2011779
                                                                                             |                                                                                           
GGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCA                                                                                             <  SRR4302535.1083245/94‑1 (MQ=60)
     AGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAA                                                                                        <  SRR4302535.3387323/94‑1 (MQ=60)
       TTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAG                                                                                      <  SRR4302535.2834986/94‑1 (MQ=60)
        TTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGG                                                                                     <  SRR4302535.19809/94‑1 (MQ=60)
           GTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGC                                                                                  <  SRR4302535.1201826/94‑1 (MQ=60)
           GTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGC                                                                                  <  SRR4302535.1227680/94‑1 (MQ=60)
           GTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGC                                                                                  >  SRR4302535.292265/1‑94 (MQ=60)
             ATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGC                                                                                >  SRR4302535.1936911/1‑94 (MQ=60)
               CAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGC                                                                              >  SRR4302535.1083270/1‑94 (MQ=60)
                    AGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACA                                                                         >  SRR4302535.2030127/1‑94 (MQ=60)
                     GTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACAT                                                                        <  SRR4302535.3240769/94‑1 (MQ=60)
                        CAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTG                                                                     <  SRR4302535.4629764/94‑1 (MQ=60)
                         AGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGG                                                                    >  SRR4302535.254039/1‑94 (MQ=60)
                           CGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCA                                                                  >  SRR4302535.4043059/1‑94 (MQ=60)
                           CGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCA                                                                  >  SRR4302535.3010198/1‑94 (MQ=60)
                              ACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATG                                                               <  SRR4302535.3123250/94‑1 (MQ=60)
                                AGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCC                                                             >  SRR4302535.2493320/1‑94 (MQ=60)
                                AGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCC                                                             <  SRR4302535.4451463/94‑1 (MQ=60)
                                  GTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTGCCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGA                                                           <  SRR4302535.296177/94‑1 (MQ=60)
                                  GTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGA                                                           <  SRR4302535.3745205/94‑1 (MQ=60)
                                    GAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAA                                                         >  SRR4302535.3446107/1‑94 (MQ=60)
                                    GAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAA                                                         >  SRR4302535.1617459/1‑94 (MQ=60)
                                    GAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAA                                                         >  SRR4302535.1477779/1‑94 (MQ=60)
                                      ATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACC                                                       <  SRR4302535.3572613/94‑1 (MQ=60)
                                      ATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACC                                                       <  SRR4302535.889045/94‑1 (MQ=60)
                                        TATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTATCATCGAAACTATCGGGCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGGAACCGT                                                     >  SRR4302535.1370957/1‑94 (MQ=60)
                                        TATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGT                                                     >  SRR4302535.166571/1‑94 (MQ=60)
                                         ATCAGGGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTC                                                    <  SRR4302535.1693747/94‑1 (MQ=60)
                                           CAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTT                                                  <  SRR4302535.7427/94‑1 (MQ=60)
                                           CAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTT                                                  <  SRR4302535.2915249/94‑1 (MQ=60)
                                           CAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTT                                                  >  SRR4302535.3266584/1‑94 (MQ=60)
                                              CGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAT                                               <  SRR4302535.3494824/94‑1 (MQ=60)
                                               GGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATT                                              >  SRR4302535.261522/1‑94 (MQ=60)
                                                GGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTC                                             <  SRR4302535.480871/94‑1 (MQ=60)
                                                     TGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCG                                        >  SRR4302535.751898/1‑94 (MQ=60)
                                                     TGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCG                                        <  SRR4302535.181893/94‑1 (MQ=60)
                                                         AGGCGAGCTTTCTCGTACCATCGAAACTATCGGACCAGTTAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAA                                    >  SRR4302535.394660/1‑94 (MQ=60)
                                                          GGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTTTGTGAAC                                   >  SRR4302535.3751908/1‑94 (MQ=60)
                                                             GAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAA                                <  SRR4302535.2971997/94‑1 (MQ=60)
                                                             GAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAA                                >  SRR4302535.64679/1‑94 (MQ=60)
                                                             GAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAA                                <  SRR4302535.3960192/94‑1 (MQ=60)
                                                              AGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAA                               >  SRR4302535.346963/1‑94 (MQ=60)
                                                               GCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAA                              >  SRR4302535.2707922/1‑94 (MQ=60)
                                                               GCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAA                              >  SRR4302535.4442850/1‑94 (MQ=60)
                                                                CTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAAC                             <  SRR4302535.2171453/94‑1 (MQ=60)
                                                                CTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAAC                             <  SRR4302535.1164718/94‑1 (MQ=60)
                                                                   TCNCGTACCATCGAAACTATCNNCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCC                          <  SRR4302535.2318520/94‑1 (MQ=60)
                                                                    CTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCC                         <  SRR4302535.1378911/94‑1 (MQ=60)
                                                                      CGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTT                       <  SRR4302535.3192620/94‑1 (MQ=60)
                                                                       GTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTC                      >  SRR4302535.1114341/1‑94 (MQ=60)
                                                                         ACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTG                    <  SRR4302535.409656/94‑1 (MQ=60)
                                                                         ACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTC                    >  SRR4302535.1841868/1‑94 (MQ=60)
                                                                           CATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCA                  >  SRR4302535.824526/1‑94 (MQ=60)
                                                                            ATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCAT                 >  SRR4302535.837768/1‑94 (MQ=60)
                                                                              CGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCG               >  SRR4302535.4330241/1‑94 (MQ=60)
                                                                               GAAACTATCGGCCCGGTAAAAGGGGCGGGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGG              <  SRR4302535.2260649/94‑1 (MQ=60)
                                                                               GAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGG              >  SRR4302535.4236013/1‑94 (MQ=60)
                                                                                 AACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAAGGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAGCCCCCTTCTGCATCGGTG            <  SRR4302535.2218600/94‑1 (MQ=60)
                                                                                  ACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGA           >  SRR4302535.4071056/1‑94 (MQ=60)
                                                                                  ACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGA           <  SRR4302535.749369/94‑1 (MQ=60)
                                                                                   CTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGAC          <  SRR4302535.2473009/94‑1 (MQ=60)
                                                                                     ATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGG        <  SRR4302535.136697/94‑1 (MQ=60)
                                                                                        GGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAA     >  SRR4302535.1603777/1‑94 (MQ=60)
                                                                                          CCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAAT   >  SRR4302535.3391593/1‑94 (MQ=60)
                                                                                          CCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAAT   >  SRR4302535.3042597/1‑94 (MQ=60)
                                                                                           CCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATC  <  SRR4302535.2609713/94‑1 (MQ=60)
                                                                                             |                                                                                           
GGAAAAGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATC  >  NC_000913/2011595‑2011779

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