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 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,306,256 T→C E127G (GAA→GGA)  yojI ← putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,306,2560TC98.5% 414.9 / ‑4.0 130E127G (GAA→GGA) yojIputative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (70/58);  total (70/60)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 142 total aligned reads displayed.

ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGC  >  NC_000913/2306166‑2306345
                                                                                          |                                                                                         
aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGGCAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         >  2:742253/1‑93 (MQ=255)
aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         <  1:3823145/93‑1 (MQ=255)
aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         <  1:510617/93‑1 (MQ=255)
aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         <  2:1585215/93‑1 (MQ=255)
aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         >  2:2443646/1‑93 (MQ=255)
aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         >  1:4478159/1‑93 (MQ=255)
 tCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccca                                                                                        >  1:3438062/1‑92 (MQ=255)
  ccAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCgg                                                                                       >  2:4049028/1‑93 (MQ=255)
    aGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCa                                                                                     >  2:1734632/1‑93 (MQ=255)
     gATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAg                                                                                    >  2:3136189/1‑93 (MQ=255)
       tAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGAc                                                                                  >  2:1327031/1‑93 (MQ=255)
       tAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGAc                                                                                  >  2:2242533/1‑93 (MQ=255)
         gCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCATACGc                                                                                >  2:4111662/1‑93 (MQ=255)
            gTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGcaca                                                                             >  1:2834021/1‑93 (MQ=255)
                ccAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAg                                                                         >  2:4686361/1‑93 (MQ=255)
                 cAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGc                                                                        >  2:666542/1‑93 (MQ=255)
                 cAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGc                                                                        <  1:710177/93‑1 (MQ=255)
                 cAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCCGCAGACGCACAACAGc                                                                        >  2:666540/1‑93 (MQ=255)
                  aGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCa                                                                       <  2:2022937/93‑1 (MQ=255)
                   gCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCaa                                                                      >  1:2323941/1‑93 (MQ=255)
                    cAACATTTTGCCCTACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAAt                                                                     <  2:2483606/93‑1 (MQ=255)
                    cAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAAt                                                                     <  1:2493234/93‑1 (MQ=255)
                    cAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAAt                                                                     <  2:2938540/93‑1 (MQ=255)
                     aaCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATg                                                                    >  1:376676/1‑93 (MQ=255)
                     aaCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATg                                                                    >  1:983259/1‑93 (MQ=255)
                     aaCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATg                                                                    <  2:3796798/93‑1 (MQ=255)
                      aCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATgg                                                                   >  2:1030072/1‑93 (MQ=255)
                      aCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATgg                                                                   <  2:2308652/93‑1 (MQ=255)
                      aCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATgg                                                                   >  2:3860454/1‑93 (MQ=255)
                       cATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGt                                                                  <  1:859047/93‑1 (MQ=255)
                         ttttGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGa                                                                <  1:511195/93‑1 (MQ=255)
                               ccGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCTCCTGCACCAGTCCCGGCAGACGCACAACAGCAATCGTGATACTgc                                                          >  1:3128275/1‑93 (MQ=255)
                               ccGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGATGCACAAAAGCAATGGTGATATTgc                                                          <  2:2572149/93‑1 (MQ=255)
                               ccGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTgc                                                          <  1:3490992/93‑1 (MQ=255)
                                  aCAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGca                                                       <  1:440924/93‑1 (MQ=255)
                                    aGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGcaca                                                     <  1:1002133/93‑1 (MQ=255)
                                    aGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGcaca                                                     >  2:2895525/1‑93 (MQ=255)
                                    aGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGcaca                                                     >  1:3011347/1‑93 (MQ=255)
                                    aGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGcaca                                                     >  1:3011346/1‑93 (MQ=255)
                                      cATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATc                                                   <  1:1024135/93‑1 (MQ=255)
                                        tCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGc                                                 >  1:139543/1‑93 (MQ=255)
                                        tCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGc                                                 >  2:1086328/1‑93 (MQ=255)
                                         ccACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCt                                                >  1:4264160/1‑93 (MQ=255)
                                         ccACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCt                                                <  1:3357946/93‑1 (MQ=255)
                                         ccACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCt                                                <  1:3311579/93‑1 (MQ=255)
                                          cacaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTg                                               <  2:3492924/93‑1 (MQ=255)
                                          cacaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTg                                               >  1:3170134/1‑93 (MQ=255)
                                           acaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTgg                                              >  1:2126275/1‑93 (MQ=255)
                                           acaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTgg                                              <  1:281696/93‑1 (MQ=255)
                                            caGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGt                                             <  2:4397691/93‑1 (MQ=255)
                                            caGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGt                                             >  1:331821/1‑93 (MQ=255)
                                             aGATACGCCGCTTAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGtt                                            >  2:3911126/1‑93 (MQ=255)
                                             aGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGtt                                            >  1:1227071/1‑93 (MQ=255)
                                             aGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGtt                                            >  2:2550690/1‑93 (MQ=255)
                                             aGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGtt                                            <  1:2675013/93‑1 (MQ=255)
                                               aTACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTaa                                          >  2:1124067/1‑93 (MQ=255)
                                               aTACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTaa                                          >  1:2362174/1‑93 (MQ=255)
                                                tACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAc                                         >  2:3201610/1‑93 (MQ=255)
                                                tACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAc                                         >  1:1812411/1‑93 (MQ=255)
                                                 aCGCCGCTGAACCGATAGTGAGAATCATCCACTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAcc                                        <  2:710087/93‑1 (MQ=255)
                                                  cgccgcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAccc                                       >  1:3992177/1‑93 (MQ=255)
                                                     cgcTGAACCGATAGTGAGAATGATCCGCTGCTCCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGc                                    <  1:2992179/93‑1 (MQ=255)
                                                     cgcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGc                                    <  1:2992180/93‑1 (MQ=255)
                                                      gcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGcc                                   >  1:3450378/1‑93 (MQ=255)
                                                      gcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGcc                                   >  1:401848/1‑93 (MQ=255)
                                                      gcTGAACCGATAGTGAGAATGATACCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGcc                                   <  2:2597145/93‑1 (MQ=255)
                                                       cTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCa                                  >  2:2365987/1‑93 (MQ=255)
                                                           aCCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCaa                              <  2:450668/93‑1 (MQ=255)
                                                              gATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGa                           <  2:4406065/93‑1 (MQ=255)
                                                              gATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGa                           >  2:4497605/1‑93 (MQ=255)
                                                              gATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGa                           >  1:1347154/1‑93 (MQ=255)
                                                              gATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGa                           <  1:3457007/93‑1 (MQ=255)
                                                               aTAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAg                          <  2:873318/93‑1 (MQ=255)
                                                                tAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAgg                         >  2:154989/1‑93 (MQ=255)
                                                                 aGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgc                        <  1:3176688/93‑1 (MQ=255)
                                                                 aGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgc                        <  2:2989131/93‑1 (MQ=255)
                                                                  gTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgcg                       >  2:522699/1‑93 (MQ=255)
                                                                  gTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgcg                       >  1:23833/1‑93 (MQ=255)
                                                                  gTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgcg                       >  2:2042630/1‑93 (MQ=255)
                                                                    gagaATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCt                     <  1:2746105/93‑1 (MQ=255)
                                                                      gaATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATTGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTAc                   >  2:873924/1‑93 (MQ=255)
                                                                      gaATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCCGCAACGACGCGCTAc                   >  2:2940947/1‑93 (MQ=255)
                                                                           aTCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgt              <  1:1411794/93‑1 (MQ=255)
                                                                           aTCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgt              <  2:1500475/93‑1 (MQ=255)
                                                                            tCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgtt             >  1:1021371/1‑93 (MQ=255)
                                                                            tCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgtt             <  1:3995154/93‑1 (MQ=255)
                                                                            accccTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgtt             <  2:3716274/92‑1 (MQ=255)
                                                                             ccccTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgttg            >  1:161164/1‑93 (MQ=255)
                                                                              cccTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgttgt           <  2:1193446/93‑1 (MQ=255)
                                                                                 tGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCa        <  2:4152075/93‑1 (MQ=255)
                                                                                 cgCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTCCGCAGATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTCCCGAGTTGTTCa        <  2:1642192/92‑1 (MQ=255)
                                                                                  gCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCaa       >  2:395918/1‑93 (MQ=255)
                                                                                   cACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAAt      <  2:1426108/93‑1 (MQ=255)
                                                                                   cACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAAt      <  2:3799162/93‑1 (MQ=255)
                                                                                     ccAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgc    >  1:804042/1‑93 (MQ=255)
                                                                                     ccAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgc    >  1:2094678/1‑93 (MQ=255)
                                                                                      cAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgcg   <  1:4578618/93‑1 (MQ=255)
                                                                                       aGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgcgc  >  2:3246277/1‑93 (MQ=255)
                                                                                       aGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgcgc  <  2:836381/93‑1 (MQ=255)
                                                                                       aGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgcgc  <  1:1286053/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGC  >  NC_000913/2306166‑2306345

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGCT  >  NC_000913/2306165‑2306346
                                                                                           |                                                                                          
CATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCC                                                                                          >  SRR4302544.4523180/1‑94 (MQ=60)
CATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCC                                                                                          <  SRR4302544.2595708/94‑1 (MQ=60)
 CTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCG                                                                                         <  SRR4302544.3948456/94‑1 (MQ=60)
 ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCG                                                                                         <  SRR4302544.515656/94‑1 (MQ=60)
 ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCG                                                                                         <  SRR4302544.3861050/94‑1 (MQ=60)
 ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCG                                                                                         >  SRR4302544.3664071/1‑94 (MQ=60)
 ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCG                                                                                         >  SRR4302544.3038873/1‑94 (MQ=60)
        TAACCGTTACCAGCAACATTTTGCCCGACAGCCACCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACG                                                                                  <  SRR4302544.4718940/94‑1 (MQ=60)
          GCCGTTCCCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCA                                                                                <  SRR4302544.1466905/94‑1 (MQ=60)
            CGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACA                                                                              >  SRR4302544.3388830/1‑94 (MQ=60)
            CGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACA                                                                              >  SRR4302544.2861403/1‑94 (MQ=60)
             GTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAA                                                                             <  SRR4302544.133188/94‑1 (MQ=60)
                  CAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCA                                                                        <  SRR4302544.717220/94‑1 (MQ=60)
                  CAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCA                                                                        >  SRR4302544.4170561/1‑94 (MQ=60)
                   AGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAA                                                                       >  SRR4302544.2345973/1‑94 (MQ=60)
                    GCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAAT                                                                      <  SRR4302544.3050303/94‑1 (MQ=60)
                     CAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATG                                                                     >  SRR4302544.992898/1‑94 (MQ=60)
                     CAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATG                                                                     <  SRR4302544.2517117/94‑1 (MQ=60)
                     CAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATG                                                                     >  SRR4302544.380517/1‑94 (MQ=60)
                        CATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTG                                                                  <  SRR4302544.867545/94‑1 (MQ=60)
                          TTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGAT                                                                <  SRR4302544.516237/94‑1 (MQ=60)
                             TGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATT                                                             <  SRR4302544.2886911/94‑1 (MQ=60)
                               CCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCTCCTGCACCAGTCCCGGCAGACGCACAACAGCAATCGTGATACTGC                                                           >  SRR4302544.3158000/1‑94 (MQ=60)
                               CCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATAATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGCGATATTGC                                                           >  SRR4302544.2483881/1‑94 (MQ=60)
                                CCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCG                                                          <  SRR4302544.3524272/94‑1 (MQ=60)
                                 CGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGC                                                         >  SRR4302544.1890367/1‑94 (MQ=60)
                                   ACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCAC                                                       <  SRR4302544.445343/94‑1 (MQ=60)
                                    CAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACA                                                      >  SRR4302544.3040208/1‑94 (MQ=60)
                                    CAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACA                                                      >  SRR4302544.3040207/1‑94 (MQ=60)
                                     AGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACAT                                                     <  SRR4302544.1012077/94‑1 (MQ=60)
                                     AGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACAT                                                     <  SRR4302544.2277575/94‑1 (MQ=60)
                                       CATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCG                                                   <  SRR4302544.1034360/94‑1 (MQ=60)
                                        ATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGC                                                  >  SRR4302544.140866/1‑94 (MQ=60)
                                        ATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGC                                                  >  SRR4302544.3239648/1‑94 (MQ=60)
                                         TCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCT                                                 >  SRR4302544.4307661/1‑94 (MQ=60)
                                          CCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTG                                                <  SRR4302544.3389991/94‑1 (MQ=60)
                                          CCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTG                                                <  SRR4302544.3343116/94‑1 (MQ=60)
                                          CCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTG                                                >  SRR4302544.3200304/1‑94 (MQ=60)
                                           CACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGG                                               >  SRR4302544.2146663/1‑94 (MQ=60)
                                            ACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGT                                              >  SRR4302544.335271/1‑94 (MQ=60)
                                            ACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGT                                              <  SRR4302544.2771085/94‑1 (MQ=60)
                                            ACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGT                                              <  SRR4302544.284631/94‑1 (MQ=60)
                                            ACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGT                                              >  SRR4302544.730966/1‑94 (MQ=60)
                                             CAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTT                                             >  SRR4302544.1239210/1‑94 (MQ=60)
                                              AGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCNCGGCAGACGCACAAAAGCAATGGNNNNATNGCGCNCANCNCTNGNNN                                            >  SRR4302544.3434232/1‑94 (MQ=60)
                                              AGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTA                                            <  SRR4302544.2700656/94‑1 (MQ=60)
                                               GATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAA                                           >  SRR4302544.2384542/1‑94 (MQ=60)
                                                ATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAC                                          >  SRR4302544.1830108/1‑94 (MQ=60)
                                                  ACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCC                                        >  SRR4302544.4031917/1‑94 (MQ=60)
                                                      CGCTGAACCGATAGTGAGAATGATCCGCTGCTCCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCC                                    <  SRR4302544.3020873/94‑1 (MQ=60)
                                                      CGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCC                                    <  SRR4302544.3020874/94‑1 (MQ=60)
                                                      CGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCC                                    >  SRR4302544.405931/1‑94 (MQ=60)
                                                      CGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCC                                    >  SRR4302544.3483331/1‑94 (MQ=60)
                                                              CGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGA                            >  SRR4302544.1360500/1‑94 (MQ=60)
                                                              CGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGA                            >  SRR4302544.1213273/1‑94 (MQ=60)
                                                               GATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAG                           <  SRR4302544.3490002/94‑1 (MQ=60)
                                                                  AGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCG                        >  SRR4302544.24065/1‑94 (MQ=60)
                                                                  AGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCG                        <  SRR4302544.3206910/94‑1 (MQ=60)
                                                                    TGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCT                      <  SRR4302544.347408/94‑1 (MQ=60)
                                                                     GAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTA                     <  SRR4302544.2772560/94‑1 (MQ=60)
                                                                            ATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTT              <  SRR4302544.1425854/94‑1 (MQ=60)
                                                                            ATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTT              >  SRR4302544.1031570/1‑94 (MQ=60)
                                                                             TCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTG             <  SRR4302544.4034921/94‑1 (MQ=60)
                                                                             TCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTG             >  SRR4302544.162735/1‑94 (MQ=60)
                                                                                     ACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGC     >  SRR4302544.2114834/1‑94 (MQ=60)
                                                                                     ACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGC     >  SRR4302544.811875/1‑94 (MQ=60)
                                                                                      CCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCG    <  SRR4302544.949579/94‑1 (MQ=60)
                                                                                       CAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGC   <  SRR4302544.4624225/94‑1 (MQ=60)
                                                                                        AGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGCT  <  SRR4302544.4459241/94‑1 (MQ=60)
                                                                                        AGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGCT  <  SRR4302544.1298720/94‑1 (MQ=60)
                                                                                        AGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGCT  >  SRR4302544.2444996/1‑94 (MQ=60)
                                                                                           |                                                                                          
CATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGCT  >  NC_000913/2306165‑2306346

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