Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,236 A→T 11.2% Q249H (CAA→CAT 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,2360AT11.2% 192.8 / 11.8 69Q249H (CAA→CATlsrAAutoinducer‑2 ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base T (5/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTGCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGG  >  NC_000913/1602087‑1602385
                                                                                                                                                     |                                                                                                                                                     
cTGCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaa                                                                                                                                                       >  1:1011984/1‑150 (MQ=255)
cTGCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaa                                                                                                                                                       >  2:1099104/1‑150 (MQ=255)
    cGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                                                                                                                   <  2:293648/150‑1 (MQ=255)
        aaTTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATg                                                                                                                                               <  1:587722/150‑1 (MQ=255)
          ttCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGc                                                                                                                                             <  2:200944/150‑1 (MQ=255)
             gCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                                                                          <  2:1095218/150‑1 (MQ=255)
                aGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCGGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                                                                       <  2:1059490/150‑4 (MQ=255)
                    tGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                                                                   <  2:210537/150‑1 (MQ=255)
                     gCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCt                                                                                                                                  >  2:1090949/1‑150 (MQ=255)
                       cGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTgg                                                                                                                                <  2:1011984/150‑1 (MQ=255)
                        gATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                                                                                                                               >  1:1223634/1‑150 (MQ=255)
                        gATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                                                                                                                               >  1:1311519/1‑150 (MQ=255)
                            gAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAAcc                                                                                                                           >  2:594352/1‑150 (MQ=255)
                              aTTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCAAGCCATAAATTATGGCTGGAGTTACCTGGTAAcctc                                                                                                                         <  2:588921/150‑3 (MQ=255)
                                tAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGccc                                                                                                                       >  1:111645/1‑150 (MQ=255)
                                tAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGccc                                                                                                                       >  2:213467/1‑150 (MQ=255)
                                    gTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                                                                   >  1:286868/1‑150 (MQ=255)
                                      gATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca                                                                                                                 <  1:101547/150‑1 (MQ=255)
                                        tGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg                                                                                                               >  1:391336/1‑150 (MQ=255)
                                           gcgACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccg                                                                                                            >  2:1255815/1‑150 (MQ=255)
                                              aCGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCgg                                                                                                        <  1:1099104/150‑1 (MQ=255)
                                                 gAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGaa                                                                                                      <  1:819645/150‑1 (MQ=255)
                                                     cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGcc                                                                                                  <  1:768891/150‑1 (MQ=255)
                                                         gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTg                                                                                              <  2:286868/150‑1 (MQ=255)
                                                          ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGc                                                                                             <  2:400727/150‑1 (MQ=255)
                                                            ttAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCATCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTg                                                                                           <  2:111645/150‑1 (MQ=255)
                                                                gCGGCAAAAACAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGacac                                                                                       >  2:949986/1‑150 (MQ=255)
                                                                 cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACt                                                                                      >  2:998017/1‑150 (MQ=255)
                                                                  ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTg                                                                                     >  2:845042/1‑150 (MQ=255)
                                                                   gCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTgg                                                                                    <  2:70783/150‑1 (MQ=255)
                                                                    cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCt                                                                                            >  2:1168689/1‑141 (MQ=255)
                                                                    cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCt                                                                                            <  1:1168689/141‑1 (MQ=255)
                                                                      aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaaa                                                                                 >  2:527910/1‑150 (MQ=255)
                                                                       aaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaaaa                                                                                <  1:594352/150‑1 (MQ=255)
                                                                         ccAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                                                                              <  1:638029/150‑1 (MQ=255)
                                                                         ccAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                                                                              <  2:1311519/150‑1 (MQ=255)
                                                                         ccAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                                                                              <  2:692617/150‑1 (MQ=255)
                                                                               aaCTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCg                                                                        >  2:1144209/1‑150 (MQ=255)
                                                                                       aCCCACGACATTATTCACGCCATCACCCCAGCGCTACGTGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGCTGCTGACACTGGAAAATCTGACCGGCGAAGGt                                                                >  2:22027/1‑150 (MQ=255)
                                                                                               cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa                                                        >  2:371843/1‑150 (MQ=255)
                                                                                               cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa                                                        >  1:960011/1‑150 (MQ=255)
                                                                                                     tCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAg                                                  <  2:1192545/150‑1 (MQ=255)
                                                                                                      cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGAGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGc                                                 >  2:758221/1‑150 (MQ=255)
                                                                                                       aGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAggc                                                <  1:949986/150‑3 (MQ=255)
                                                                                                        ggCCATCCCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTCGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCt                                               >  2:138884/1‑150 (MQ=255)
                                                                                                          ccATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGa                                             <  1:845042/150‑1 (MQ=255)
                                                                                                              cACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGc                                                                                             >  2:1311589/1‑98 (MQ=255)
                                                                                                              cACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGc                                                                                             <  1:1311589/98‑1 (MQ=255)
                                                                                                                ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCa                                       <  2:1399692/150‑1 (MQ=255)
                                                                                                                ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCa                                       >  1:633491/1‑150 (MQ=255)
                                                                                                                 cccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCaa                                      >  1:1025285/1‑150 (MQ=255)
                                                                                                                   cAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATg                                    >  2:1220071/1‑150 (MQ=255)
                                                                                                                    aGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGc                                   >  2:474188/1‑150 (MQ=255)
                                                                                                                           cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGaa                            <  1:1297502/150‑1 (MQ=255)
                                                                                                                           cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGaa                            >  1:401807/1‑150 (MQ=255)
                                                                                                                            gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGaaa                           >  1:603843/1‑150 (MQ=255)
                                                                                                                            gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGaaa                           >  1:92392/1‑150 (MQ=255)
                                                                                                                                    tCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGc                   <  1:1317089/150‑1 (MQ=255)
                                                                                                                                          tctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggct             <  2:1223634/150‑1 (MQ=255)
                                                                                                                                           actccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggctg            >  1:846150/4‑150 (MQ=255)
                                                                                                                                            tGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggctgg           >  2:237991/1‑150 (MQ=255)
                                                                                                                                            ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggc              >  1:309135/3‑147 (MQ=255)
                                                                                                                                            ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCctggc              <  2:309135/145‑1 (MQ=255)
                                                                                                                                             tccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctggg          >  1:61471/2‑150 (MQ=255)
                                                                                                                                             gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGcc                  <  1:257002/142‑1 (MQ=255)
                                                                                                                                             gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctggg          >  1:74995/1‑150 (MQ=255)
                                                                                                                                             gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACAGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGcc                  >  2:257002/1‑142 (MQ=255)
                                                                                                                                              ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctgggc         >  1:367094/1‑150 (MQ=255)
                                                                                                                                              ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTggctgggc         >  2:1420210/1‑150 (MQ=255)
                                                                                                                                                 gccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGgctgggctgg      >  2:646604/1‑150 (MQ=255)
                                                                                                                                                   caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCtggtg    <  2:391336/150‑1 (MQ=255)
                                                                                                                                                     aaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTggg  <  1:646604/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCCGGAAATTCGCCAGATTGCCGATCGAATTAGCGTGATGCGCGACGGAACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAGCCTGACGCTCAATGCCGGAGAAATTCTGGGCCTGGCTGGGCTGGTGGG  >  NC_000913/1602087‑1602385

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: