Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,631,642 G→C 6.4% D304H (GAC→CAC) ‡ slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,6420GC 6.4% 317.2 / 9.6 93D304H (GAC→CAC) ‡sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (46/41);  new base C (2/4);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCAGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCG  >  NC_000913/4631494‑4631791
                                                                                                                                                    |                                                                                                                                                     
gCAGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGa                                                                                                                                                      >  1:995388/1‑150 (MQ=255)
 cAGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGAc                                                                                                                                                     >  1:663681/1‑150 (MQ=255)
  aGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACg                                                                                                                                                    <  1:872921/150‑1 (MQ=255)
   ggATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGc                                                                                                                                                   >  2:1524777/1‑150 (MQ=255)
    gATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGcc                                                                                                                                                  <  1:441412/150‑1 (MQ=255)
      tGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCAt                                                                                                                                                >  2:1092140/1‑150 (MQ=255)
       gCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCAtt                                                                                                                                               <  1:1313939/150‑1 (MQ=255)
        cTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTa                                                                                                                                              >  2:122584/1‑150 (MQ=255)
         tGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGTCATCGCGCCAtttc                                                                                                                                             <  1:869339/150‑3 (MQ=255)
         tGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTAt                                                                                                                                             >  1:1355824/1‑150 (MQ=255)
          gagaATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGTCATCGCGCCAtttcg                                                                                                                                            <  2:1310231/150‑4 (MQ=255)
          gagaATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGTCATCGCGCCAtttcg                                                                                                                                            <  2:1098894/150‑4 (MQ=255)
          gagaATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATg                                                                                                                                            <  2:394050/150‑1 (MQ=255)
            gaATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGTCATCGCGCCAtttcgcc                                                                                                                                          <  2:663388/150‑6 (MQ=255)
               tGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTc                                                                                                                                       >  2:84187/1‑150 (MQ=255)
                 cACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGc                                                                                                                                     >  1:228721/1‑150 (MQ=255)
                  aCGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCa                                                                                                                                    <  1:158249/150‑1 (MQ=255)
                      cTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGGCGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCt                                                                                                                                <  1:1631770/150‑1 (MQ=255)
                              cccATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCt                                                                                                                        >  2:1254349/1‑150 (MQ=255)
                              cccATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCt                                                                                                                        <  2:1447694/150‑1 (MQ=255)
                              cccATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCt                                                                                                                        <  2:867255/150‑1 (MQ=255)
                                cATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTa                                                                                                                      <  2:1251966/150‑1 (MQ=255)
                                         cccAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGcc                                                                                                             >  1:673158/1‑150 (MQ=255)
                                         cccAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGAGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGACAAATCGCGCCATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGcc                                                                                                             >  2:607058/1‑150 (MQ=255)
                                          ccAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCg                                                                                                            <  1:554893/150‑1 (MQ=255)
                                          ccAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCg                                                                                                            >  1:1006141/1‑150 (MQ=255)
                                          ccAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCg                                                                                                            <  1:526164/150‑1 (MQ=255)
                                           cAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCgt                                                                                                           >  2:462635/1‑150 (MQ=255)
                                            aGGCGCAGCAGCTTAATGAAGATCAGGTTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCgtg                                                                                                          >  2:1578367/1‑150 (MQ=255)
                                            aGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCgtg                                                                                                          <  1:462635/150‑1 (MQ=255)
                                                gcagcaGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACg                                                                                                      >  1:723563/1‑150 (MQ=255)
                                                      gcTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGgc                                                                                                >  2:492782/1‑150 (MQ=255)
                                                        ttAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGgcgc                                                                                              >  2:926481/1‑150 (MQ=255)
                                                         tAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCACCGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGACCGCCGTGTACGAATGGCGCt                                                                                             >  2:1527367/1‑150 (MQ=255)
                                                         tAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCt                                                                                             <  1:926481/150‑1 (MQ=255)
                                                           aTGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTg                                                                                           <  1:1481194/150‑1 (MQ=255)
                                                            tGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAAc                                                                                                                >  1:416534/1‑128 (MQ=255)
                                                            tGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAAc                                                                                                                <  2:416534/128‑1 (MQ=255)
                                                            tGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTgg                                                                                          <  2:1225170/150‑1 (MQ=255)
                                                            tGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTgg                                                                                          >  2:839427/1‑150 (MQ=255)
                                                             gAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGc                                                                                         >  2:705173/1‑150 (MQ=255)
                                                             gAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGc                                                                                         <  2:995388/150‑1 (MQ=255)
                                                               aaaTCACATACAGTAGCTGAACCATATCGTAGACTGTCGTTTGATGTGCAACGATTTAACCTACCACCAACCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCAc                                                                                       >  2:661992/3‑150 (MQ=255)
                                                                gATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCAcc                                                                                      >  2:268029/1‑150 (MQ=255)
                                                                         cAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGAtcgt                                                                             >  2:764100/1‑150 (MQ=255)
                                                                           ggAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGAtcgtcg                                                                           >  1:1608078/1‑150 (MQ=255)
                                                                            gAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTcgc                                                                          <  1:222546/150‑1 (MQ=255)
                                                                             aGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTcgcg                                                                         <  2:1320961/150‑1 (MQ=255)
                                                                                tGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGcc                                                                      >  2:225079/1‑150 (MQ=255)
                                                                                tGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGcc                                                                      >  1:91374/1‑150 (MQ=255)
                                                                                tGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGcc                                                                      <  1:764100/150‑1 (MQ=255)
                                                                                 gcgcGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCt                                                                     >  1:1136411/1‑150 (MQ=255)
                                                                                    cgATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGaa                                                                  <  2:323614/150‑1 (MQ=255)
                                                                                     gATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAAt                                                                 <  1:1092140/150‑1 (MQ=255)
                                                                                     gATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAAt                                                                 <  1:1035579/150‑1 (MQ=255)
                                                                                     gATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAAt                                                                 >  2:29944/1‑150 (MQ=255)
                                                                                         tcgtcgCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACct                                                             <  1:122584/150‑1 (MQ=255)
                                                                                             cgCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTTAATACctggct                                                         <  2:1006141/150‑1 (MQ=255)
                                                                                              gCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTcgcg                                                                         <  1:39032/133‑1 (MQ=255)
                                                                                              gCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTcgcg                                                                         >  2:39032/1‑133 (MQ=255)
                                                                                              gCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACctggctg                                                        <  2:663681/150‑1 (MQ=255)
                                                                                                  ggCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcgc                                                    >  1:864803/1‑150 (MQ=255)
                                                                                                   gCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                   <  2:276038/150‑1 (MQ=255)
                                                                                                   gCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                   <  1:1524777/150‑1 (MQ=255)
                                                                                                     gTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTc                                                 >  2:494027/1‑150 (MQ=255)
                                                                                                       ttGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                               >  1:1477708/1‑150 (MQ=255)
                                                                                                        tGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                               <  1:101138/149‑1 (MQ=255)
                                                                                                        tGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                               >  2:95427/1‑149 (MQ=255)
                                                                                                        tGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                               <  1:95427/149‑1 (MQ=255)
                                                                                                        tGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                               >  2:101138/1‑149 (MQ=255)
                                                                                                         gATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                             >  1:125270/1‑150 (MQ=255)
                                                                                                             ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATg                                         <  1:1254349/150‑1 (MQ=255)
                                                                                                              gCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                        >  1:682572/1‑150 (MQ=255)
                                                                                                                  cGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGc                                    <  2:723563/150‑1 (MQ=255)
                                                                                                                        cACCGACGAGGAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTGGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGGGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGTAGCgtaaga                              <  2:633658/150‑1 (MQ=255)
                                                                                                                        cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaaga                              >  2:446131/1‑150 (MQ=255)
                                                                                                                         aCCGACGAGCAGGCGAAATGGCGCGATGACGCGATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagag                             >  1:550807/1‑150 (MQ=255)
                                                                                                                             acgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                         <  2:1534406/150‑1 (MQ=255)
                                                                                                                             acgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                         <  1:54635/150‑1 (MQ=255)
                                                                                                                                aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATg                      >  2:658789/1‑150 (MQ=255)
                                                                                                                                aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATg                      <  2:1355824/150‑1 (MQ=255)
                                                                                                                                 gCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGa                     >  2:86584/1‑150 (MQ=255)
                                                                                                                                   gcgcataaTGGCGTCATCGCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAAt                   >  1:895672/7‑150 (MQ=255)
                                                                                                                                   aGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAAt                   >  2:211545/1‑150 (MQ=255)
                                                                                                                                    ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATg                  <  2:914179/150‑1 (MQ=255)
                                                                                                                                     gCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATgg                 >  2:743475/1‑150 (MQ=255)
                                                                                                                                      cataaTGGCGTCATCGCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGc                >  1:1173386/4‑150 (MQ=255)
                                                                                                                                      cGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGc                >  2:858189/1‑150 (MQ=255)
                                                                                                                                      cGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGc                <  2:91374/150‑1 (MQ=255)
                                                                                                                                              cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTgg        <  2:82592/150‑1 (MQ=255)
                                                                                                                                              cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTgg        >  1:275666/1‑150 (MQ=255)
                                                                                                                                                cgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGTATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCa      <  1:492782/150‑1 (MQ=255)
                                                                                                                                                   tGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGc   >  2:628503/1‑150 (MQ=255)
                                                                                                                                                    gACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGTGTTACTGGCAGGCg  <  1:446131/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTCAGGAGCTGCGCGATATCGTCGCCTGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCG  >  NC_000913/4631494‑4631791

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: