Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,017,926 Δ1 bp 5.1% coding (1373/1374 nt) fliI → flagellum‑specific ATP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,017,9260A. 5.1% 236.9 / 5.7 78coding (1373/1374 nt)fliIflagellum‑specific ATP synthase
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base . (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CGCGTATGCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTG  >  NC_000913/2017777‑2018075
                                                                                                                                                     |                                                                                                                                                     
cgcgTATGCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATa                                                                                                                                                       <  2:353499/150‑1 (MQ=255)
   gTATGCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAAcc                                                                                                                                                    <  1:788202/150‑1 (MQ=255)
   gTATGCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAAcc                                                                                                                                                    <  1:1206132/150‑1 (MQ=255)
    tATGCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAAccc                                                                                                                                                   <  2:251670/150‑1 (MQ=255)
    tATGCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAAccc                                                                                                                                                   <  2:1335572/150‑1 (MQ=255)
     aTGCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCa                                                                                                                                                  >  1:1276827/1‑150 (MQ=255)
       gCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAgg                                                                                                                                                >  1:320035/1‑150 (MQ=255)
        ccAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGga                                                                                                                                               >  2:743765/1‑150 (MQ=255)
        ccAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGga                                                                                                                                               >  2:439765/1‑150 (MQ=255)
              gCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAAc                                                                                                                                         <  2:217509/150‑1 (MQ=255)
              gCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGCAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTACCCGACAGTGTCATAACCCAGGAGATAAc                                                                                                                                         <  2:879480/150‑1 (MQ=255)
                 gCGAGCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGGCTCTCCCGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGc                                                                                                                                      <  2:252789/150‑1 (MQ=255)
                   gATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCAAATTGCAACAAGCCATTTTAGAACGCGCGGACTCGGAAGCCGCTCTCCACGGCCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGAGAACGGCAg                                                                                                                                    >  1:1261475/1‑150 (MQ=255)
                   gATCCGATGCTCGATAAAGCAATCGCCCTGTGGCGGCAGCTGGAGGGCTAATAGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCCTCTCTCCAGGGGCTGGAGCGTATTTTCCCGCCAGTGTCATAACCCAGGAGAAAACGGCac                                                                                                                                    >  1:1266691/1‑149 (MQ=255)
                    aTCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATACCCAGGAGATAACGGCAGAt                                                                                                                                  >  2:793349/1‑150 (MQ=255)
                    aTCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATACCCAGGAGATAACGGCAGAt                                                                                                                                  >  2:703562/1‑150 (MQ=255)
                      ccGATGCTCGATAAAGACATCGCCCTGTGGCCGGAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGTAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATg                                                                                                                                 <  1:249549/150‑1 (MQ=255)
                        gATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGc                                                                                                                               <  2:1150717/150‑1 (MQ=255)
                        gATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGc                                                                                                                               <  2:373024/150‑1 (MQ=255)
                        gATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCAAAACACAGGAGATAACGGCAGATGGc                                                                                                                               >  2:845324/1‑150 (MQ=255)
                           gCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGa                                                                                                                            <  2:529000/150‑1 (MQ=255)
                              cGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACa                                                                                                                         <  2:22278/150‑1 (MQ=255)
                                  aaaGCCATGGCCCTGGGGCCGCAGCGGGAGGGCTTTTTGCAACAAGCCATTTGTCAACGCGCGGACTGGGAAGCGGCTCTCCAGGGGCTGGAGCGTATTTACCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGt                                                                                                                     <  2:230960/150‑1 (MQ=255)
                                   aaGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACACGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTg                                                                                                                    <  2:260792/150‑1 (MQ=255)
                                   aaGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTg                                                                                                                    <  1:500642/150‑1 (MQ=255)
                                    aGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGTAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTgc                                                                                                                   >  1:275418/1‑150 (MQ=255)
                                      ccATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGTGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACTGCAGATGGCAGAACATGGTgcgc                                                                                                                 >  1:648881/1‑150 (MQ=255)
                                      ccATCGCCCTGTGGCCGAAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTgcgc                                                                                                                 >  1:747639/1‑150 (MQ=255)
                                        aTCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGGCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTg                                                                                                               <  1:599329/150‑1 (MQ=255)
                                                  ggCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGAGCGCTGGCGACCCTGa                                                                                                     >  1:1085158/1‑150 (MQ=255)
                                                    ccGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGaaa                                                                                                   <  1:745762/150‑1 (MQ=255)
                                                    ccGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGaaa                                                                                                   >  1:315154/1‑150 (MQ=255)
                                                         gctgtcgggCTATTTGCAACAAGGCATTTTTCAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCt                                                                                              <  2:255552/144‑1 (MQ=255)
                                                          cTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTg                                                                                             >  1:924108/1‑150 (MQ=255)
                                                            ggAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGc                                                                                           <  1:1071463/150‑1 (MQ=255)
                                                             gAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCa                                                                                          <  2:1204135/150‑1 (MQ=255)
                                                                 gCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCATGAAAGATCTGGCAGaaa                                                                                      >  2:1009205/1‑150 (MQ=255)
                                                                      ttGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAagag                                                                                 <  1:585696/150‑1 (MQ=255)
                                                                       tGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAgg                                                                                >  2:374064/1‑150 (MQ=255)
                                                                          aacaaGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTag                                                                             <  2:1276827/150‑1 (MQ=255)
                                                                           acaaGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGGGACCCTGAAAGATCTGGCAGAAAAAGAGGTtga                                                                            >  2:443168/1‑150 (MQ=255)
                                                                           acaaGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTaga                                                                            >  1:354650/1‑150 (MQ=255)
                                                                             aaGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGATAAAGAGGTAGAgg                                                                          >  1:1305113/1‑150 (MQ=255)
                                                                             aaGGAATTTTTCAACGCGCGCCATGGGACGCGGCTCTACAGTGTCTGGAGCGTATTTTCCCGACCGTGTCATCACCCCGGAGATAACGGCAGATGCCACAACATGGCGCGCTGGCGACCCTGAAAGTTCTGGCAGAAAAAGAGGTAGAgg                                                                          <  2:197924/150‑1 (MQ=255)
                                                                              aGGCATTTTTGAATGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGa                                                                         >  1:749609/1‑150 (MQ=255)
                                                                              aGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGa                                                                         <  1:439765/150‑1 (MQ=255)
                                                                               ggCTTTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGAGAGTGTCATAACCCATGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAACGAGGTAGAGGAt                                                                        <  2:1279614/150‑1 (MQ=255)
                                                                                       tGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCcgcgcg                                                                >  1:733035/1‑150 (MQ=255)
                                                                                         aaCGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGcc                                                              >  1:585141/1‑150 (MQ=255)
                                                                                             cgcgGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCctgct                                                          <  1:1189754/150‑1 (MQ=255)
                                                                                              gcgGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGAATTTTCCCGACAGTGTCATAACCCAGTAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCctgctg                                                         <  1:639366/150‑1 (MQ=255)
                                                                                               cgGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTgg                                                        >  2:628256/1‑150 (MQ=255)
                                                                                                    tGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGaa                                                   >  2:1088408/1‑150 (MQ=255)
                                                                                                      ggAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAAt                                                 <  2:1179922/150‑1 (MQ=255)
                                                                                                       gAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATg                                                <  1:1068472/150‑1 (MQ=255)
                                                                                                             gctctcCAGGGGCTGTAGCGAATTATCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGTCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTcgc                                          <  2:914927/149‑1 (MQ=255)
                                                                                                              ctctcCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTcgcg                                         >  2:1070271/1‑150 (MQ=255)
                                                                                                               tctcCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCgg                                        >  1:288742/1‑150 (MQ=255)
                                                                                                                    tgtgtcTGTAGCGTATTTTCCCTACAGTGTCATAACCTTGGAGTTAACGGGATATGTCATAACATTGTGCGCTGTCGACCCTGAAAGATCTCGAAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCTGATGTc                                   <  2:665391/145‑1 (MQ=255)
                                                                                                                    agtggCTGGAGCGTATTTTCCCGACTGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGTTCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTc                                   <  2:630780/147‑1 (MQ=255)
                                                                                                                    aGGGGCGGGAGCGTATTTTCCCGACAGTGTCATACCCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTc                                   <  1:1145222/150‑1 (MQ=255)
                                                                                                                      gggCTGGAGCGTATTTTCCCGACAGTGTCATACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTcagc                                <  1:703562/150‑1 (MQ=255)
                                                                                                                      gggCTGGAGCGTATTTTCCCGACAGTGTCATACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTcagc                                <  1:793349/150‑1 (MQ=255)
                                                                                                                             agagTATTTTCCCGTCAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCgg                          >  2:775090/4‑150 (MQ=255)
                                                                                                                             aGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGCGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGAGTCGCGGATGTCAGCAGGCgg                          >  2:664590/1‑150 (MQ=255)
                                                                                                                             aGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCgg                          >  2:564093/1‑150 (MQ=255)
                                                                                                                                     ttCCCGACAGTGTCTTAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTTGCAGAAAACGTGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAg                  <  1:845324/150‑1 (MQ=255)
                                                                                                                                     ttCCCGACAGTGTCATAACCCATGAGATAAAGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCTCGCCTGCTGGGTGAAATGCGTCGCGGATCTCAGCAGGCGGAAGAACAg                  >  1:71455/1‑150 (MQ=255)
                                                                                                                                      tCCCGACAGTGTCATAACCCAGGAGATACCGTCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGc                 <  2:370472/150‑1 (MQ=255)
                                                                                                                                           aCAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCaaa            <  2:320035/150‑1 (MQ=255)
                                                                                                                                              gtgtCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATg         >  2:392665/1‑150 (MQ=255)
                                                                                                                                                gtCATAACCCAGGAGCTAACGGCCGATGGCAGAACATGGTGCGCGGGCGACCCTGCAAGATCGGGCAGAAAACGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCt       <  2:466865/150‑1 (MQ=255)
                                                                                                                                                gtCATAACCCAGGAGATAACGGCAGATGGCAGAAGATGGTGCGCTGGAGACCGTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCt       <  2:1093470/150‑1 (MQ=255)
                                                                                                                                                gtCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCAGAAGAACAGCTCAAAATGCt       <  2:1047328/150‑1 (MQ=255)
                                                                                                                                                 tCATAACCCGGGAGAAAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCtg      <  2:405693/150‑1 (MQ=255)
                                                                                                                                                  cATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCtga     >  2:917657/1‑150 (MQ=255)
                                                                                                                                                  cATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCtga     >  1:608122/1‑150 (MQ=255)
                                                                                                                                                  cATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCtga     >  2:201288/1‑150 (MQ=255)
                                                                                                                                                   aTAACCCAGGAGATAACGTCAGATGGCAGATCATGCTGCGCTTGCGACACTGAAAGATCTGGCATAAAAATATGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCtgat    <  1:1004766/150‑1 (MQ=255)
                                                                                                                                                     aaCCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCtgattg  <  2:430190/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGTATGCCAAAGGCAGCGATCCGATGCTCGATAAAGCCATCGCCCTGTGGCCGCAGCTGGAGGGCTATTTGCAACAAGGCATTTTTGAACGCGCGGACTGGGAAGCGTCTCTCCAGGGGCTGGAGCGTATTTTCCCGACAGTGTCATAACCCAGGAGATAACGGCAGATGGCAGAACATGGTGCGCTGGCGACCCTGAAAGATCTGGCAGAAAAAGAGGTAGAGGATGCCGCGCGCCTGCTGGGTGAAATGCGTCGCGGATGTCAGCAGGCGGAAGAACAGCTCAAAATGCTGATTG  >  NC_000913/2017777‑2018075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: