Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,155 A→G 6.8% pseudogene (1124/1358 nt) gatC ← galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,1550AG 6.8% 257.8 / 6.8 88pseudogene (1124/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base G (2/4);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGG  >  NC_000913/2173006‑2173304
                                                                                                                                                     |                                                                                                                                                     
ccGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATa                                                                                                                                                       >  1:119920/1‑150 (MQ=255)
ccGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATa                                                                                                                                                       <  2:1082026/150‑1 (MQ=255)
 cGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAg                                                                                                                                                      <  2:1613149/150‑1 (MQ=255)
   ataataATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTATTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGtt                                                                                                                                                    >  1:363275/1‑150 (MQ=255)
     aataatTAAAACGCGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGTTAATGTGAGAACCGCCCTGATCAATTCAAGCCACCATAACAACGTCTATTAAATCAACAGCATTATCCGCCAGTTGGCTGTGCAAACAAATAGttta                                                                                                                                                  >  2:1830846/1‑149 (MQ=255)
      ataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGc                                                                                                                                                 <  2:669146/150‑1 (MQ=255)
       taatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCCGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCg                                                                                                                                                <  1:1065263/150‑1 (MQ=255)
        aatGAAACCGGGAATATTTTGCGGGGAGAAACCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcac                                                                                                                                               <  1:1899311/150‑3 (MQ=255)
        aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGt                                                                                                                                               >  2:318841/1‑150 (MQ=255)
        aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGt                                                                                                                                               <  1:298794/150‑1 (MQ=255)
         attaaaaCGGGAATATTTTGCGGGGAGAAAACCTGAAAAAGTCACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTc                                                                                                                                              >  2:30135/4‑150 (MQ=255)
         atGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcaca                                                                                                                                              <  1:83539/150‑4 (MQ=255)
          tGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcacac                                                                                                                                             <  1:315366/150‑5 (MQ=255)
           gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcacacc                                                                                                                                            <  1:1698568/150‑6 (MQ=255)
           gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgcg                                                                                                                                           <  2:829049/150‑1 (MQ=255)
            aaaCCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgcg                                                                                                                                           >  1:911575/1‑150 (MQ=255)
              aCCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGAt                                                                                                                                         <  1:385802/150‑1 (MQ=255)
                  ggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCcac                                                                                                                                     >  2:949851/1‑150 (MQ=255)
                  ggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCcac                                                                                                                                     >  2:1983188/1‑150 (MQ=255)
                           tGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATg                                                                                                                            >  1:92855/1‑150 (MQ=255)
                              ggggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTc                                                                                                                         >  2:1584901/1‑150 (MQ=255)
                               aggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCa                                                                                                                        >  2:1907105/2‑150 (MQ=255)
                                ggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCAt                                                                                                                       >  2:1045289/1‑150 (MQ=255)
                                    aaaaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatg                                                                                                                   <  2:209297/150‑1 (MQ=255)
                                      aaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgat                                                                                                                 >  2:1489804/1‑150 (MQ=255)
                                      aaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgat                                                                                                                 <  2:102814/150‑1 (MQ=255)
                                      aaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgat                                                                                                                 <  1:407968/150‑1 (MQ=255)
                                       aCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                                >  1:1054549/1‑150 (MQ=255)
                                       aCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                                <  2:1808905/150‑1 (MQ=255)
                                       aCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                                <  1:628590/150‑1 (MQ=255)
                                         cTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGac                                                                                                              >  2:1411156/1‑150 (MQ=255)
                                           gAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAACATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGacac                                                                                                            >  2:1182705/1‑150 (MQ=255)
                                                aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCgag                                                                                                       <  1:1830846/150‑1 (MQ=255)
                                                aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCgac                                                                                                       >  1:222197/1‑149 (MQ=255)
                                                 gTAACCAGGTAATGGGAGAACCGCCCTGATCAATTGAAGCCACCATACCCCCGGCTTTTAACGCACCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCgaga                                                                                                      >  2:1509407/1‑150 (MQ=255)
                                                  tAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGAt                                                                                                     <  1:1489804/150‑1 (MQ=255)
                                                    aCCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTa                                                                                                   >  2:1608162/1‑150 (MQ=255)
                                                     ccAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTaa                                                                                                  >  2:1853021/1‑150 (MQ=255)
                                                      cAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAg                                                                                                 <  1:594673/150‑1 (MQ=255)
                                                         gTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTg                                                                                              <  2:1054549/150‑1 (MQ=255)
                                                               ggAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAAc                                                                                        <  1:1435450/150‑1 (MQ=255)
                                                                gagaACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACa                                                                                       <  1:318841/150‑1 (MQ=255)
                                                                 agaACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAg                                                                                      <  2:627139/150‑1 (MQ=255)
                                                                    aCCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGAtt                                                                                   >  2:1958909/1‑150 (MQ=255)
                                                                       gCCCTGATCCATTGAAGCCACCATACCCCGTGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGGCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTcc                                                                                <  2:363275/150‑1 (MQ=255)
                                                                         ccTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCAc                                                                              <  1:1584901/150‑1 (MQ=255)
                                                                         ccTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCAc                                                                              <  1:235673/150‑1 (MQ=255)
                                                                            gATCCATTGAAGCCACCATACCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgcg                                                                                                                                           >  1:769524/1‑85 (MQ=255)
                                                                            gATCCATTGAAGCCACCATACCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgcg                                                                                                                                           <  2:769524/85‑1 (MQ=255)
                                                                              tCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGc                                                                         <  2:1118864/150‑1 (MQ=255)
                                                                              tCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGc                                                                         <  2:402005/150‑1 (MQ=255)
                                                                              tCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGc                                                                         >  1:402005/1‑150 (MQ=255)
                                                                              tCCATTGAAGCCACCATACCCCCGGATTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGc                                                                         >  1:1601198/1‑150 (MQ=255)
                                                                                    gaagaaaccaTACCCCCTGCTTTTAACGCCCCAGCATTAGCCGCCACTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGCGTGATGCTCATAATGATGACACCCGACATAAAGGTGCGGAACAGATTTCAACGATGCACGGcc                                                                   >  2:1409512/7‑149 (MQ=255)
                                                                                    gAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCg                                                                   >  2:1956985/1‑150 (MQ=255)
                                                                                       gCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGAcc                                                                <  1:30135/150‑1 (MQ=255)
                                                                                       gCCACCATACCCCCGGCTTTTAACGCCACAGCATTAGCCGACAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGAcc                                                                >  2:967596/1‑150 (MQ=255)
                                                                                          accaTACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGcc                                                             >  2:1430234/1‑150 (MQ=255)
                                                                                          accaTACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGcc                                                             <  2:1639958/150‑1 (MQ=255)
                                                                                           ccaTACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCa                                                            >  2:231826/1‑150 (MQ=255)
                                                                                              tACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCg                                                         <  2:198392/150‑1 (MQ=255)
                                                                                                  cccGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCAc                                                     <  1:1290038/150‑1 (MQ=255)
                                                                                                    cGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACga                                                   <  2:92855/150‑1 (MQ=255)
                                                                                                          ttAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGc                                             <  1:1907105/150‑1 (MQ=255)
                                                                                                           tAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGcc                                            >  1:1775998/1‑150 (MQ=255)
                                                                                                           tAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGcc                                            <  1:1983188/150‑1 (MQ=255)
                                                                                                            aaCGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCg                                           >  2:698224/1‑150 (MQ=255)
                                                                                                                ccccAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAtgg                                       >  2:1047392/1‑150 (MQ=255)
                                                                                                                 cccAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAtggt                                      <  2:445869/150‑1 (MQ=255)
                                                                                                                    aGCATTAGCCGCTAGTTGGGTGTGCAAACCAATAGTTTGCGTTGTTATTAACTGGTTGATGCTCATAATGTTTATACCCGAGATTAATGTGTGTAACAGATTTCCACGATTCAGTGCGACCTCCATCGCCACGAAGAAGCCGATGGTGGc                                   <  2:222197/150‑1 (MQ=255)
                                                                                                                    aGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGc                                   >  1:1489672/1‑150 (MQ=255)
                                                                                                                     gCATTAGCCGCCAGTTGGGTGTGCAAACCACTAGTTTGCGTCGCGATCCACAGGGGGATGCTCATAATCTTGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGTAGAAGCCGATGGTGGCa                                  <  2:885158/150‑1 (MQ=255)
                                                                                                                     gCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCa                                  >  1:1484744/1‑150 (MQ=255)
                                                                                                                           gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATTATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCg                            >  1:932740/1‑150 (MQ=255)
                                                                                                                           gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCg                            <  1:1958909/150‑1 (MQ=255)
                                                                                                                           gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCg                            <  2:1775998/150‑1 (MQ=255)
                                                                                                                             cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGcc                          <  1:1608162/150‑1 (MQ=255)
                                                                                                                                 aGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAc                      >  2:1234029/1‑150 (MQ=255)
                                                                                                                                    tGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgc                   <  1:1045289/150‑1 (MQ=255)
                                                                                                                                        tgtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagca               <  1:233160/150‑1 (MQ=255)
                                                                                                                                        gacgCAAACTATTGGTTTGCGTCGCGATCCATAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagca               >  2:160884/4‑150 (MQ=255)
                                                                                                                                         gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCAc              <  1:235671/150‑1 (MQ=255)
                                                                                                                                         acgCAAACTATTGGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCAc              >  2:460550/3‑150 (MQ=255)
                                                                                                                                           gCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCt            >  1:1438940/1‑150 (MQ=255)
                                                                                                                                           gCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCt            <  1:1411156/150‑1 (MQ=255)
                                                                                                                                            cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTg           <  1:698224/150‑1 (MQ=255)
                                                                                                                                            cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTg           >  1:919117/1‑150 (MQ=255)
                                                                                                                                              aaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGAt         <  2:919117/150‑1 (MQ=255)
                                                                                                                                                 cAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTAc      >  2:894899/1‑150 (MQ=255)
                                                                                                                                                   aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTAccc    <  1:20922/150‑1 (MQ=255)
                                                                                                                                                     aGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGACAAGATCGCCAAACGGCAGCACCTGATTACACgg  >  2:873441/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGG  >  NC_000913/2173006‑2173304

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