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,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,2560TC99.2% 392.9 / ‑4.7 125E127G (GAA→GGA) yojIputative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (66/58);  total (66/59)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 132 total aligned reads displayed.

ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGC  >  NC_000913/2306166‑2306345
                                                                                          |                                                                                         
aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         >  1:4105124/1‑93 (MQ=255)
aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         >  1:4338425/1‑93 (MQ=255)
aTCCAGATACCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTccc                                                                                         <  2:753173/93‑1 (MQ=255)
 tCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCg                                                                                        >  1:582572/1‑93 (MQ=255)
   cAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGc                                                                                      >  1:2070732/1‑93 (MQ=255)
     gATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAg                                                                                    >  2:454299/1‑93 (MQ=255)
     gATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAg                                                                                    >  2:404948/1‑93 (MQ=255)
     gATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCACCTGCACCAGTCCCGGCAg                                                                                    >  2:3270866/1‑93 (MQ=255)
     aaTAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAg                                                                                    <  2:1898198/92‑1 (MQ=255)
      aTAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGa                                                                                   <  2:2112331/93‑1 (MQ=255)
      aTAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGa                                                                                   <  1:737199/93‑1 (MQ=255)
       tAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGAc                                                                                  <  2:3484839/93‑1 (MQ=255)
       tAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGAc                                                                                  >  1:2670740/1‑93 (MQ=255)
        aGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACg                                                                                 <  2:1136963/93‑1 (MQ=255)
         gCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGc                                                                                <  2:1295451/93‑1 (MQ=255)
           cGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGcac                                                                              <  1:3504332/93‑1 (MQ=255)
             ttACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACaa                                                                            >  2:3092852/1‑93 (MQ=255)
               aCCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACaaaa                                                                          <  2:2991574/93‑1 (MQ=255)
                ccAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAg                                                                         >  1:4384489/1‑93 (MQ=255)
                ccAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAg                                                                         >  1:2397711/1‑93 (MQ=255)
                ccAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAg                                                                         >  2:2768456/1‑93 (MQ=255)
                ccAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAg                                                                         <  1:3561103/93‑1 (MQ=255)
                 cAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGc                                                                        <  2:2785110/93‑1 (MQ=255)
                  aGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCa                                                                       >  1:4304750/1‑93 (MQ=255)
                   ggaaCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGACCCCCTGCTCCAGTCCCGGCAGACGCACAAAAGCaa                                                                      <  2:1629955/91‑1 (MQ=255)
                   gCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCaa                                                                      <  2:2565879/93‑1 (MQ=255)
                     aaCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATg                                                                    >  2:2802173/1‑93 (MQ=255)
                     aaCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATg                                                                    >  2:3101500/1‑93 (MQ=255)
                     aaCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATg                                                                    <  2:2565593/93‑1 (MQ=255)
                     aaCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATg                                                                    >  2:4228041/1‑93 (MQ=255)
                      aCATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATgg                                                                   >  1:3415105/1‑93 (MQ=255)
                         ttttGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGa                                                                <  1:2913073/93‑1 (MQ=255)
                         ttttGCCCAACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGa                                                                <  2:1646784/93‑1 (MQ=255)
                          tttGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGat                                                               >  2:1489668/1‑93 (MQ=255)
                            tGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGatat                                                             <  1:4273198/93‑1 (MQ=255)
                            tGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGatat                                                             >  2:3330452/1‑93 (MQ=255)
                            tGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGatat                                                             >  1:1640257/1‑93 (MQ=255)
                             gCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATAtt                                                            >  2:3889455/1‑93 (MQ=255)
                               ccGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTgc                                                          >  1:1340138/1‑93 (MQ=255)
                                cGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTgcg                                                         >  1:2929505/1‑93 (MQ=255)
                                   cAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGcac                                                      >  2:3263166/1‑93 (MQ=255)
                                    aGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGcaca                                                     >  2:1469025/1‑93 (MQ=255)
                                       aTCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCg                                                  >  2:1202474/1‑93 (MQ=255)
                                        tCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGc                                                 <  1:3106887/93‑1 (MQ=255)
                                        tCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGc                                                 <  1:2603912/93‑1 (MQ=255)
                                        tCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGc                                                 >  2:554188/1‑93 (MQ=255)
                                         ccACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCt                                                >  1:2225828/1‑93 (MQ=255)
                                          cacaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTg                                               >  1:2128066/1‑93 (MQ=255)
                                          cacaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTg                                               >  1:1733614/1‑93 (MQ=255)
                                           acaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTgg                                              <  1:1311466/93‑1 (MQ=255)
                                           acaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTgg                                              >  1:4652784/1‑93 (MQ=255)
                                           acaGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTgg                                              >  1:118866/1‑93 (MQ=255)
                                            caGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGt                                             >  1:4483971/1‑93 (MQ=255)
                                            caGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGt                                             >  1:3450566/1‑93 (MQ=255)
                                             aGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAATGCAATGGTGATATTGCGCACATCGCTGGtt                                            <  2:2784842/93‑1 (MQ=255)
                                              gATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTa                                           <  1:1785345/93‑1 (MQ=255)
                                                tACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAc                                         >  1:1655560/1‑93 (MQ=255)
                                                tACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAc                                         >  1:1368721/1‑93 (MQ=255)
                                                  cgccgcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAccc                                       <  2:741567/93‑1 (MQ=255)
                                                  cgccgcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAccc                                       <  1:3533947/93‑1 (MQ=255)
                                                  cgccgcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAccc                                       <  2:2973171/93‑1 (MQ=255)
                                                  cgccgcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAccc                                       <  1:339156/93‑1 (MQ=255)
                                                  cgccgcTGAACCGATAGTGAGAATGACCCCCTGCACCAGCCCCGGCAGACGCACAAAAGCAATGGTGATATTGCTCACATCGCTGGTTAAccc                                       <  2:1802344/93‑1 (MQ=255)
                                                      gcTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGcc                                   <  2:786534/93‑1 (MQ=255)
                                                        tGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAg                                 >  1:233022/1‑93 (MQ=255)
                                                         gAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGc                                >  2:385715/1‑93 (MQ=255)
                                                           aCCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCaa                              <  1:2182846/93‑1 (MQ=255)
                                                           aCCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCaa                              >  2:3968383/1‑93 (MQ=255)
                                                           aCCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCaa                              >  2:3518587/1‑93 (MQ=255)
                                                            ccGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAAc                             >  1:4205155/1‑93 (MQ=255)
                                                              gATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGa                           <  2:412095/93‑1 (MQ=255)
                                                              gATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGa                           >  1:3777532/1‑93 (MQ=255)
                                                              gATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGACAGCAACGa                           >  1:4243693/1‑93 (MQ=255)
                                                               aTAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAg                          >  1:1565069/1‑93 (MQ=255)
                                                               aTAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAg                          <  2:4229838/93‑1 (MQ=255)
                                                               aTAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAg                          >  1:4335965/1‑93 (MQ=255)
                                                               aTAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAg                          >  2:3398587/1‑93 (MQ=255)
                                                                tAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAgg                         >  1:1004968/1‑93 (MQ=255)
                                                                tAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAgg                         >  2:617663/1‑93 (MQ=255)
                                                                tAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAgg                         <  1:2997258/93‑1 (MQ=255)
                                                                tAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAgg                         <  2:3415098/93‑1 (MQ=255)
                                                                  gTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgcg                       >  2:4154015/1‑93 (MQ=255)
                                                                   tGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgcgc                      >  1:1155927/1‑93 (MQ=255)
                                                                   tGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgcgc                      <  1:4015838/93‑1 (MQ=255)
                                                                   tGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGgcgc                      <  2:1480514/93‑1 (MQ=255)
                                                                       aaTGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTAcc                  >  1:2533119/1‑93 (MQ=255)
                                                                       aaTGATCCCCTCCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTAcc                  <  2:310607/93‑1 (MQ=255)
                                                                         tGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGa                <  1:4150243/93‑1 (MQ=255)
                                                                          gATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAg               >  1:1047210/1‑93 (MQ=255)
                                                                           aTCCCCTGCACCAGTCCCGGCAGACGCACGAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgt              <  2:4601734/93‑1 (MQ=255)
                                                                            tCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgtt             <  1:440727/93‑1 (MQ=255)
                                                                               ccTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAgttgtt          <  1:4133404/93‑1 (MQ=255)
                                                                                 tGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCa        >  2:2649013/1‑93 (MQ=255)
                                                                                  gCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCaa       <  1:3980812/93‑1 (MQ=255)
                                                                                   cACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAAt      <  1:1867487/93‑1 (MQ=255)
                                                                                    aCCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATg     <  2:1349729/93‑1 (MQ=255)
                                                                                     ccAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgc    <  2:3731299/93‑1 (MQ=255)
                                                                                      cAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgcg   >  2:565090/1‑93 (MQ=255)
                                                                                       aGTCCCGGCAGACGCCCAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgcgc  <  2:4333018/93‑1 (MQ=255)
                                                                                       aGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATgcgc  <  2:2103418/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGC  >  NC_000913/2306166‑2306345

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGCT  >  NC_000913/2306164‑2306346
                                                                                            |                                                                                          
CCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCC                                                                                           >  SRR4302535.3973560/1‑94 (MQ=60)
 CATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCC                                                                                          >  SRR4302535.4368093/1‑94 (MQ=60)
 CATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCC                                                                                          >  SRR4302535.4133229/1‑94 (MQ=60)
  ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCG                                                                                         >  SRR4302535.586657/1‑94 (MQ=60)
  ATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCG                                                                                         >  SRR4302535.4640535/1‑94 (MQ=60)
    CCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGC                                                                                       >  SRR4302535.2085656/1‑94 (MQ=60)
        ATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGAC                                                                                   <  SRR4302535.742383/94‑1 (MQ=60)
        ATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGAC                                                                                   >  SRR4302535.2689849/1‑94 (MQ=60)
             CGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACA                                                                              <  SRR4302535.3528683/94‑1 (MQ=60)
               TTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAA                                                                            >  SRR4302535.2201210/1‑94 (MQ=60)
                 ACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAG                                                                          >  SRR4302535.2414950/1‑94 (MQ=60)
                 ACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAG                                                                          >  SRR4302535.4414487/1‑94 (MQ=60)
                 ACCAACAACATGTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAG                                                                          <  SRR4302535.2032070/94‑1 (MQ=60)
                  CCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGC                                                                         >  SRR4302535.3608455/1‑94 (MQ=60)
                  CCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGC                                                                         <  SRR4302535.3585795/94‑1 (MQ=60)
                  CCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGC                                                                         >  SRR4302535.3519374/1‑94 (MQ=60)
                   CAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCA                                                                        >  SRR4302535.4334202/1‑94 (MQ=60)
                       AACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGG                                                                    >  SRR4302535.2775638/1‑94 (MQ=60)
                       AACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGG                                                                    >  SRR4302535.3438862/1‑94 (MQ=60)
                           TTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGAT                                                                <  SRR4302535.2933652/94‑1 (MQ=60)
                            TTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATA                                                               <  SRR4302535.3045038/94‑1 (MQ=60)
                             TTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATAT                                                              >  SRR4302535.1652116/1‑94 (MQ=60)
                             TTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATAT                                                              <  SRR4302535.1269604/94‑1 (MQ=60)
                              TGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATT                                                             <  SRR4302535.4302448/94‑1 (MQ=60)
                                CCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGC                                                           >  SRR4302535.1349700/1‑94 (MQ=60)
                                 CCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCG                                                          >  SRR4302535.2950144/1‑94 (MQ=60)
                                     CANNNTNCACAGATACGCCGCTNAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACA                                                      <  SRR4302535.1820372/94‑1 (MQ=60)
                                         ATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGC                                                  >  SRR4302535.4216467/1‑94 (MQ=60)
                                          TCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCT                                                 <  SRR4302535.3128517/94‑1 (MQ=60)
                                          TCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCT                                                 <  SRR4302535.4326625/94‑1 (MQ=60)
                                          TCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCT                                                 <  SRR4302535.2622605/94‑1 (MQ=60)
                                          TCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCT                                                 >  SRR4302535.2241663/1‑94 (MQ=60)
                                           CCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTG                                                >  SRR4302535.1746144/1‑94 (MQ=60)
                                           CCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTG                                                >  SRR4302535.2143412/1‑94 (MQ=60)
                                            CACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGG                                               >  SRR4302535.119698/1‑94 (MQ=60)
                                            CACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGG                                               >  SRR4302535.1728634/1‑94 (MQ=60)
                                            CACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGG                                               >  SRR4302535.4684287/1‑94 (MQ=60)
                                             ACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGT                                              >  SRR4302535.4514751/1‑94 (MQ=60)
                                             ACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGT                                              <  SRR4302535.1320838/94‑1 (MQ=60)
                                             ACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGT                                              >  SRR4302535.3474550/1‑94 (MQ=60)
                                                GATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAA                                           <  SRR4302535.1798265/94‑1 (MQ=60)
                                                 ATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAC                                          >  SRR4302535.1667534/1‑94 (MQ=60)
                                                 ATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAC                                          >  SRR4302535.1380116/1‑94 (MQ=60)
                                                 ATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAAC                                          >  SRR4302535.1378504/1‑94 (MQ=60)
                                                    CGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCC                                       <  SRR4302535.341506/94‑1 (MQ=60)
                                                    CGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCC                                       <  SRR4302535.3558481/94‑1 (MQ=60)
                                                     GCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCTCACATCGCTGGTTAACCCCG                                      >  SRR4302535.2494141/1‑94 (MQ=60)
                                                         CTGAACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAG                                  >  SRR4302535.234632/1‑94 (MQ=60)
                                                             ACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAAC                              <  SRR4302535.2198411/94‑1 (MQ=60)
                                                             ACCGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAAC                              >  SRR4302535.4233934/1‑94 (MQ=60)
                                                               CGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGA                            >  SRR4302535.1963241/1‑94 (MQ=60)
                                                               CGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGA                            >  SRR4302535.3803571/1‑94 (MQ=60)
                                                               CGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGA                            >  SRR4302535.4099971/1‑94 (MQ=60)
                                                               CGATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGACAGCAACGA                            >  SRR4302535.4272739/1‑94 (MQ=60)
                                                                GATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAG                           >  SRR4302535.1576369/1‑94 (MQ=60)
                                                                GATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAG                           >  SRR4302535.4365621/1‑94 (MQ=60)
                                                                 ATAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGG                          >  SRR4302535.1012110/1‑94 (MQ=60)
                                                                  TAGTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGC                         <  SRR4302535.3018276/94‑1 (MQ=60)
                                                                    GTGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGC                       >  SRR4302535.1164186/1‑94 (MQ=60)
                                                                     TGAGAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCT                      <  SRR4302535.4043306/94‑1 (MQ=60)
                                                                        GAATGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACC                   >  SRR4302535.2551369/1‑94 (MQ=60)
                                                                           TGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAG                <  SRR4302535.4178663/94‑1 (MQ=60)
                                                                           TGATCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAG                >  SRR4302535.1054635/1‑94 (MQ=60)
                                                                              TCCCCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTG             <  SRR4302535.443812/94‑1 (MQ=60)
                                                                                 CCTGCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTC          <  SRR4302535.4161698/94‑1 (MQ=60)
                                                                                    GCACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAAT       <  SRR4302535.4008007/94‑1 (MQ=60)
                                                                                     CACCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATG      <  SRR4302535.1880947/94‑1 (MQ=60)
                                                                                      TCCAGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGC     <  SRR4302535.4405800/94‑1 (MQ=60)
                                                                                         AGTCCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATGGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGCT  >  SRR4302535.2533545/1‑94 (MQ=60)
                                                                                            |                                                                                          
CCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACATCGCTGGTTAACCCCGCCAGCAACGAGGCGCTACCGAGTTGTTCAATGCGCT  >  NC_000913/2306164‑2306346

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