Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,236 A→T 6.7% Q249H (CAA→CAT 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,2360AT 6.7% 258.9 / 5.7 89Q249H (CAA→CATlsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (50/33);  new base T (3/3);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGG  >  NC_000913/1602089‑1602385
                                                                                                                                                   |                                                                                                                                                     
gCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaa                                                                                                                                                     >  2:582598/1‑150 (MQ=255)
  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:1817742/1‑150 (MQ=255)
             cAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                                                                        <  2:570586/150‑1 (MQ=255)
             cAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                                                                        >  2:141149/1‑150 (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:1636227/150‑1 (MQ=255)
                                                                ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                     >  2:475160/1‑150 (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                                       <  2:1454633/150‑1 (MQ=255)
                                                                                                              ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCa                                       <  1:1388968/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:586135/150‑1 (MQ=255)
                                                                                                                                       ctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggc              <  2:594037/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)
                                                                                                                                                   |                                                                                                                                                     
GCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGG  >  NC_000913/1602089‑1602385

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: