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,208,197 A→C T360P (ACC→CCC)  b2122 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,208,1970AC99.2% 371.8 / NA 119T360P (ACC→CCC) b2122hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (76/42);  minor base T (1/0);  total (77/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GGAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTACCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTG  >  NC_000913/2208107‑2208286
                                                                                          |                                                                                         
ggAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTccc                                                                                         >  2:1914282/1‑93 (MQ=255)
ggAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTccc                                                                                         >  2:2266975/1‑93 (MQ=255)
ggAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTccc                                                                                         <  1:1953866/93‑1 (MQ=255)
   gaAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGAt                                                                                      >  2:2291147/1‑93 (MQ=255)
   gaAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGAt                                                                                      >  2:4670454/1‑93 (MQ=255)
     aGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATcc                                                                                    >  2:53640/1‑93 (MQ=255)
      gCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATccc                                                                                   >  2:2829831/1‑93 (MQ=255)
       cAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCg                                                                                  >  2:2514790/1‑93 (MQ=255)
        aaaGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTtcccgatcccga                                                                                 >  1:303819/1‑93 (MQ=255)
        aaaGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTtcccgatcccga                                                                                 <  1:2375961/93‑1 (MQ=255)
         aaGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGac                                                                                >  1:1286749/1‑93 (MQ=255)
          aGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGaca                                                                               >  2:3267067/1‑93 (MQ=255)
          aGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGaca                                                                               >  2:4038233/1‑93 (MQ=255)
          aGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGaca                                                                               >  2:605146/1‑93 (MQ=255)
           gCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGacac                                                                              <  2:2792703/93‑1 (MQ=255)
            cTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACAcc                                                                             <  2:141874/93‑1 (MQ=255)
             ttGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCg                                                                            >  2:1524041/1‑93 (MQ=255)
             ttGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCg                                                                            >  1:3145099/1‑93 (MQ=255)
               gATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGc                                                                          >  2:3956928/1‑93 (MQ=255)
               gATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGc                                                                          >  2:922134/1‑93 (MQ=255)
                 tcatcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAg                                                                        >  1:703569/1‑93 (MQ=255)
                  catcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGt                                                                       <  2:1263225/93‑1 (MQ=255)
                  catcTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGt                                                                       >  1:1527786/1‑93 (MQ=255)
                     cTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAtt                                                                    <  2:3672231/93‑1 (MQ=255)
                     cTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAtt                                                                    <  2:4663291/93‑1 (MQ=255)
                     cTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAtt                                                                    <  1:3950217/93‑1 (MQ=255)
                      tGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAttt                                                                   >  1:4573446/1‑93 (MQ=255)
                      tGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAttt                                                                   >  2:86212/1‑93 (MQ=255)
                       gCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTAtttt                                                                  <  1:1327718/93‑1 (MQ=255)
                             ggTGTCTTTAGGTATGCGTTGCTAGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCt                                                            <  1:654490/93‑1 (MQ=255)
                              gtgtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCtt                                                           >  1:3575389/1‑93 (MQ=255)
                               tgtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCttt                                                          <  1:3029641/93‑1 (MQ=255)
                               tgtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCttt                                                          <  1:3869784/93‑1 (MQ=255)
                                gtCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATCGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTc                                                         >  2:1989154/1‑93 (MQ=255)
                                 tCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCTGATCCCGACACCGGCAGTATTTTGCACCTTTcg                                                        <  2:3392952/93‑1 (MQ=255)
                                  cTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcgc                                                       >  2:1052099/1‑93 (MQ=255)
                                   tttAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcgcg                                                      >  2:636751/1‑93 (MQ=255)
                                   tttAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTcgcg                                                      >  2:2901413/1‑93 (MQ=255)
                                     tAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCa                                                    <  1:635318/93‑1 (MQ=255)
                                       ggTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGt                                                  >  1:247949/1‑93 (MQ=255)
                                       ggTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGt                                                  <  2:3715690/93‑1 (MQ=255)
                                        gTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGtt                                                 >  1:2749196/1‑93 (MQ=255)
                                         tATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTg                                                >  1:3199707/1‑93 (MQ=255)
                                         tATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTg                                                >  1:2785375/1‑93 (MQ=255)
                                         tATGCGTTGCTGGCAGGATATTGAGCATTATGGGTTACGCATCTGGTTTTCCGATACCGACACCGGCAGTATTTTGCACCTTTCACGCAGTTg                                                >  2:1584919/1‑93 (MQ=255)
                                          aTGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGtttg                                               >  2:2636371/1‑91 (MQ=255)
                                          aTGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTgg                                               >  2:3369686/1‑93 (MQ=255)
                                           tGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGc                                              >  1:117805/1‑93 (MQ=255)
                                              gTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCcgc                                           >  1:4309552/1‑93 (MQ=255)
                                              gTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCcgc                                           >  1:3451720/1‑93 (MQ=255)
                                                tGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGa                                         <  2:745096/93‑1 (MQ=255)
                                                   tGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGt                                      >  1:2568027/1‑93 (MQ=255)
                                                     gCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGa                                    >  1:1620598/1‑93 (MQ=255)
                                                       aGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAAc                                  <  2:1511155/93‑1 (MQ=255)
                                                        ggATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACa                                 >  2:1288605/1‑93 (MQ=255)
                                                        ggATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACa                                 <  2:3868365/93‑1 (MQ=255)
                                                           tatTGAGCATTATGGTTTACGCATCTGGTTTCCCGGTCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCCAAGTGAACAGGa                              >  2:887384/1‑93 (MQ=255)
                                                           tatTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGa                              >  1:3653734/1‑93 (MQ=255)
                                                           tatTGAGCATTACGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGa                              >  1:313669/1‑93 (MQ=255)
                                                           tatTAAGCATTATAGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGa                              <  1:411579/93‑1 (MQ=255)
                                                            atTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGCTGGCCGCGAAGTGAACAGGaa                             >  2:244291/1‑93 (MQ=255)
                                                             ttGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGaaa                            >  1:1437608/1‑93 (MQ=255)
                                                             ttGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGaaa                            >  1:3849848/1‑93 (MQ=255)
                                                              tGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGaaaa                           <  2:257101/93‑1 (MQ=255)
                                                              tGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCCCGAGGGGAACAGGaaaa                           >  1:322778/1‑93 (MQ=255)
                                                                aGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACt                         <  1:3908638/93‑1 (MQ=255)
                                                                 gCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTc                        >  1:4021894/1‑93 (MQ=255)
                                                                  cATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCa                       >  2:3828115/1‑93 (MQ=255)
                                                                  cATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATCTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCa                       <  1:3039436/93‑1 (MQ=255)
                                                                   aTTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAc                      <  2:2456914/93‑1 (MQ=255)
                                                                   aTTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAc                      <  1:1032066/93‑1 (MQ=255)
                                                                   aTTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAc                      >  2:4531670/1‑93 (MQ=255)
                                                                    ttATGGTTTACGCATCTGGTTTCCCGATCCCGGCACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAcc                     >  2:1433190/1‑93 (MQ=255)
                                                                    ttATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAcc                     <  1:3619432/93‑1 (MQ=255)
                                                                    ttATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAcc                     >  2:4632171/1‑93 (MQ=255)
                                                                    ttATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAcc                     >  1:2878028/1‑93 (MQ=255)
                                                                     tATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGCAGTGCACAGGAAAACTCACCg                    >  2:1165857/1‑93 (MQ=255)
                                                                     tATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCg                    >  2:752230/1‑93 (MQ=255)
                                                                          tttACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGGACCTTTCGCGCAGCTGGCCGCGAAGTGAACAAGAAAACTCACCGGCAGc               >  1:824962/1‑93 (MQ=255)
                                                                          tttACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGc               >  2:2573565/1‑93 (MQ=255)
                                                                          tttACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGc               >  2:286/1‑93 (MQ=255)
                                                                            tACGCATCTGGTTTCCCGCTCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTa             <  1:2140632/93‑1 (MQ=255)
                                                                            tACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTa             >  2:3602996/1‑93 (MQ=255)
                                                                             aCGCATCTGGTTTCCCGAGCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAc            <  2:3094914/93‑1 (MQ=255)
                                                                             aCGCATCTAGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAAATCACCGGCAGCTAc            >  2:4560451/1‑93 (MQ=255)
                                                                               gCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAcgc          >  1:3347913/1‑93 (MQ=255)
                                                                                cATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAcgcg         >  2:2822836/1‑93 (MQ=255)
                                                                                cATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAcgcg         >  2:2006159/1‑93 (MQ=255)
                                                                                cATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAcgcg         >  2:3073029/1‑93 (MQ=255)
                                                                                 aTCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGt        <  1:2243683/93‑1 (MQ=255)
                                                                                 aTCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCACAGTTCGCCGCGAAGTGAACAGGAAAAATCAGCGGCAGCTAcgcgc        >  2:3157927/1‑92 (MQ=255)
                                                                                 aTCTGGCTTCCCGATCCCGACACCGGCAGGATTTTGTACCTTTCGCGCAGTTCCCCGCGTAGTGACCAGCAAAACTCACCGGCAGCTACGCGt        >  2:3038540/1‑93 (MQ=255)
                                                                                  tCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTc       <  1:2881056/93‑1 (MQ=255)
                                                                                  tCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTAGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAcgcggc       >  1:2510131/1‑91 (MQ=255)
                                                                                    tGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCgg     >  1:3446352/1‑93 (MQ=255)
                                                                                    tGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCgg     >  1:3339663/1‑93 (MQ=255)
                                                                                      gTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCt   <  1:2650309/93‑1 (MQ=255)
                                                                                       tttCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTg  >  1:97953/1‑93 (MQ=255)
                                                                                       tttCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTg  >  1:2848567/1‑93 (MQ=255)
                                                                                       tttCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTg  >  1:2921380/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTACCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTG  >  NC_000913/2208107‑2208286

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
GGAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTACCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTA  >  NC_000913/2208107‑2208290
                                                                                          |                                                                                             
GGAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCG                                                                                            <  SRR4302535.142560/94‑1 (MQ=60)
GGAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCG                                                                                            <  SRR4302535.1967928/94‑1 (MQ=60)
       CAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGA                                                                                     >  SRR4302535.305927/1‑94 (MQ=60)
        AAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGAC                                                                                    <  SRR4302535.2393030/94‑1 (MQ=60)
        AAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGAC                                                                                    >  SRR4302535.1295917/1‑94 (MQ=60)
           GCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTGCGCATCTGGTTTCCCGATCCCGACACC                                                                                 <  SRR4302535.644748/94‑1 (MQ=60)
            CTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCG                                                                                >  SRR4302535.3166971/1‑94 (MQ=60)
                ATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAG                                                                            >  SRR4302535.3935206/1‑94 (MQ=60)
                ATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAG                                                                            >  SRR4302535.708503/1‑94 (MQ=60)
                 TCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGT                                                                           >  SRR4302535.1538792/1‑94 (MQ=60)
                     CTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTT                                                                       <  SRR4302535.3977186/94‑1 (MQ=60)
                     CTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTT                                                                       >  SRR4302535.4604621/1‑94 (MQ=60)
                       GCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTG                                                                     <  SRR4302535.1337192/94‑1 (MQ=60)
                           CTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACC                                                                 >  SRR4302535.1784559/1‑94 (MQ=60)
                             GGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTT                                                               >  SRR4302535.3600176/1‑94 (MQ=60)
                             GGTGTCTTTAGGTATGCGTTGCTAGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTT                                                               <  SRR4302535.659060/94‑1 (MQ=60)
                               TGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTC                                                             <  SRR4302535.3896382/94‑1 (MQ=60)
                               TGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTC                                                             <  SRR4302535.3050849/94‑1 (MQ=60)
                                    TTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCA                                                        >  SRR4302535.3429128/1‑94 (MQ=60)
                                     TAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAG                                                       <  SRR4302535.639746/94‑1 (MQ=60)
                                      AGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGT                                                      >  SRR4302535.249652/1‑94 (MQ=60)
                                      AGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGT                                                      >  SRR4302535.3656623/1‑94 (MQ=60)
                                       GGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTT                                                     >  SRR4302535.2768795/1‑94 (MQ=60)
                                        GTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTG                                                    >  SRR4302535.3221859/1‑94 (MQ=60)
                                        GTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTG                                                    >  SRR4302535.2805195/1‑94 (MQ=60)
                                          ATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGC                                                  >  SRR4302535.118632/1‑94 (MQ=60)
                                           TGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCC                                                 >  SRR4302535.261008/1‑94 (MQ=60)
                                             CGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGC                                               >  SRR4302535.4339034/1‑94 (MQ=60)
                                             CGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGC                                               >  SRR4302535.3475712/1‑94 (MQ=60)
                                              GTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCG                                              <  SRR4302535.4115392/94‑1 (MQ=60)
                                                 GCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCTGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAG                                           <  SRR4302535.2078140/94‑1 (MQ=60)
                                                 GCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAG                                           >  SRR4302535.2920417/1‑94 (MQ=60)
                                                  CTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGT                                          >  SRR4302535.2586499/1‑94 (MQ=60)
                                                    GGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGA                                        >  SRR4302535.1632290/1‑94 (MQ=60)
                                                          ATATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGA                                  >  SRR4302535.3679031/1‑94 (MQ=60)
                                                          ATATTGAGCATTACGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGA                                  >  SRR4302535.315840/1‑94 (MQ=60)
                                                           TATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAA                                 >  SRR4302535.3228893/1‑94 (MQ=60)
                                                           TATTAAGCATTATAGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAA                                 <  SRR4302535.414469/94‑1 (MQ=60)
                                                            ATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAA                                >  SRR4302535.3876341/1‑94 (MQ=60)
                                                            ATTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAA                                >  SRR4302535.1447882/1‑94 (MQ=60)
                                                             TTGAGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCCCGAGGGGAACAGGAAAA                               >  SRR4302535.325010/1‑94 (MQ=60)
                                                                AGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTC                            <  SRR4302535.3935409/94‑1 (MQ=60)
                                                                AGCATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTC                            >  SRR4302535.4049414/1‑94 (MQ=60)
                                                                  CATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATCTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCAC                          <  SRR4302535.3060723/94‑1 (MQ=60)
                                                                   ATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACC                         <  SRR4302535.1039396/94‑1 (MQ=60)
                                                                   ATTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACC                         >  SRR4302535.2898409/1‑94 (MQ=60)
                                                                    TTATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCG                        <  SRR4302535.3644524/94‑1 (MQ=60)
                                                                     TATGGTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGG                       <  SRR4302535.846351/94‑1 (MQ=60)
                                                                         GTTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGGACCTTTCGCGCAGCTGGCCGCGAAGTGAACAAGAAAACTCACCGGCAGC                   >  SRR4302535.830745/1‑94 (MQ=60)
                                                                          TTTACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCT                  <  SRR4302535.2813693/94‑1 (MQ=60)
                                                                            TACGCATCTGGTTTCCCGCTCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAC                <  SRR4302535.2156042/94‑1 (MQ=60)
                                                                            TACGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTAC                <  SRR4302535.3101201/94‑1 (MQ=60)
                                                                              CGCATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGC              >  SRR4302535.3371234/1‑94 (MQ=60)
                                                                                CATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGT            >  SRR4302535.2380774/1‑94 (MQ=60)
                                                                                 ATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTC           <  SRR4302535.2259638/94‑1 (MQ=60)
                                                                                 ATCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTAGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGGC           >  SRR4302535.2528225/1‑94 (MQ=60)
                                                                                  TCTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCG          <  SRR4302535.2901448/94‑1 (MQ=60)
                                                                                   CTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGG         >  SRR4302535.3470319/1‑94 (MQ=60)
                                                                                   CTGGTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGG         >  SRR4302535.3362931/1‑94 (MQ=60)
                                                                                      GTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTG      >  SRR4302535.2941989/1‑94 (MQ=60)
                                                                                      GTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTG      >  SRR4302535.2868753/1‑94 (MQ=60)
                                                                                      GTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTG      <  SRR4302535.2669285/94‑1 (MQ=60)
                                                                                      GTTTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTG      >  SRR4302535.98646/1‑94 (MQ=60)
                                                                                        TTCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTT    >  SRR4302535.467108/1‑94 (MQ=60)
                                                                                         TCCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTT   <  SRR4302535.1872795/94‑1 (MQ=60)
                                                                                          CCCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTA  <  SRR4302535.2640425/94‑1 (MQ=60)
                                                                                          |                                                                                             
GGAGAAGCAAAGCTTGATCATCTGCGGCTGGTGTCTTTAGGTATGCGTTGCTGGCAGGATATTGAGCATTATGGTTTACGCATCTGGTTTACCGATCCCGACACCGGCAGTATTTTGCACCTTTCGCGCAGTTGGCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTA  >  NC_000913/2208107‑2208290

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