Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,239 A→T 6.7% K250N (AAA→AAT lsrA → autoinducer 2 import ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,602,2390AT 6.7% 256.9 / 6.6 89K250N (AAA→AATlsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (50/32);  new base T (3/3);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

CGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGGGGC  >  NC_000913/1602091‑1602388
                                                                                                                                                    |                                                                                                                                                     
cGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCACAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCGTCACCCCAGCGGTACGGGAAACATCGCTCTCTGCAAGCCTAAAAt                                                                                                                                                      >  1:95411/1‑150 (MQ=255)
  gAAATTCGCCAGATTGCCTATCGAATTCGCGGGATGCGCGACGGAACCATCGACTTAAGCGGCAAAACCAGCGAAGTTTCTTCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTa                                                                                                                                                    <  1:187339/150‑1 (MQ=255)
     aTTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATgg                                                                                                                                                 <  1:834144/150‑1 (MQ=255)
     aTTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATgg                                                                                                                                                 <  1:1669057/150‑1 (MQ=255)
        cGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAACACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTCCGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTg                                                                                                                                              >  1:1621666/1‑150 (MQ=255)
        cGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAATTATGGCTgg                                                                                                                                             <  1:1374354/150‑1 (MQ=255)
         gCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCCCCCCAGCGGTACGGGAAAAAGCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                                                                             >  1:1361428/1‑150 (MQ=255)
           cAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                                                                           >  2:141149/1‑150 (MQ=255)
           cAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                                                                           >  2:1817742/1‑150 (MQ=255)
           cAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                                                                           <  2:570586/150‑1 (MQ=255)
            tgtttGCAGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                                                                          <  2:344264/147‑4 (MQ=255)
             gATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcaga                                                                                                                                         <  1:933410/150‑5 (MQ=255)
               ttGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcagaga                                                                                                                                       <  1:228194/150‑7 (MQ=255)
                  ccGATCGAATTAGCGTGATGCTCGACGGAAGAATCGCCTTAAGCGGCAAACCCAGCGAACTGGCTACCGAAGACATTATTCAGGCGATCACCCGAGCGGTACGGGAAAAAGCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTg                                                                                                                                    <  1:101020/150‑1 (MQ=255)
                  ccGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTg                                                                                                                                    >  2:507761/1‑150 (MQ=255)
                   cGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTgg                                                                                                                                   >  1:586135/1‑150 (MQ=255)
                   cGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACAAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTgg                                                                                                                                   >  1:594037/1‑150 (MQ=255)
                    gATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCCCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                                                                                                                                  >  1:1333039/1‑150 (MQ=255)
                    gATCGAATTAGCGTGATGCGCGACGGAACAATCGCCTTAAGAGGCAAAACCAGAGAACTGTCTACCGACGTGATTATTCAGGACATCACCCCAGCGGTACGTGATCCATCGCTCTCTGCCCGCCAAAAAGTATGGCGCGAGTTACCTTGt                                                                                                                                  >  1:290736/1‑150 (MQ=255)
                        gAAATACCGTGATGCGCCACCGAACAATCGCCTTAAGCGGCAAAACCAGTCAACAGTCTACAGACGACAATAATCAGGCCATCACACCATCGGTACGGGAAAAATAGCTCTCTGCCAGACAAAAATTAAGCCTGGAGTTACCTGGTAAcc                                                                                                                              >  2:497207/1‑150 (MQ=255)
                              gCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCTGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcac                                                                                                                        >  2:386833/1‑150 (MQ=255)
                                gTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCCGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                                                                      >  1:1627354/1‑150 (MQ=255)
                                    tGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCCCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg                                                                                                                  >  2:707992/1‑150 (MQ=255)
                                       gcgACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccg                                                                                                               >  1:1454633/1‑150 (MQ=255)
                                        cgACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGAGGTACGGGAAAAATCGCTCTCTGCCAGGCAAAAATTATGGATGGAGTTACCTGGTAACCGACACAACATgccgcc                                                                                                             <  2:474210/150‑1 (MQ=255)
                                           cGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACATCATGCCGCCgg                                                                                                           >  2:188194/1‑150 (MQ=255)
                                            ggAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCAACCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGa                                                                                                          >  1:1630919/1‑150 (MQ=255)
                                             gAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGCCATTATTCAGGCCATCCCCCCAGCGGAACGGGAAAAATCGCTCTCTGCCAGCCAAAAAATATGGCTGGAGTTACCTGGTAACCGGCCACAACATGCCGCGGGaa                                                                                                         >  2:1109795/1‑150 (MQ=255)
                                               aCCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACg                                                                                                       >  2:1398629/1‑150 (MQ=255)
                                                   tCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGCCGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAACATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTCCCTGGTAACCGCCCACAACATGCCGCCGGAACGCCgg                                                                                                   >  1:1452414/1‑150 (MQ=255)
                                                   tCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCTGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCgg                                                                                                   >  1:849862/1‑150 (MQ=255)
                                                    cGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGACGCCGGAACGCCGGt                                                                                                  >  1:1363450/1‑150 (MQ=255)
                                                     gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTg                                                                                                 >  1:1426689/1‑150 (MQ=255)
                                                      ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGc                                                                                                <  1:1605618/150‑1 (MQ=255)
                                                          aatcGGCAAAACCATCGAACTGTCAACCGAAGACATTATTCAGCCCATCACCCCACCGTTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTTCCTGTTAACCGCCCACAACATGCCCCCGGAACGACGGTGCTGac                                                                                            <  2:1643121/147‑1 (MQ=255)
                                                              ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGACCGCCGGTGCTGACACTg                                                                                        >  1:869238/1‑150 (MQ=255)
                                                              ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                        >  2:475160/1‑150 (MQ=255)
                                                              ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                        <  2:1636227/150‑1 (MQ=255)
                                                              ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                        >  2:669418/1‑150 (MQ=255)
                                                                      cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                                                                                >  2:1388968/1‑150 (MQ=255)
                                                                      cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGAAACGCCGGTGCTGACACTGGAAAATCt                                                                                <  1:176792/150‑1 (MQ=255)
                                                                        gCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACTCTGGAAAATCTgc                                                                              >  1:505297/1‑149 (MQ=255)
                                                                           aaCTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCg                                                                           <  2:604882/150‑1 (MQ=255)
                                                                                tCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGaa                                                                      >  1:574187/1‑150 (MQ=255)
                                                                                tCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAg                                                                     <  2:460636/150‑1 (MQ=255)
                                                                                 cTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAg                                                                     >  1:983721/1‑150 (MQ=255)
                                                                                  tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAgg                                                                    >  1:40639/1‑150 (MQ=255)
                                                                                     cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACATGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt                                                                 >  2:833134/1‑150 (MQ=255)
                                                                                       acgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa                                                               <  1:141149/150‑1 (MQ=255)
                                                                                       acgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa                                                               <  1:1817742/150‑1 (MQ=255)
                                                                                        cgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg                                                              >  1:546219/1‑150 (MQ=255)
                                                                                        cgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg                                                              <  2:1333039/150‑1 (MQ=255)
                                                                                        atccATTCTTCAGTCAATAACACCATATGTACCGGAAAAATATCTCTCTTCCATCCAAAAATTCTGGCTCGAGTTACCTGGTAAACGCAGAGAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg                                                              <  2:290736/147‑1 (MQ=255)
                                                                                         gacATTATTCAGGCCATCACCCCAGCTGTACGGGAAAAATCGCTCTCTGCCAGACAAACATTATGTCTGGAGTAACCTGGTAACCGCCCACACCATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAgg                                                             <  1:469412/150‑1 (MQ=255)
                                                                                         gacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAgg                                                             >  1:1341934/1‑150 (MQ=255)
                                                                                           cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa                                                           <  2:942374/150‑1 (MQ=255)
                                                                                               attCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTc                                                       >  1:183952/1‑150 (MQ=255)
                                                                                                ttCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCa                                                      >  2:1506504/1‑150 (MQ=255)
                                                                                                     gCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACAGGCGAAGGTTTCAGGAATGTCAGCCTg                                                 >  1:1778487/1‑150 (MQ=255)
                                                                                                      ccATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAATTCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGa                                                <  2:1783039/150‑1 (MQ=255)
                                                                                                       cATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGAc                                               >  1:1342543/1‑150 (MQ=255)
                                                                                                            ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCa                                          <  1:1388968/150‑1 (MQ=255)
                                                                                                            ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCa                                          <  2:1454633/150‑1 (MQ=255)
                                                                                                              ccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAAt                                        <  1:56007/150‑1 (MQ=255)
                                                                                                               cAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAAAGACAGCCTGACGCTCAATg                                       >  1:435453/1‑150 (MQ=255)
                                                                                                                  cGGTACGGGAAACATCGCGCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCg                                    <  1:1720027/150‑1 (MQ=255)
                                                                                                                        gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGCCGCTCAATGCCGGAGaaa                              >  1:1593664/1‑150 (MQ=255)
                                                                                                                        gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGaaa                              >  2:1332281/1‑150 (MQ=255)
                                                                                                                         ggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGTAAATCTGACCGGCGCAGGTGTCAGGAATGTCAGCCTGCCGCTCAATGCCGGAGAAAt                             >  1:296897/1‑150 (MQ=255)
                                                                                                                          gAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAAtt                            >  2:4528/1‑150 (MQ=255)
                                                                                                                             aaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTg                         <  2:986638/150‑1 (MQ=255)
                                                                                                                              aaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGCTATTCTgg                        >  2:95890/1‑150 (MQ=255)
                                                                                                                                tCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGc                      <  2:842388/150‑1 (MQ=255)
                                                                                                                                 cGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACATCATGCCGCCGGAACGCCGGTGCTGACACTGGATAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCCATGCCGGAGAAATTCTGGGcc                     >  1:1617777/1‑150 (MQ=255)
                                                                                                                                 cGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACACGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGcc                     <  1:67648/150‑1 (MQ=255)
                                                                                                                                   ggtaactccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGCCCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctg                   >  2:1280331/8‑150 (MQ=255)
                                                                                                                                   ctctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctg                   >  2:1635831/1‑150 (MQ=255)
                                                                                                                                   ctctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCAGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctg                   <  2:1350123/150‑1 (MQ=255)
                                                                                                                                     ctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGCCACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGCCGCTCAATGCCGGAGAAATTCTGGGCctggc                 >  2:256945/1‑150 (MQ=255)
                                                                                                                                     ctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggc                 <  2:594037/150‑1 (MQ=255)
                                                                                                                                     ctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggc                 <  2:586135/150‑1 (MQ=255)
                                                                                                                                      aactccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGACGCCGGAACGACGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggct                >  2:1076927/5‑150 (MQ=255)
                                                                                                                                         gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctggg             >  2:1445450/1‑150 (MQ=255)
                                                                                                                                          ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGCCGCTCAATGCCGGAGAAATTCTGGGCCTggctgggc            >  1:1026156/1‑150 (MQ=255)
                                                                                                                                          ccagccaAAACTTATGGCGGGAGTTACCTGGTAACCGCCCACAACCTGCCGCCGGAACGCCGGGGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAACTTCTGGGCCTggctgggc            <  2:1268626/150‑1 (MQ=255)
                                                                                                                                            agccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctgggctg          <  1:514398/150‑1 (MQ=255)
                                                                                                                                             gccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGCCGCTCAATGCCGGAGAAATTCTGGGCCTGgctgggatgg         >  1:1702121/1‑150 (MQ=255)
                                                                                                                                             gccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGgctgggctgg         >  1:1415399/1‑150 (MQ=255)
                                                                                                                                               caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGACTGGGCtggtg       <  1:1196890/150‑1 (MQ=255)
                                                                                                                                                 aaaaTTATGGCTGGAGTTACCTGGTAACCGCACACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTggg     >  2:1232447/1‑150 (MQ=255)
                                                                                                                                                 aaaaTTATGGCTGGAGTTACATGGTAACCGCCCACAACATGCCGCCGGAACGCCGGAGCTGACACTGGAATATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTggg     >  2:110246/1‑150 (MQ=255)
                                                                                                                                                  aaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTgggg    <  2:1426689/150‑1 (MQ=255)
                                                                                                                                                   aaTTATGGCTGGAGTTACCTGGTAACCGCCCTCAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTggggg   <  1:582598/150‑1 (MQ=255)
                                                                                                                                                    aTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGGGGc  >  1:1563797/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGGGGC  >  NC_000913/1602091‑1602388

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: