Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,242 A→T 8.4% L251F (TTA→TTT lsrA → Autoinducer‑2 ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,602,2420AT 8.4% 240.6 / 10.9 83L251F (TTA→TTTlsrAAutoinducer‑2 ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (49/27);  new base T (2/5);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

GCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGGGG  >  NC_000913/1602100‑1602387
                                                                                                                                              |                                                                                                                                                 
gCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                                                                            <  2:592845/150‑1 (MQ=255)
 ccAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCTGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                                                                           <  2:1568371/150‑2 (MQ=255)
 ccAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                                                                           <  2:1083089/150‑2 (MQ=255)
 ccAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                                                                           <  1:927519/150‑2 (MQ=255)
   aGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                                                                         <  2:27690/150‑4 (MQ=255)
   aGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                                                                         <  1:169845/150‑4 (MQ=255)
   aGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGt                                                                                                                                         <  2:837111/150‑1 (MQ=255)
     aTTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTa                                                                                                                                       <  1:1602326/150‑1 (MQ=255)
       tGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                                                                     >  1:1455996/1‑150 (MQ=255)
        gCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCt                                                                                                                                    <  2:813095/150‑1 (MQ=255)
          cGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTgg                                                                                                                                  >  2:643736/1‑150 (MQ=255)
           gATCGAATTAGCGTGATGCGCGATGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                                                                                                                                 >  1:1042892/1‑150 (MQ=255)
           gATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                                                                                                                                 <  2:32984/150‑1 (MQ=255)
             tCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCGACCGACGAGATTATTCAGGCCATCACCCGAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCTAAAGTTATGGCTGGAGTTACCTGGTac                                                                                                                               >  2:863608/1‑149 (MQ=255)
                      cGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcaca                                                                                                                      >  2:568112/1‑150 (MQ=255)
                                gACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                                                                            >  2:568281/1‑150 (MQ=255)
                                       ccATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACACCATGCCGCCGGAACGc                                                                                                     >  1:440466/1‑150 (MQ=255)
                                            gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTg                                                                                                >  1:709300/1‑150 (MQ=255)
                                                 aaGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGac                                                                                           >  1:608939/1‑150 (MQ=255)
                                                  aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGaca                                                                                          >  2:220538/1‑150 (MQ=255)
                                                  aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGaca                                                                                          >  1:428398/1‑150 (MQ=255)
                                                  aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGaca                                                                                          >  1:1340898/1‑150 (MQ=255)
                                                    cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACt                                                                                        >  2:1584686/1‑150 (MQ=255)
                                                     ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                       <  2:1042892/150‑1 (MQ=255)
                                                     ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                       <  1:335057/150‑1 (MQ=255)
                                                      gCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTgg                                                                                      >  2:1581368/1‑150 (MQ=255)
                                                      gCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTgg                                                                                      >  1:446621/1‑150 (MQ=255)
                                                       cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGa                                                                                     <  1:568112/150‑1 (MQ=255)
                                                        aaaaCCAGCGAACTGTTTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaa                                                                                    <  1:1523043/150‑1 (MQ=255)
                                                        aaaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaa                                                                                    <  2:1340898/150‑1 (MQ=255)
                                                            ccAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                                                                                <  2:1233963/150‑1 (MQ=255)
                                                            ccAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                                                                                >  1:758816/1‑150 (MQ=255)
                                                             cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                                                                               >  1:642274/1‑150 (MQ=255)
                                                             cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                                                                               <  1:1584686/150‑1 (MQ=255)
                                                              aGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                                                                              <  2:712585/150‑1 (MQ=255)
                                                              aGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                                                                              <  1:1469603/150‑1 (MQ=255)
                                                                 gAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGAcc                                                                           <  1:998559/150‑1 (MQ=255)
                                                                   aCTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCgg                                                                         <  1:568281/150‑1 (MQ=255)
                                                                     tGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACTTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCg                                                                       >  1:549190/1‑150 (MQ=255)
                                                                            cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt                                                                >  1:1154084/1‑150 (MQ=255)
                                                                            cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt                                                                >  1:567553/1‑150 (MQ=255)
                                                                             gacgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACTTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc                                                               <  2:549190/150‑1 (MQ=255)
                                                                               cgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg                                                             <  1:1212073/150‑1 (MQ=255)
                                                                                 acATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGa                                                           >  2:1309752/1‑150 (MQ=255)
                                                                                  cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa                                                          >  1:1444828/1‑150 (MQ=255)
                                                                                    ttattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg                                                        <  2:608939/150‑1 (MQ=255)
                                                                                        tCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAg                                                    >  2:295498/1‑150 (MQ=255)
                                                                                         cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCTAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGc                                                   >  1:1508914/1‑150 (MQ=255)
                                                                                         cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTACCCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGc                                                   >  2:1573141/1‑150 (MQ=255)
                                                                                         cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGc                                                   >  1:151448/1‑150 (MQ=255)
                                                                                         cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGc                                                   >  2:1226553/1‑150 (MQ=255)
                                                                                         cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGc                                                   >  2:1171419/1‑150 (MQ=255)
                                                                                         cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGc                                                   >  1:623508/1‑150 (MQ=255)
                                                                                           ggCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCt                                                 >  1:591332/1‑150 (MQ=255)
                                                                                              cATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCGGAc                                              >  1:196440/1‑150 (MQ=255)
                                                                                                   ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCa                                         >  1:110047/1‑150 (MQ=255)
                                                                                                   ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCa                                         >  1:456630/1‑150 (MQ=255)
                                                                                                    cccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCaa                                        <  1:643736/150‑1 (MQ=255)
                                                                                                      cAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATg                                      <  2:1157187/150‑1 (MQ=255)
                                                                                                          ggTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCgg                                  <  2:424784/150‑1 (MQ=255)
                                                                                                            tACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCCGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGgag                                <  1:1570879/150‑1 (MQ=255)
                                                                                                                 gAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAAtt                           >  2:1624110/1‑150 (MQ=255)
                                                                                                                 gAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAAtt                           >  2:1629541/1‑150 (MQ=255)
                                                                                                                    aaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTg                        <  1:295498/150‑1 (MQ=255)
                                                                                                                    aaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTg                        >  2:128036/1‑150 (MQ=255)
                                                                                                                    aaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTg                        >  1:896816/1‑150 (MQ=255)
                                                                                                                     aaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTgg                       >  2:312294/1‑150 (MQ=255)
                                                                                                                       tCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCTGCCGGAACGCCGCTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGc                     >  2:895714/1‑150 (MQ=255)
                                                                                                                          ctctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctg                  >  2:334079/1‑150 (MQ=255)
                                                                                                                          ctctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctg                  <  2:1290276/150‑1 (MQ=255)
                                                                                                                               ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggctgg             >  1:134049/3‑150 (MQ=255)
                                                                                                                                gccagccaATAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGCCGCTCAATGCCGGAGAAATTCTGGGCCTggctggg            >  2:404735/1‑150 (MQ=255)
                                                                                                                                gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctggg            >  2:224703/1‑150 (MQ=255)
                                                                                                                                 gatctcaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGCAAAATCTGACCGGCGAAGGTTTCAGGAATGTCTGCGTGACGCTCAATGCCGTAGAAATTCTGGGCCTggctgggc           <  1:1545796/145‑1 (MQ=255)
                                                                                                                                 ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctgggc           >  1:1143533/1‑150 (MQ=255)
                                                                                                                                    gccaAAAATTATGGCTGGAGTTACCTGGTAACAGCCCACAACATGGCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGcctgggctgg        >  1:1572051/1‑150 (MQ=255)
                                                                                                                                      caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCtggtg      >  2:509386/1‑150 (MQ=255)
                                                                                                                                      caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCtggtg      >  2:130960/1‑150 (MQ=255)
                                                                                                                                      caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTCACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCCGAGAAATTCTGGGCCTGGCTGGGCtggtg      >  2:833060/1‑150 (MQ=255)
                                                                                                                                      caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTCCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCtggtg      >  2:1258018/1‑150 (MQ=255)
                                                                                                                                      caAAAATTATCGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctgggcagatc      >  2:1545796/1‑145 (MQ=255)
                                                                                                                                       aaaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCtggtgg     >  1:430186/1‑150 (MQ=255)
                                                                                                                                          aaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTggggg  <  2:151448/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGGGG  >  NC_000913/1602100‑1602387

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