Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,957,451 T→C L140L (CTA→CTG b1875 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,957,4510TC99.2% 388.5 / ‑4.5 121L140L (CTA→CTGb1875hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (52/68);  total (52/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 126 total aligned reads displayed.

GTTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGTAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCA  >  NC_000913/1957362‑1957540
                                                                                         |                                                                                         
gttgttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAgg                                                                                        <  2:2967759/93‑1 (MQ=255)
 ttgttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGc                                                                                       <  2:64930/93‑1 (MQ=255)
 ttgttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGc                                                                                       >  2:39195/1‑93 (MQ=255)
  tgttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCt                                                                                      <  1:2104143/93‑1 (MQ=255)
   gttTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTg                                                                                     <  1:4684128/93‑1 (MQ=255)
       tGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTc                                                                                 >  1:1034072/1‑93 (MQ=255)
          cATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCtt                                                                              <  2:4562444/93‑1 (MQ=255)
            tCGGTGACCGCCAGCGTCGGATTCGTCAATCGCTCATACTCACCTTTTGGCGAACTGGTTTTCACTGAGGTCCCCGGCAGGCTGCGTCCTtca                                                                            >  1:3024212/1‑93 (MQ=255)
             cGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTtcat                                                                           <  2:1734212/93‑1 (MQ=255)
              ggTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTtcatc                                                                          >  1:4246701/1‑93 (MQ=255)
               gTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTtcatca                                                                         <  1:3095458/93‑1 (MQ=255)
                ggACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAg                                                                        <  1:1015916/92‑1 (MQ=255)
                   ccGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGCGTCCCTCATCaaaac                                                                     >  1:4364912/1‑92 (MQ=255)
                   ccGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGaaa                                                                     >  2:4719213/1‑93 (MQ=255)
                    cGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAg                                                                    <  1:614819/93‑1 (MQ=255)
                     gCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGc                                                                   >  1:1664453/1‑93 (MQ=255)
                     gCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGc                                                                   <  2:1748843/93‑1 (MQ=255)
                      ccAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGATCTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCa                                                                  <  2:4225515/93‑1 (MQ=255)
                      ccAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCa                                                                  <  1:3237533/93‑1 (MQ=255)
                      ccAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCa                                                                  >  1:4293623/1‑93 (MQ=255)
                         gCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATg                                                               >  1:490141/1‑93 (MQ=255)
                         gCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATg                                                               <  1:2153779/93‑1 (MQ=255)
                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGTAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             <  2:4528385/93‑1 (MQ=255)
                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             >  1:4628388/1‑93 (MQ=255)
                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             >  2:1301316/1‑93 (MQ=255)
                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             >  2:1443200/1‑93 (MQ=255)
                           gTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGcc                                                             >  1:429442/1‑93 (MQ=255)
                            tCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCa                                                            >  1:79341/1‑93 (MQ=255)
                             cGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAg                                                           >  1:133143/1‑93 (MQ=255)
                             cGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAg                                                           <  1:3769985/93‑1 (MQ=255)
                              ggATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAgc                                                          <  2:2414564/93‑1 (MQ=255)
                              ggATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAgc                                                          <  1:2900169/93‑1 (MQ=255)
                               gATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAgcg                                                         <  2:4292087/93‑1 (MQ=255)
                               gATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAgcg                                                         >  2:308896/1‑93 (MQ=255)
                                 ttCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCa                                                       <  2:3833406/93‑1 (MQ=255)
                                 ttCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGGTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCa                                                       >  1:837982/1‑93 (MQ=255)
                                   cGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACg                                                     <  2:1605852/93‑1 (MQ=255)
                                   cGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACg                                                     <  2:1992544/93‑1 (MQ=255)
                                    tgTAATCGCTCATGCTCGCCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACgg                                                    <  1:4503829/92‑1 (MQ=255)
                                      tAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGc                                                  >  1:1706820/1‑93 (MQ=255)
                                       aaTCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCa                                                 >  1:3754189/1‑93 (MQ=255)
                                        aTCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAt                                                >  2:217703/1‑93 (MQ=255)
                                        aTCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAt                                                >  1:3755024/1‑93 (MQ=255)
                                        aTCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAt                                                <  2:2605299/93‑1 (MQ=255)
                                         tCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAtt                                               >  2:3518642/1‑93 (MQ=255)
                                         tCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAtt                                               <  1:4539629/93‑1 (MQ=255)
                                         tCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAtt                                               <  1:4456395/93‑1 (MQ=255)
                                          cGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTc                                              <  2:431289/93‑1 (MQ=255)
                                          cGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTc                                              <  1:1028051/93‑1 (MQ=255)
                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:674478/1‑93 (MQ=255)
                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:2876774/1‑93 (MQ=255)
                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:527187/1‑93 (MQ=255)
                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:4284682/1‑93 (MQ=255)
                                            cTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAg                                            >  2:1851113/1‑93 (MQ=255)
                                             tCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           >  1:95069/1‑93 (MQ=255)
                                             tCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           >  2:2799564/1‑93 (MQ=255)
                                             tCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           <  1:2305083/93‑1 (MQ=255)
                                             tCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           <  2:1404076/93‑1 (MQ=255)
                                             tCATGCTCACCTTATGGCGCACTGATTTTCACTGAGATCCCAGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCaaa                                           <  1:4204240/93‑3 (MQ=255)
                                             ccATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGc                                           <  2:258029/92‑1 (MQ=255)
                                              cATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCg                                          >  2:4245752/1‑93 (MQ=255)
                                               aTGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGt                                         <  1:4474688/93‑1 (MQ=255)
                                                 gCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGttt                                       >  1:3110444/1‑93 (MQ=255)
                                                  cTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTc                                      <  2:2425760/93‑1 (MQ=255)
                                                  cTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTc                                      >  1:2628597/1‑93 (MQ=255)
                                                  cTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTc                                      <  1:3029149/93‑1 (MQ=255)
                                                    cACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTg                                    >  2:1600358/1‑93 (MQ=255)
                                                    cACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTg                                    <  1:178356/93‑1 (MQ=255)
                                                    cACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTg                                    >  1:4447886/1‑93 (MQ=255)
                                                     aCCTTTTGGCGAACTGGTGTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCAGCAGAAAGCAATGCCAGCCCACGGGCCTTCAGCGTTTCTgg                                   <  1:3282186/93‑1 (MQ=255)
                                                      ccTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGa                                  <  1:4127034/93‑1 (MQ=255)
                                                       cTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGAt                                 <  2:1544020/93‑1 (MQ=255)
                                                           tGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCAcc                             <  2:618313/93‑1 (MQ=255)
                                                            ggCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCAccc                            >  2:4729562/1‑93 (MQ=255)
                                                            ggCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCAccc                            <  1:3314390/93‑1 (MQ=255)
                                                              cGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCgg                          >  2:2930549/1‑93 (MQ=255)
                                                              cGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCGTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGGTTCTGGATCACCCgg                          <  2:2655000/93‑1 (MQ=255)
                                                                aaCTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCa                        >  1:19416/1‑93 (MQ=255)
                                                                  cTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCCGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGa                      >  1:904258/1‑93 (MQ=255)
                                                                   tGGTTTTCACTGAGATCCCCCGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGaa                     >  1:2627893/1‑93 (MQ=255)
                                                                     gTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACg                   <  2:4342029/93‑1 (MQ=255)
                                                                     gTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGAGCATTCAGCGTTTCTGGATCACCCGGCAGAACg                   >  1:4574797/1‑93 (MQ=255)
                                                                      ttttCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGa                  >  1:2580459/1‑93 (MQ=255)
                                                                        ttCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGAtt                <  1:39200/93‑1 (MQ=255)
                                                                        ttCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGAtt                <  2:4644577/93‑1 (MQ=255)
                                                                         tCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGAttt               <  2:348641/93‑1 (MQ=255)
                                                                         tCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGAttt               <  2:1941981/93‑1 (MQ=255)
                                                                            cTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCaa            >  2:2347939/1‑93 (MQ=255)
                                                                             tGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGGCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAat           <  1:4492594/93‑1 (MQ=255)
                                                                              gagaTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAata          <  1:3883157/93‑1 (MQ=255)
                                                                                gaTCCCCGGCAGGCTGGGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATg        <  2:3504463/93‑1 (MQ=255)
                                                                                gaTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATg        <  1:2148302/93‑1 (MQ=255)
                                                                                gaTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATg        <  2:2946104/93‑1 (MQ=255)
                                                                                gaTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATg        >  2:310716/1‑93 (MQ=255)
                                                                                 aTCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGt       >  2:3349517/1‑93 (MQ=255)
                                                                                  tCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGtt      <  1:3359763/93‑1 (MQ=255)
                                                                                  tCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGtt      >  2:4541971/1‑93 (MQ=255)
                                                                                   ccccGGCAGGCTGAGTCCTTCATCAGAAAGCAAACCCCACGCACGGGCCTTCAGCGGTTCTGGATCAGCCCGCAGACAGATTTCAATATGTTc     >  2:1123741/1‑93 (MQ=255)
                                                                                     ccGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTccc   <  2:721553/93‑1 (MQ=255)
                                                                                      cGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCa  <  2:3186250/93‑1 (MQ=255)
                                                                                         |                                                                                         
GTTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGTAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCA  >  NC_000913/1957362‑1957540

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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
TAGTTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGTAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCT  >  NC_000913/1957360‑1957544
                                                                                           |                                                                                             
TAGTTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAG                                                                                             <  SRR4302544.1859657/94‑1 (MQ=60)
 AGTTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGG                                                                                            >  SRR4302544.141438/1‑94 (MQ=60)
   TTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCT                                                                                          <  SRR4302544.2523666/94‑1 (MQ=60)
    TGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTG                                                                                         <  SRR4302544.2124385/94‑1 (MQ=60)
     GTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGA                                                                                        <  SRR4302544.4730429/94‑1 (MQ=60)
        TTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTC                                                                                     >  SRR4302544.1044357/1‑94 (MQ=60)
             ATCGGTGACCGCCAGCGTCGGATTCGTCAATCGCTCATACTCACCTTTTGGCGAACTGGTTTTCACTGAGGTCCCCGGCAGGCTGCGTCCTTCA                                                                                >  SRR4302544.3053191/1‑94 (MQ=60)
             ATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCA                                                                                >  SRR4302544.809454/1‑94 (MQ=60)
               CGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATC                                                                              <  SRR4302544.3853527/94‑1 (MQ=60)
               CGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATC                                                                              >  SRR4302544.4289994/1‑94 (MQ=60)
                 GTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAG                                                                            <  SRR4302544.3124935/94‑1 (MQ=60)
                  TGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGA                                                                           >  SRR4302544.4725059/1‑94 (MQ=60)
                  GGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGA                                                                           <  SRR4302544.1026051/94‑1 (MQ=60)
                    ACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGCGTCCCTCATCAaaac                                                                         >  SRR4302544.4409192/1‑90 (MQ=60)
                      CGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGC                                                                       <  SRR4302544.620853/94‑1 (MQ=60)
                      CGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGC                                                                       >  SRR4302544.1680847/1‑94 (MQ=60)
                       GCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCA                                                                      >  SRR4302544.4337461/1‑94 (MQ=60)
                        CCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAA                                                                     <  SRR4302544.3268520/94‑1 (MQ=60)
                          AGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATG                                                                   >  SRR4302544.494951/1‑94 (MQ=60)
                           GCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGC                                                                  <  SRR4302544.2174263/94‑1 (MQ=60)
                            CGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCC                                                                 >  SRR4302544.4674276/1‑94 (MQ=60)
                            CGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCC                                                                 >  SRR4302544.433798/1‑94 (MQ=60)
                             GTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCA                                                                >  SRR4302544.80107/1‑94 (MQ=60)
                              TCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAG                                                               >  SRR4302544.134379/1‑94 (MQ=60)
                              TCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAG                                                               <  SRR4302544.4248778/94‑1 (MQ=60)
                               CGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGC                                                              <  SRR4302544.3807364/94‑1 (MQ=60)
                                GGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCG                                                             <  SRR4302544.2928095/94‑1 (MQ=60)
                                  ATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGGTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCA                                                           >  SRR4302544.846245/1‑94 (MQ=60)
                                      TGTAATCGCTCATGCTCGCCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGG                                                       <  SRR4302544.4549002/94‑1 (MQ=60)
                                       GTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGC                                                      >  SRR4302544.1723639/1‑94 (MQ=60)
                                        TAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCA                                                     >  SRR4302544.3791464/1‑94 (MQ=60)
                                         AATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCAT                                                    >  SRR4302544.3792303/1‑94 (MQ=60)
                                           TCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTC                                                  <  SRR4302544.4585008/94‑1 (MQ=60)
                                           TCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTC                                                  <  SRR4302544.4501270/94‑1 (MQ=60)
                                            CGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCA                                                 <  SRR4302544.1038302/94‑1 (MQ=60)
                                              CTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGC                                               >  SRR4302544.95944/1‑94 (MQ=60)
                                               TCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCG                                              <  SRR4302544.2326862/94‑1 (MQ=60)
                                               TCATGCTCACCTTATGGCGCACTGATTTTCACTGAGATCCCAGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAaag                                              <  SRR4302544.4246358/94‑4 (MQ=60)
                                                 ATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTT                                            <  SRR4302544.4519687/94‑1 (MQ=60)
                                                  TGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTT                                           >  SRR4302544.3140039/1‑94 (MQ=60)
                                                   GCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTC                                          >  SRR4302544.2653797/1‑94 (MQ=60)
                                                    CTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCT                                         <  SRR4302544.3058169/94‑1 (MQ=60)
                                                     TCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTG                                        >  SRR4302544.4492700/1‑94 (MQ=60)
                                                      CACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGG                                       <  SRR4302544.180120/94‑1 (MQ=60)
                                                        CCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGAT                                     <  SRR4302544.4167964/94‑1 (MQ=60)
                                                              GGCGAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCG                               <  SRR4302544.3345953/94‑1 (MQ=60)
                                                                 GAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCA                            >  SRR4302544.19609/1‑94 (MQ=60)
                                                                 GAACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCA                            >  SRR4302544.2085428/1‑94 (MQ=60)
                                                                   ACTGGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCCGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGA                          >  SRR4302544.913210/1‑94 (MQ=60)
                                                                    CTGGTTTTCACTGAGATCCCCCGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAA                         >  SRR4302544.2653089/1‑94 (MQ=60)
                                                                    CCGGTTTCCACTGAGATCCCCGGCAGGCTGAGGCCTTCACCAGAAAGCAATGCCAGCGCACGGGCAGCCAGCGGGTCAGGATCACCCGGCAGAA                         >  SRR4302544.3647503/1‑94 (MQ=48)
                                                                      GGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACG                       <  SRR4302544.3013888/94‑1 (MQ=60)
                                                                      GGTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGAGCATTCAGCGTTTCTGGATCACCCGGCAGAACG                       >  SRR4302544.4620382/1‑94 (MQ=60)
                                                                       GTTTTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGA                      >  SRR4302544.2605256/1‑94 (MQ=60)
                                                                          TTCACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTT                   <  SRR4302544.39588/94‑1 (MQ=60)
                                                                            CACTGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCA                 <  SRR4302544.3228628/94‑1 (MQ=60)
                                                                               TGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGGCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATA              <  SRR4302544.4537694/94‑1 (MQ=60)
                                                                               TGAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATA              <  SRR4302544.4537695/94‑1 (MQ=60)
                                                                                GAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATAT             <  SRR4302544.4442532/94‑1 (MQ=60)
                                                                                GAGATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATAT             <  SRR4302544.3921728/94‑1 (MQ=60)
                                                                                  GATCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGT           <  SRR4302544.2168760/94‑1 (MQ=60)
                                                                                    TCCCCGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTC         <  SRR4302544.3391825/94‑1 (MQ=60)
                                                                                        CGGCAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAA     >  SRR4302544.1542440/1‑94 (MQ=60)
                                                                                           CAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCT  >  SRR4302544.4736072/1‑94 (MQ=60)
                                                                                           |                                                                                             
TAGTTGTTTTGCCATCGGTGACCGCCAGCGTCGGATTCGGTAATCGCTCATGCTCACCTTTTGGCGAACTGGTTTTCACTGAGATCCCCGGTAGGCTGAGTCCTTCATCAGAAAGCAATGCCAGCGCACGGGCATTCAGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCT  >  NC_000913/1957360‑1957544

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