Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,363,452 G→A 5.1% A61V (GCC→GTC) ‡ yfaY ← inactive PncC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,363,4520GA 5.1% 340.8 / 5.5 98A61V (GCC→GTC) ‡yfaYinactive PncC family protein
Reads supporting (aligned to +/- strand):  ref base G (39/54);  new base A (2/3);  total (41/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAACACT  >  NC_000913/2363303‑2363598
                                                                                                                                                     |                                                                                                                                                  
ttACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCgg                                                                                                                                                    >  2:1114607/1‑150 (MQ=255)
ttACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCgg                                                                                                                                                    >  2:1114639/1‑150 (MQ=255)
 tACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGc                                                                                                                                                   <  1:1026155/150‑1 (MQ=255)
    tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                                >  1:994413/1‑150 (MQ=255)
    tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                                <  1:854454/150‑1 (MQ=255)
    tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                                >  2:1399400/1‑150 (MQ=255)
    tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                                >  1:1241798/1‑150 (MQ=255)
      gTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCt                                                                                                                                              <  1:143557/150‑1 (MQ=255)
        ccACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGa                                                                                                                                            <  2:823247/150‑1 (MQ=255)
         cACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGAc                                                                                                                                           >  1:665800/1‑150 (MQ=255)
          aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatc                                                                                                                                          <  1:1168457/150‑3 (MQ=255)
          aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatc                                                                                                                                          <  1:1168460/150‑3 (MQ=255)
          aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACt                                                                                                                                          <  2:860210/150‑1 (MQ=255)
              tCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatcgtta                                                                                                                                      <  2:348882/150‑7 (MQ=255)
                gTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTc                                                                                                                                    >  2:350061/1‑150 (MQ=255)
                gTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTc                                                                                                                                    >  2:582752/1‑150 (MQ=255)
                   gAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTcgcg                                                                                                                                 <  1:741764/150‑1 (MQ=255)
                   gAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTcgcg                                                                                                                                 <  1:565586/150‑1 (MQ=255)
                        agggTTCCATCTCTTTGAGGGAGGCTTCTTGCAGCACCCGGCGGTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCCTCGCTGGTCGGTCCCAGCCCGCCGTTAACGCTCTGCACATCGGCGTGGTGACTACGTTCGCgaagaa                                                                                                                            <  2:1235182/147‑1 (MQ=255)
                        aGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGGCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCgaagaa                                                                                                                            <  2:1235181/150‑1 (MQ=255)
                          cGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATg                                                                                                                          >  1:1168967/1‑150 (MQ=255)
                            ttCCATCTCTTTGTGACTGGCTTCTTGCTGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGTGTGCGCTTAAATAATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTTCTGACTACGTTCGCGAAGAATGGt                                                                                                                        <  2:1392887/150‑1 (MQ=255)
                            ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                        <  1:825332/150‑1 (MQ=255)
                            ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                        <  2:11674/150‑1 (MQ=255)
                            ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                        <  2:520665/150‑1 (MQ=255)
                             tCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTg                                                                                                                       >  2:731164/1‑150 (MQ=255)
                             tCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTg                                                                                                                       >  2:379870/1‑150 (MQ=255)
                              ccATCTCTTTGAGCCAGGCTTCATGCAGCACCCGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGa                                                                                                                      <  2:915593/150‑1 (MQ=255)
                              ccATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGa                                                                                                                      <  2:590028/150‑1 (MQ=255)
                              ccATCTCTTTGAGCCAGGCGTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTCAATCATCGCTGGTCGTTCCCATCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGTAGAATGGTGa                                                                                                                      <  2:775599/150‑1 (MQ=255)
                                 tctcTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTa                                                                                                                   >  1:527910/1‑150 (MQ=255)
                                     tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                               >  2:808440/1‑150 (MQ=255)
                                     tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                               <  2:567138/150‑1 (MQ=255)
                                     tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                               >  1:720437/1‑150 (MQ=255)
                                      ttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtca                                                                                                              <  2:34870/150‑1 (MQ=255)
                                        gcgCCAGGCTTCCTGCCGCCCCCGGCCTGCACCTCGCGCTGGCGCAGCGGCGAGTGCGCATAACCCATCGCTGGACGGTCCCAGCCCGCCGAGACCGACCAGCACAGCGGCGTGCTGACAACGTTCGCGAAGAATGGTGACTAAGtcatc                                                                                                            <  2:1052265/148‑1 (MQ=255)
                                         aGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtcatca                                                                                                           >  1:1367445/1‑150 (MQ=255)
                                           ccAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAg                                                                                                         <  1:855756/150‑1 (MQ=255)
                                            cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                        >  2:448726/1‑150 (MQ=255)
                                            cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                        >  2:1133360/1‑150 (MQ=255)
                                            cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                        >  2:104605/1‑150 (MQ=255)
                                            cAGGCTTCATGCAGCACCACGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                        <  1:396888/150‑1 (MQ=255)
                                             aGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGt                                                                                                       >  1:442416/1‑150 (MQ=255)
                                             aGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGt                                                                                                       >  1:442380/1‑150 (MQ=255)
                                                 ttCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATc                                                                                                   <  1:299907/150‑1 (MQ=255)
                                                   cATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATccc                                                                                                 <  1:135384/150‑1 (MQ=255)
                                                   cATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATccc                                                                                                 <  1:135321/150‑1 (MQ=255)
                                                      gcagcaCCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCa                                                                                              <  2:51981/150‑1 (MQ=255)
                                                              aGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTg                                                                                      >  1:938893/1‑150 (MQ=255)
                                                                gCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcg                                                                                    >  2:1189958/1‑150 (MQ=255)
                                                                 ccTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcgg                                                                                   <  1:1160985/150‑1 (MQ=255)
                                                                 ccTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcgg                                                                                   <  1:1160931/150‑1 (MQ=255)
                                                                  cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  >  2:848208/1‑150 (MQ=255)
                                                                  cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  >  1:1042251/1‑150 (MQ=255)
                                                                  cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  <  1:339699/150‑1 (MQ=255)
                                                                  cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  <  1:504653/150‑1 (MQ=255)
                                                                  cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                  >  1:584008/1‑150 (MQ=255)
                                                                   ttaaCCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGACGTTAACGATCAGCACATCGGCGTGCTGACTCCGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgaggcg                                                                                 >  1:1266772/4‑150 (MQ=255)
                                                                   ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                 >  1:528212/1‑150 (MQ=255)
                                                                   ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                 >  2:225994/1‑150 (MQ=255)
                                                                   ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                 >  1:1266775/1‑150 (MQ=255)
                                                                          tttGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAAt                                                                          <  1:797783/150‑1 (MQ=255)
                                                                                gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                    <  1:641663/150‑1 (MQ=255)
                                                                                gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                    <  1:641594/150‑1 (MQ=255)
                                                                                gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                    <  1:641575/150‑1 (MQ=255)
                                                                                    cAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCt                                                                <  1:509872/150‑1 (MQ=255)
                                                                                    cAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCt                                                                <  1:509913/150‑1 (MQ=255)
                                                                                      gcggcgAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCtga                                                              >  1:974584/1‑150 (MQ=255)
                                                                                      gcggcgAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCtga                                                              <  2:448035/150‑1 (MQ=255)
                                                                                           gAGTGCGCTTATATCATCGCTGGTCGGTCCCAGCGCGCCGTTAACGATCAGCTCATCGGCGTGCTGACTACGTTCGCGAATAATGGTGACTAAGTCATCTAGGTTATCCCCCACCTTATTGCGGCGATATAATGGCAACCCCTGATgaaa                                                         <  2:793551/150‑1 (MQ=255)
                                                                                           gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  1:574012/150‑1 (MQ=255)
                                                                                           gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  1:1226747/150‑1 (MQ=255)
                                                                                           gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  2:1129755/150‑1 (MQ=255)
                                                                                           gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  1:874063/150‑1 (MQ=255)
                                                                                           gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                         <  1:236516/150‑1 (MQ=255)
                                                                                                    tAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATc                                                <  1:1206813/150‑1 (MQ=255)
                                                                                                     aaaTCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCg                                               <  2:759619/150‑1 (MQ=255)
                                                                                                      aaTCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCgg                                              <  2:1313329/150‑1 (MQ=255)
                                                                                                          atcGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGgccag                                          >  1:1258251/1‑150 (MQ=255)
                                                                                                             gCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGgccagcca                                       <  2:1128063/150‑1 (MQ=255)
                                                                                                                  tcggtcCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAg                                  >  2:227570/1‑150 (MQ=255)
                                                                                                                  tcggtcCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAg                                  >  2:1124956/1‑150 (MQ=255)
                                                                                                                    ggtcCCTGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTAGGTTCGCTATTATTGGTGACTAAGTCATCAAGGTTATCCCCCACCTTTTTGCGGCGATTTAATGGCAACCCCTTATGAAATAAAAAATCGGCCAGCCAGGCAGCg                                <  1:897027/150‑1 (MQ=255)
                                                                                                                       cccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTa                             >  2:731377/1‑150 (MQ=255)
                                                                                                                        ccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAg                            >  1:1236040/1‑150 (MQ=255)
                                                                                                                        ccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATAGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAg                            >  1:1236050/1‑150 (MQ=255)
                                                                                                                         cAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAgt                           <  2:146587/150‑1 (MQ=255)
                                                                                                                            cccGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGACAAAATCGGCCAGCCAGGCAGCGTTAGTGTc                        >  2:729875/1‑150 (MQ=255)
                                                                                                                                   ttAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGagt                 >  2:1060865/1‑148 (MQ=255)
                                                                                                                                     aaCGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTg               >  1:748539/1‑150 (MQ=255)
                                                                                                                                       tagtCAGCACGCCGATGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGcc             >  2:597425/4‑150 (MQ=255)
                                                                                                                                         gtCAGCACGCCGATGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCg           >  1:498339/2‑150 (MQ=255)
                                                                                                                                         aTCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCTTATTGCGGCGAGATATTGGAAACCCATGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCg           <  1:416598/150‑1 (MQ=255)
                                                                                                                                          tccgCACATGGGCGTGCTGAGTCGGTTCGCGACGAAGGGTGACTAAGTCATCTCGGTTATCCCCCACGGTATTGGGGCGTGTTAATGGCAACCCCTTATTAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCgt          <  1:551194/147‑1 (MQ=255)
                                                                                                                                             gCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTa       >  1:1375087/1‑150 (MQ=255)
                                                                                                                                              cacaTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTaa      >  1:1144033/1‑150 (MQ=255)
                                                                                                                                                caTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAaca    <  2:183249/150‑1 (MQ=255)
                                                                                                                                                 aTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAacac   <  2:1291705/150‑1 (MQ=255)
                                                                                                                                                 aTCGGCGTGCTGACGACGTTCGCGAAGACTGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCTACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGTTTTGCACGTGTAacac   <  2:1291684/150‑1 (MQ=255)
                                                                                                                                                  tCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAACACt  <  2:731528/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAACACT  >  NC_000913/2363303‑2363598

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