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 F3 I160 R1 220 104.1 6051714 95.9% 5803593 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,215,724 T→G W657G (TGG→GGG)  metH → methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,215,7240TG99.1% 423.2 / NA 113W657G (TGG→GGG) metHmethionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (41/71);  minor base A (1/0);  total (42/71)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAG  >  NZ_CP009273/4215639‑4215810
                                                                                     |                                                                                      
cGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGa                                                                                    >  1:1024750/1‑90 (MQ=255)
cGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGa                                                                                    <  1:442724/90‑1 (MQ=255)
cGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGa                                                                                    >  2:134432/1‑90 (MQ=255)
cGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGa                                                                                    <  2:487005/90‑1 (MQ=255)
 gAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGaa                                                                                   >  2:1493015/1‑90 (MQ=255)
 gAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGaa                                                                                   <  2:479768/90‑1 (MQ=255)
 gAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGaa                                                                                   >  2:2988450/1‑90 (MQ=255)
  aGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAg                                                                                  <  1:1409275/90‑1 (MQ=255)
  aGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAg                                                                                  <  1:2739095/90‑1 (MQ=255)
  aGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAg                                                                                  <  1:649345/90‑1 (MQ=255)
     gTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGa                                                                               >  1:2036416/1‑90 (MQ=255)
       ttACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAAt                                                                             >  1:1575689/1‑90 (MQ=255)
         aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           <  2:779571/90‑1 (MQ=255)
         aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           <  1:1493015/90‑1 (MQ=255)
         aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           <  1:677778/90‑1 (MQ=255)
         aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           <  2:2111162/90‑1 (MQ=255)
         aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           <  2:1027143/90‑1 (MQ=255)
               gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           <  2:1640332/84‑1 (MQ=255)
               gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           >  1:1640332/1‑84 (MQ=255)
               gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCt                                                                     >  2:2486925/1‑90 (MQ=255)
               gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCt                                                                     >  1:1555995/1‑90 (MQ=255)
                 ttGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCtag                                                                   >  2:70779/1‑88 (MQ=255)
                 ttGCCGAGAAATATCGAGGCAGCAAAACCGACGAGACGGCCAACGCGGAGGAGGCGGAGAGGCGGTCGAGAGAAGAGAAAAAACCTCggg                                                                   >  1:2723567/1‑90 (MQ=255)
                  tGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGa                                                                  >  1:2111190/1‑90 (MQ=255)
                  tGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGa                                                                  >  1:2487383/1‑90 (MQ=255)
                   gCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGaa                                                                 <  2:1776241/90‑1 (MQ=255)
                   gCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGaa                                                                 <  2:2309522/90‑1 (MQ=255)
                    ccGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAAt                                                                <  2:1957437/90‑1 (MQ=255)
                      gagaAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATAc                                                              <  1:134432/90‑1 (MQ=255)
                          aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          <  1:2460924/90‑1 (MQ=255)
                          aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          <  2:1472742/90‑1 (MQ=255)
                          aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          <  1:2261788/90‑1 (MQ=255)
                          aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          <  2:80623/90‑1 (MQ=255)
                          aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          <  1:1968295/90‑1 (MQ=255)
                           atatCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCt                                                         <  2:1547341/90‑1 (MQ=255)
                           atatCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCt                                                         <  2:1257999/90‑1 (MQ=255)
                           atatCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCt                                                         <  2:1167734/90‑1 (MQ=255)
                           atatCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCt                                                         <  1:1277297/90‑1 (MQ=255)
                           atatCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCt                                                         <  1:1390853/90‑1 (MQ=255)
                            tatCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTg                                                        >  1:2008406/1‑90 (MQ=255)
                                gcgGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCa                                                    >  1:1292638/1‑90 (MQ=255)
                                gcgGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCa                                                    >  2:2727800/1‑90 (MQ=255)
                                 cgGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaa                                                   <  1:2937583/90‑1 (MQ=255)
                                  ggCAGCAAAACCGACGACCCCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  <  2:2020782/90‑1 (MQ=255)
                                  ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  <  2:1549556/90‑1 (MQ=255)
                                  ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  <  1:2479597/90‑1 (MQ=255)
                                  ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  <  2:196726/90‑1 (MQ=255)
                                  ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  >  2:2416570/1‑90 (MQ=255)
                                  ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  <  1:2914261/90‑1 (MQ=255)
                                     agcaAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGc                                               <  2:1575689/90‑1 (MQ=255)
                                           aCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTAcc                                         <  2:2413363/90‑1 (MQ=255)
                                           aCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTAcc                                         <  1:281399/90‑1 (MQ=255)
                                           aCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTAcc                                         <  2:1581921/90‑1 (MQ=255)
                                           aCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTAcc                                         <  1:1565327/90‑1 (MQ=255)
                                              gacgacACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAg                                      >  2:2106661/1‑90 (MQ=255)
                                              gacgacACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAg                                      >  2:2011570/1‑90 (MQ=255)
                                                cgacACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGtt                                    >  2:2948938/1‑90 (MQ=255)
                                                   cacCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTAt                                 <  1:403655/90‑1 (MQ=255)
                                                     ccGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCg                               <  1:2948938/90‑1 (MQ=255)
                                                     ccGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCg                               >  2:699281/1‑90 (MQ=255)
                                                     ccGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCg                               >  2:369083/1‑90 (MQ=255)
                                                     ccGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCg                               >  2:2939957/1‑90 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGATGTGAATaaa                                                                          <  2:1528340/44‑1 (MQ=38)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGATGTGAATaaa                                                                          >  1:1528340/1‑44 (MQ=38)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaaa                                                                          >  1:204669/1‑44 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaaa                                                                          <  2:204669/44‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  1:2024252/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  2:1651901/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  1:2178731/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  1:2661432/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  2:2026752/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  2:715123/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  2:1423484/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  2:2180996/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTAGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  1:70779/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGAAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  1:1929332/90‑1 (MQ=255)
                                                         cAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCa                           <  1:2930280/90‑1 (MQ=255)
                                                             gCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGAt                       <  1:2727800/90‑1 (MQ=255)
                                                               ccAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATAc                     <  2:2533886/90‑1 (MQ=255)
                                                                 agcagGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          >  1:2010850/1‑51 (MQ=255)
                                                                 agcagGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCg                   >  2:1927251/1‑90 (MQ=255)
                                                                     ggCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          >  2:745639/1‑47 (MQ=255)
                                                                      gCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAg                          >  2:1971340/1‑78 (MQ=255)
                                                                      gCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCgaagaa              >  1:2406326/1‑90 (MQ=255)
                                                                          ggTGGCGCTCGGGGGAAGTAAATAAACGTCAGGAATCCTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGccc          <  1:634064/89‑1 (MQ=255)
                                                                          aGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATCCTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGccc          <  1:2326859/90‑1 (MQ=255)
                                                                          aGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGGGTTTATCGAGCAGGATACCGAAGAAGccc          <  1:908756/90‑1 (MQ=255)
                                                                          aGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            >  2:2322820/1‑72 (MQ=255)
                                                                          aGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGccc          <  2:1849949/90‑1 (MQ=255)
                                                                          aGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGccc          <  1:667964/90‑1 (MQ=255)
                                                                           gTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCg         <  2:43018/90‑1 (MQ=255)
                                                                           gTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCg         >  2:2430332/1‑90 (MQ=255)
                                                                             ggCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAg                             <  2:1334312/68‑1 (MQ=255)
                                                                             ggCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAg                             >  1:1334312/1‑68 (MQ=255)
                                                                               cgcTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  >  1:541381/1‑45 (MQ=255)
                                                                               cgcTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  <  2:541381/45‑1 (MQ=255)
                                                                               cgcTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcag     >  2:1204207/1‑90 (MQ=255)
                                                                               cgcTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcag     >  1:2509137/1‑90 (MQ=255)
                                                                                  tCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcagcag  >  1:1349770/1‑90 (MQ=255)
                                                                                  tCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcagcag  <  2:2104873/90‑1 (MQ=255)
                                                                                     |                                                                                      
CGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAG  >  NZ_CP009273/4215639‑4215810

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GACGATGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCC  >  NZ_CP009273/4215628‑4215821
                                                                                                |                                                                                                 
GACGATGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGG                                                                                                <  SRR3722030.2422845/100‑1 (MQ=60)
GACGATGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGG                                                                                                <  SRR3722030.155404/100‑1 (MQ=60)
 ACGATGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGA                                                                                               >  SRR3722030.1040213/1‑100 (MQ=60)
   GATGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAG                                                                                             <  SRR3722030.924752/100‑1 (MQ=60)
   GATGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAG                                                                                             <  SRR3722030.575948/100‑1 (MQ=60)
     TGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTG                                                                                           <  SRR3722030.2050342/100‑1 (MQ=60)
     TGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTG                                                                                           <  SRR3722030.1291695/100‑1 (MQ=60)
      GGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGA                                                                                          >  SRR3722030.2072445/1‑100 (MQ=60)
        CACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAAT                                                                                        >  SRR3722030.1601196/1‑100 (MQ=60)
           CGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAA                                                                                     <  SRR3722030.367074/100‑1 (MQ=60)
           CGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAA                                                                                     <  SRR3722030.449450/100‑1 (MQ=60)
             AGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACG                                                                                   <  SRR3722030.2789685/100‑1 (MQ=60)
             AGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACG                                                                                   <  SRR3722030.1431509/100‑1 (MQ=60)
             AGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACG                                                                                   <  SRR3722030.659404/100‑1 (MQ=60)
                GTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCT                                                                                >  SRR3722030.1581166/1‑100 (MQ=60)
                GTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCT                                                                                >  SRR3722030.1666955/1‑100 (MQ=60)
                   TACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGA                                                                             >  SRR3722030.2532879/1‑100 (MQ=60)
                   TACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGA                                                                             >  SRR3722030.2148847/1‑100 (MQ=60)
                    ccaccgagatctacactggccatgtcgtcggcagcgtcagatgtgtataagagacagGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAA                                                                            <  SRR3722030.1107514/43‑1 (MQ=60)
                    ACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAA                                                                            <  SRR3722030.1516939/100‑1 (MQ=60)
                    ACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAA                                                                            <  SRR3722030.688270/100‑1 (MQ=60)
                         caccgagatctacactggccatgtcgtcggcagcgtcagatgtgtataagagaCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTC                                                                       <  SRR3722030.1500252/47‑1 (MQ=60)
                             TGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTG                                                                   >  SRR3722030.2043766/1‑100 (MQ=60)
                                 GAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCA                                                               >  SRR3722030.1312218/1‑100 (MQ=60)
                                 GAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCA                                                               <  SRR3722030.136383/100‑1 (MQ=60)
                                     tcgtcggcagcgtcagatgtgtataagagacagACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGG                                                           <  SRR3722030.757125/67‑1 (MQ=60)
                                     AATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGG                                                           <  SRR3722030.2505891/100‑1 (MQ=60)
                                     AATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGG                                                           <  SRR3722030.2302681/100‑1 (MQ=60)
                                     AATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGG                                                           <  SRR3722030.2002746/100‑1 (MQ=60)
                                      ATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGC                                                          <  SRR3722030.1296497/100‑1 (MQ=60)
                                      ATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGC                                                          <  SRR3722030.1412701/100‑1 (MQ=60)
                                            CGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACC                                                    <  SRR3722030.2991727/100‑1 (MQ=60)
                                             GGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCG                                                   <  SRR3722030.2968010/100‑1 (MQ=60)
                                             GGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCG                                                   <  SRR3722030.2524905/100‑1 (MQ=60)
                                                      ACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCG                                          <  SRR3722030.285487/100‑1 (MQ=60)
                                                      ACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCG                                          <  SRR3722030.1590665/100‑1 (MQ=60)
                                                         GACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGATGTGAATAAACGTCTGGAATctgtctcttatacacatctgacgctgccgacgacat                                       >  SRR3722030.1552997/1‑64 (MQ=60)
                                                         GACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATctgtctcttatacacatctgacgctgccgacgacat                                       >  SRR3722030.207582/1‑64 (MQ=60)
                                                              CACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGAT                                  <  SRR3722030.409730/100‑1 (MQ=60)
                                                                CCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATAC                                <  SRR3722030.3003316/100‑1 (MQ=60)
                                                                  GCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCtgtctcttatacacatctgacgctgccg                              >  SRR3722030.2046256/1‑72 (MQ=60)
                                                                   agagacagCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                             <  SRR3722030.2364975/92‑1 (MQ=60)
                                                                   CCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                             <  SRR3722030.2059984/100‑1 (MQ=60)
                                                                   CCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                             <  SRR3722030.1440858/100‑1 (MQ=60)
                                                                   CCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                             <  SRR3722030.2710464/100‑1 (MQ=60)
                                                                   CCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                             <  SRR3722030.2217837/100‑1 (MQ=60)
                                                                   CCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTAGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                             <  SRR3722030.71857/100‑1 (MQ=60)
                                                                   CCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGAAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                             <  SRR3722030.1962883/100‑1 (MQ=60)
                                                                    CAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAA                            <  SRR3722030.2984313/100‑1 (MQ=60)
                                                                     AGCGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAG                           <  SRR3722030.2005862/100‑1 (MQ=60)
                                                                       CGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAA                         >  SRR3722030.2450169/1‑100 (MQ=60)
                                                                        GCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAG                        <  SRR3722030.2778178/100‑1 (MQ=60)
                                                                              CAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGctgtctcttata                  >  SRR3722030.1354861/1‑88 (MQ=60)
                                                                                GGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGctgtctcttatacacatctgacgctgccgacgaca                >  SRR3722030.549734/1‑65 (MQ=60)
                                                                                GGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAG                >  SRR3722030.2555050/1‑100 (MQ=60)
                                                                                   GGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAG             >  SRR3722030.1370695/1‑100 (MQ=60)
                                                                                    GAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGG            >  SRR3722030.1014569/1‑100 (MQ=60)
                                                                                    GAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGG            >  SRR3722030.463342/1‑100 (MQ=60)
                                                                                    GAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGG            >  SRR3722030.1235573/1‑100 (MQ=60)
                                                                                     GGTGGCGCTCGGGGGAAGTAAATAAACGTCAGGAATCCTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGC           <  SRR3722030.643909/100‑1 (MQ=60)
                                                                                     AGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATCCTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGC           <  SRR3722030.2369087/100‑1 (MQ=60)
                                                                                     AGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGGGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGC           <  SRR3722030.922696/100‑1 (MQ=60)
                                                                                     AGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGC           <  SRR3722030.678315/100‑1 (MQ=60)
                                                                                         GCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACG       >  SRR3722030.1184613/1‑100 (MQ=60)
                                                                                             TCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCC   <  SRR3722030.2740520/100‑1 (MQ=60)
                                                                                              CGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCC  >  SRR3722030.248147/1‑100 (MQ=60)
                                                                                              CGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCC  >  SRR3722030.773047/1‑100 (MQ=60)
                                                                                                |                                                                                                 
GACGATGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCC  >  NZ_CP009273/4215628‑4215821

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