Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,964,425 G→C 8.9% A209G (GCC→GGC)  flhA ← flagellar biosynthesis protein FlhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,964,4250GC 8.9% 256.3 / 12.3 79A209G (GCC→GGC) flhAflagellar biosynthesis protein FlhA
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base C (3/4);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGG  >  NC_000913/1964276‑1964570
                                                                                                                                                     |                                                                                                                                                 
ccAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATgg                                                                                                                                                   >  1:928621/1‑150 (MQ=255)
 cAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGc                                                                                                                                                  <  2:1028637/150‑1 (MQ=255)
  aGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCa                                                                                                                                                 <  1:1437657/150‑1 (MQ=255)
     gcCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATcg                                                                                                                                              <  2:1070338/150‑1 (MQ=255)
         gAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccgg                                                                                                                                          <  1:908964/150‑2 (MQ=255)
          aaTTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccggg                                                                                                                                         <  1:773573/150‑3 (MQ=255)
           aTTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATCGCCGGGa                                                                                                                                        <  1:461045/150‑1 (MQ=255)
             ttGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTGGGGCGCGTGTCCCTTGCTGTTGCCATGTTGCAGCACGCCGACCAGCAACCCTCCGACAATGTTAATGAGCATGATGAGGATCCCGGCGATGGCATCGCCGCGAAc                                                                                                                                      <  2:269387/150‑1 (MQ=255)
              tgtgCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACa                                                                                                                                     >  1:435310/1‑150 (MQ=255)
              tgtgCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACAAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATCGCCGGGatcc                                                                                                                                     <  2:1215456/150‑4 (MQ=255)
               gtgCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACaa                                                                                                                                    >  1:735720/1‑150 (MQ=255)
               gtgCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACaa                                                                                                                                    >  2:55992/1‑150 (MQ=255)
                 gCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAac                                                                                                                                  <  2:206486/150‑1 (MQ=255)
                 gCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAac                                                                                                                                  <  2:29139/150‑1 (MQ=255)
                  ccaccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAact                                                                                                                                 <  1:317295/150‑1 (MQ=255)
                     ccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAacttac                                                                                                                              >  2:373571/1‑150 (MQ=255)
                         aCCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGc                                                                                                                          <  2:858162/150‑1 (MQ=255)
                          ccGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGcc                                                                                                                         >  1:445068/1‑150 (MQ=255)
                           cGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGccc                                                                                                                        <  1:1224126/150‑1 (MQ=255)
                                ccAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTcc                                                                                                                   >  2:1231921/1‑150 (MQ=255)
                                ccAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTcc                                                                                                                   <  1:838392/150‑1 (MQ=255)
                                       tCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGc                                                                                                            <  1:554161/150‑1 (MQ=255)
                                               gTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTaaaaa                                                                                                    >  2:371355/1‑150 (MQ=255)
                                                    aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAACCTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGc                                                                                               >  2:697976/1‑150 (MQ=255)
                                                    aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGc                                                                                               <  1:55992/150‑1 (MQ=255)
                                                     cTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCt                                                                                              >  1:1288770/1‑150 (MQ=255)
                                                     cTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCt                                                                                              >  1:1283456/1‑150 (MQ=255)
                                                          cgccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCTCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTg                                                                                         <  2:257121/150‑1 (MQ=255)
                                                            ccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAg                                                                                       <  2:235015/150‑1 (MQ=255)
                                                                    cccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCg                                                                               >  2:1056783/1‑150 (MQ=255)
                                                                    cccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCg                                                                               >  1:1228601/1‑150 (MQ=255)
                                                                     ccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCgg                                                                              <  2:127239/150‑1 (MQ=255)
                                                                     ccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCgg                                                                              >  2:752078/1‑150 (MQ=255)
                                                                        tGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAgc                                                                           >  1:716694/1‑150 (MQ=255)
                                                                            cATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAgcggcg                                                                       <  2:527877/150‑1 (MQ=255)
                                                                             aTGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGt                                                                      >  1:1325819/1‑150 (MQ=255)
                                                                                ccATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtttt                                                                   <  2:435310/150‑1 (MQ=255)
                                                                                     ttGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGc                                                              >  1:679925/1‑150 (MQ=255)
                                                                                     ttGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGc                                                              <  1:371355/150‑1 (MQ=255)
                                                                                      tGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGcc                                                             <  1:373571/150‑1 (MQ=255)
                                                                                      tGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGcc                                                             <  2:1288770/150‑1 (MQ=255)
                                                                                      tGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGcc                                                             <  2:1283456/150‑1 (MQ=255)
                                                                                      tGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGcc                                                             <  2:521478/150‑1 (MQ=255)
                                                                                       gcagcaCGCCGACCAGCACCCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGGACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCt                                                            >  2:684403/1‑150 (MQ=255)
                                                                                       gcagcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCt                                                            >  2:1317497/1‑150 (MQ=255)
                                                                                       gcagcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCt                                                            >  2:315784/1‑150 (MQ=255)
                                                                                        cagcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTc                                                           <  2:928621/150‑1 (MQ=255)
                                                                                         agcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCa                                                          >  2:861955/1‑150 (MQ=255)
                                                                                         agcaCCCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCa                                                          >  1:1307291/1‑150 (MQ=255)
                                                                                             cGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCtt                                                      >  1:1442888/1‑150 (MQ=255)
                                                                                              gCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTc                                                     <  1:898628/150‑1 (MQ=255)
                                                                                                 gACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCAcc                                                  <  1:315784/150‑1 (MQ=255)
                                                                                                   ccAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACcaa                                                <  1:1403840/150‑1 (MQ=255)
                                                                                                   ccAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACcaa                                                <  1:931872/150‑1 (MQ=255)
                                                                                                         aCCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCaatcaatc                                          <  1:957657/150‑1 (MQ=255)
                                                                                                          cccGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATcc                                         >  1:220885/1‑150 (MQ=255)
                                                                                                              ccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCGTCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCt                                                            <  2:14848/127‑1 (MQ=255)
                                                                                                              ccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCGTCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCt                                                            >  1:14848/1‑127 (MQ=255)
                                                                                                                  cAATGTTAATGACCATGATGAGGATCGCGGCGATGGCATCGCCGCGAAGAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTaa                                 <  1:462376/150‑1 (MQ=255)
                                                                                                                  cAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTAGCCTCATCTTCACCAATCAATCCGGCGTTaa                                 >  2:1406124/1‑150 (MQ=255)
                                                                                                                     tGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGt                              >  1:1435217/1‑150 (MQ=255)
                                                                                                                       ttAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCg                            >  1:825732/1‑150 (MQ=255)
                                                                                                                        tAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCgg                           >  2:836840/1‑150 (MQ=255)
                                                                                                                        tAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCgg                           >  2:842319/1‑150 (MQ=255)
                                                                                                                          aTGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCg                         <  1:257208/150‑1 (MQ=255)
                                                                                                                           tGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGt                        <  2:1325819/150‑1 (MQ=255)
                                                                                                                           tGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGt                        <  1:1500791/150‑1 (MQ=255)
                                                                                                                               cATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCa                                                    <  1:691171/118‑1 (MQ=255)
                                                                                                                               cATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCa                                                    >  2:691171/1‑118 (MQ=255)
                                                                                                                                       ggATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCt            >  1:881048/1‑150 (MQ=255)
                                                                                                                                       ggATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCt            >  2:541992/1‑150 (MQ=255)
                                                                                                                                        gttcgcGGCGATGCCATCGCCGGGAACAAACTTACTTGCCCCCTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACATCGGAGCGGCGTATTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAACGACGCCGTCAATCGACATCTg           >  2:960870/6‑150 (MQ=255)
                                                                                                                                        gATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGCAGCGGCGTTTTTTCGCCTCATGTTAACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTg           >  1:471686/1‑150 (MQ=255)
                                                                                                                                           cgcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCGTCACCAATCAATCCGGCGGTAAGGTCGGCGTCAATCGCCATCTGCtt        >  2:964339/3‑150 (MQ=255)
                                                                                                                                           cccGGCGATGGCATCGCCGCGAACAAACTTACTTGCCGCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCtt        >  1:781989/1‑150 (MQ=255)
                                                                                                                                           cccGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCtt        >  1:14799/1‑150 (MQ=255)
                                                                                                                                            ccGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCCTTAAGCTCGGCGTCAATCGCCATCTGCTTa       >  2:1085186/1‑150 (MQ=255)
                                                                                                                                               gCGATGCCATCGCCGGGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTAccc    >  2:1444674/1‑150 (MQ=255)
                                                                                                                                                cGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCg   <  1:453551/150‑1 (MQ=255)
                                                                                                                                                 gATGGCATCGCCACGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCgg  <  2:699282/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGG  >  NC_000913/1964276‑1964570

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