Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,242 A→T 6.5% L251F (TTA→TTT 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,2420AT 6.5% 256.3 / 6.7 92L251F (TTA→TTTlsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (52/33);  new base T (3/3);  total (55/37)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

GAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGGGGC  >  NC_000913/1602093‑1602388
                                                                                                                                                     |                                                                                                                                                  
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:669418/1‑150 (MQ=255)
                                                            ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                        >  2:475160/1‑150 (MQ=255)
                                                            ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                        <  2:1636227/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                  
GAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGGGGC  >  NC_000913/1602093‑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: