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 F41 I0 R1 212 10.4 395026 92.3% 364608 148.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,741 G→C 78.7% A168A (GCG→GCC talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,7410GC78.7% 139.8 / 36.6 61A168A (GCG→GCCtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base G (6/7);  new base C (27/21);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-02

GATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAG  >  CP009273/8598‑8882
                                                                                                                                               |                                                                                                                                             
gATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGGGCGGAAGCg                                                                                                                                         <  2:9711/150‑1 (MQ=255)
     gggcatcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCgt                                                                                                                                    >  1:126244/8‑150 (MQ=255)
      ggcatcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCgtc                                                                                                                                   >  1:145644/7‑149 (MQ=255)
        catcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCtt                                                                                                                                 >  1:156430/5‑150 (MQ=255)
         atcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTc                                                                                                                                >  1:194250/4‑150 (MQ=255)
            aGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTg                                                                                                                             <  1:1544/150‑1 (MQ=255)
                 cGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGAtttc                                                                                                                        >  2:132185/1‑150 (MQ=255)
                  gATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCg                                                                                                                       <  2:59341/150‑1 (MQ=255)
                    tCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGcc                                                                                                                     <  1:29056/150‑1 (MQ=255)
                     cgcaTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCg                                                                                                                    >  2:115272/4‑150 (MQ=255)
                        aTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGttc                                                                                                                 <  1:10952/150‑2 (MQ=255)
                              aTCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTgg                                                                                                            <  1:88225/149‑1 (MQ=255)
                               tCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTGCGCTCAGGCTCGCGCTTGCACCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGcc                                                                                                          <  2:166432/150‑1 (MQ=255)
                                cAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCg                                                                                                         >  1:145374/1‑150 (MQ=255)
                                cAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGATGGCCg                                                                                                         >  2:179167/1‑150 (MQ=255)
                                    cTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTAtc                                                                                                     <  1:132185/150‑2 (MQ=255)
                                       gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCtc                                                                                                  <  1:37766/144‑2 (MQ=255)
                                           ccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACt                                                                                              >  1:140175/3‑150 (MQ=255)
                                            caCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTg                                                                                             <  1:42370/149‑1 (MQ=255)
                                              ccTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACt                                                                                              <  2:87062/147‑1 (MQ=255)
                                              ccTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGt                                                                                           <  1:169713/150‑1 (MQ=255)
                                               cTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTa                                                                                          >  2:101798/1‑150 (MQ=255)
                                                tGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTa                                                                                          >  1:88846/1‑149 (MQ=255)
                                                tGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTa                                                                                          >  2:114722/1‑149 (MQ=255)
                                                 ggCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTAc                                                                                         >  2:190167/1‑149 (MQ=255)
                                                  gCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACa                                                                                        <  1:4235/149‑1 (MQ=255)
                                                        tATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGaa                                                                                 <  2:15725/150‑1 (MQ=255)
                                                            cGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATAcc                                                                             >  1:4973/1‑150 (MQ=255)
                                                                 cgctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGac                                                                          >  2:169550/5‑147 (MQ=255)
                                                                      aaCAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaaga                                                                      >  2:5890/1‑147 (MQ=255)
                                                                          gCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTa                                                               <  2:19050/150‑1 (MQ=255)
                                                                              gAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGca                                                           >  2:80813/1‑149 (MQ=255)
                                                                              gAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGca                                                           <  1:134773/150‑2 (MQ=255)
                                                                               aaaaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGcac                                                          >  1:52724/1‑148 (MQ=255)
                                                                                aaaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCAcc                                                         >  2:191553/1‑150 (MQ=255)
                                                                                 aaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCg                                                        >  1:67635/1‑150 (MQ=255)
                                                                                   agaagGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCgg                                                       >  2:49892/1‑149 (MQ=255)
                                                                                       ggCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCag                                                    >  1:107487/1‑148 (MQ=255)
                                                                                         cATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCagaaga                                                <  1:191553/150‑1 (MQ=255)
                                                                                           ttaaCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATc                                              <  2:194250/148‑1 (MQ=255)
                                                                                           tCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATc                                              <  1:38708/150‑1 (MQ=255)
                                                                                               cTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCggg                                          >  2:192567/1‑150 (MQ=255)
                                                                                               cTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCggg                                          <  1:158248/150‑1 (MQ=255)
                                                                                                     ccTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGAGTGGTACAAAGCCAACTCCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGt                                    <  1:179167/150‑1 (MQ=255)
                                                                                                      cTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGtt                                   <  1:8436/150‑1 (MQ=255)
                                                                                                      cTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGtt                                   >  1:168965/1‑150 (MQ=255)
                                                                                                       tGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGttt                                  >  2:62865/1‑150 (MQ=255)
                                                                                                       tGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGttt                                  >  1:172228/1‑150 (MQ=255)
                                                                                                           ccTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGt                              <  1:101798/150‑1 (MQ=255)
                                                                                                            ctgctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGTGCGTGGTTTCCGtc                             <  2:145374/150‑2 (MQ=255)
                                                                                                            ctgctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGtc                             <  1:49892/150‑2 (MQ=255)
                                                                                                             tgctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGCTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGt                              >  1:149479/1‑148 (MQ=255)
                                                                                                                  ttCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGa                        >  2:160980/1‑149 (MQ=255)
                                                                                                                       cTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTa                  <  1:160980/150‑1 (MQ=255)
                                                                                                                        ttCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTAc                 >  1:62505/1‑150 (MQ=255)
                                                                                                                            cTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGt             <  1:142050/150‑1 (MQ=255)
                                                                                                                             tCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGta            >  1:142210/1‑150 (MQ=255)
                                                                                                                                  cTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGtactac        >  2:147552/1‑149 (MQ=255)
                                                                                                                                    cgcgCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACaag     >  1:124849/4‑149 (MQ=255)
                                                                                                                                      cgCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACaag     >  2:162050/2‑147 (MQ=255)
                                                                                                                                       gCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACaagcag  <  1:162850/150‑5 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAG  >  CP009273/8598‑8882

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTAT  >  CP009273/8603‑8891
                                                                                                                                          |                                                                                                                                                      
gggcatcAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTC                                                                                                                                            >  SRR12117902.128545/8‑151 (MQ=60)
 GGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCT                                                                                                                                           >  SRR12117902.148309/1‑151 (MQ=60)
   catcAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTC                                                                                                                                         >  SRR12117902.159284/5‑151 (MQ=60)
    ATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCC                                                                                                                                        >  SRR12117902.197769/1‑151 (MQ=60)
      CAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTG                                                                                                                                      <  SRR12117902.1567/151‑1 (MQ=60)
              ATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCC                                                                                                                              <  SRR12117902.29571/151‑1 (MQ=60)
                  CATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTC                                                                                                                          <  SRR12117902.11138/151‑1 (MQ=60)
                        GATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGG                                                                                                                     <  SRR12117902.89837/150‑1 (MQ=60)
                           CAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGT                                                                                                                 >  SRR12117902.148035/1‑151 (MQ=60)
                              ACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATC                                                                                                              <  SRR12117902.134602/151‑1 (MQ=60)
                                 GGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATcctc                                                                                                           <  SRR12117902.38448/151‑5 (MQ=60)
                                      CCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTG                                                                                                      >  SRR12117902.142734/1‑151 (MQ=60)
                                      CCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTG                                                                                                      <  SRR12117902.43150/151‑1 (MQ=60)
                                        ACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGT                                                                                                    <  SRR12117902.172812/151‑1 (MQ=60)
                                           TGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTAC                                                                                                  >  SRR12117902.90464/1‑150 (MQ=60)
                                            ggcagggcattcgcgccgctGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACA                                                                                                 <  SRR12117902.4293/130‑1 (MQ=60)
                                                  gcattcgcgccgctGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCA                                                                                           >  SRR12117902.48021/15‑150 (MQ=60)
                                                    attcgcgccgctGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCcaac                                                                                         >  SRR12117902.66536/13‑146 (MQ=60)
                                                       CGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCG                                                                                     >  SRR12117902.5041/1‑151 (MQ=60)
                                                                        GGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCA                                                                    <  SRR12117902.137230/151‑1 (MQ=60)
                                                                          AAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACC                                                                  >  SRR12117902.53683/1‑151 (MQ=60)
                                                                            AAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGG                                                                >  SRR12117902.68853/1‑151 (MQ=60)
                                                                                  GGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGA                                                            >  SRR12117902.109447/1‑149 (MQ=60)
                                                                                   GCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGA                                                         <  SRR12117902.195034/151‑1 (MQ=60)
                                                                                     ATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATC                                                       <  SRR12117902.39427/151‑1 (MQ=60)
                                                                                         ACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGG                                                   <  SRR12117902.161135/151‑1 (MQ=60)
                                                                                               ACCTGAGGCTGCTGTTCTCCTTCGCTCAGGCACGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGAGTGGTACAAAGCCAACTCCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGT                                             <  SRR12117902.182408/151‑1 (MQ=60)
                                                                                                CCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTT                                            <  SRR12117902.8591/151‑1 (MQ=60)
                                                                                                 CTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTT                                           >  SRR12117902.172055/1‑151 (MQ=60)
                                                                                                  TGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTC                                          >  SRR12117902.175353/1‑151 (MQ=60)
                                                                                                     CCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGT                                       <  SRR12117902.103654/151‑1 (MQ=60)
                                                                                                      GCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtc                                      <  SRR12117902.50811/151‑5 (MQ=60)
                                                                                                        TGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGCTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtc                                      >  SRR12117902.152197/1‑145 (MQ=60)
                                                                                                                 CCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTA                           <  SRR12117902.163914/151‑1 (MQ=60)
                                                                                                                   TTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacg                         >  SRR12117902.63642/1‑134 (MQ=60)
                                                                                                                      GCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagt                      <  SRR12117902.144632/151‑21 (MQ=60)
                                                                                                                        TCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtat                    >  SRR12117902.144796/1‑129 (MQ=60)
                                                                                                                               CGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtattacaagc             >  SRR12117902.127127/1‑122 (MQ=60)
                                                                                                                                        GCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtattacaagcagcacggct    <  SRR12117902.64004/151‑39 (MQ=60)
                                                                                                                                          CGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCcgtcagcgaaatctacgagtattacaagcagcacggctat  >  SRR12117902.151021/1‑111 (MQ=60)
                                                                                                                                          |                                                                                                                                                      
TGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTAT  >  CP009273/8603‑8891

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