Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,336 C→G 11.4% pseudogene (945/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,3360CG11.4% 226.4 / 16.3 85pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (43/31);  new base G (6/4);  total (50/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTG  >  NC_000913/1466192‑1466484
                                                                                                                                                |                                                                                                                                                      
tatCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                                                                                                                   >  1:1243005/1‑150 (MQ=255)
 ttCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTTGCACCGGCACCCAGAtt                                                                                                                                                  <  2:321111/149‑1 (MQ=255)
  tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTa                                                                                                                                                 >  1:906773/1‑150 (MQ=255)
   cACCGGTGACAGCGCGACGGTGGATAACAATGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGCTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                                                                                                                >  2:1571007/1‑150 (MQ=255)
   cACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTaa                                                                                                                                                >  1:1161642/1‑150 (MQ=255)
   cACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTaa                                                                                                                                                <  2:1185778/150‑1 (MQ=255)
     ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAATGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                              <  1:2039725/150‑1 (MQ=255)
      cGGTGACAGCTCGACGTTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGCAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                                                                                                             <  2:1837340/150‑1 (MQ=255)
       ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAACCTCAATAACGAAGGCGAGAGCACCATCACCACCGCTGGCACCGGCACCCAGATTAATGGt                                                                                                                                            >  2:874593/1‑150 (MQ=255)
        gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                           <  1:1626815/150‑1 (MQ=255)
        gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGACAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                           <  2:2040078/150‑1 (MQ=255)
          gACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga‑t                                                                                                                                        >  2:1146517/1‑150 (MQ=255)
            cAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga‑tga                                                                                                                                      >  2:695706/1‑150 (MQ=255)
            cAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCAATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGAACCAGCACCAACGGTGGCACCGGCACCCAGATTAATGGtga‑tga                                                                                                                                      >  2:87488/1‑150 (MQ=255)
             aGCGCGACGGTGGATAACAAGGTCACCATGACCGTCACCGATCAGGATTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGG‑TGAc                                                                                                                                     <  1:329877/150‑1 (MQ=255)
             aGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGAc                                                                                                                                     <  2:73362/150‑1 (MQ=255)
              gcgcGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACg                                                                                                                                    >  2:772461/1‑150 (MQ=255)
                gcgACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGcc                                                                                                                                  >  1:1748073/1‑150 (MQ=255)
                          aTAACAAGGGCACCATGACCGTCACCGGTCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAaac‑a                                                                                                                       <  1:1634790/150‑1 (MQ=255)
                          aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAGGGGTATCCAGAGCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAaac‑a                                                                                                                       <  1:1901554/150‑1 (MQ=255)
                            aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGCCGCCACGGCAaac‑aac                                                                                                                     >  1:1757154/1‑150 (MQ=255)
                            aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAaac‑aac                                                                                                                     >  2:570116/1‑150 (MQ=255)
                                 gggCACCATGACCGTCACCGATCCGGAGTCCATGCGTATCCAGATCGACGGTGACAAGGCCATCGGCAATAACGCAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACgg                                                                                                                <  1:335649/150‑1 (MQ=255)
                                   gCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCa                                                                                                              <  1:984178/150‑1 (MQ=255)
                                   gCACCATGACCGTCACCGATCCGGAGTCCAT‑GGTTTACACATATAAGGTGACAAGGCCATATTAAATAAAGAAGGCGACCGCACCATCCCCAAAGGTTGCCCCCGCACCCATATTAATGGTGATTCACCCCAC‑GCAAACAAACAACgcaa                                                                                                              >  2:1972634/1‑147 (MQ=12)
                                      ccATGGCCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGACAGCACCATCACCAACTGTGGCACCGTCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCaaaa                                                                                                           <  1:321140/150‑1 (MQ=255)
                                         tGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAaccc                                                                                                        >  1:1976136/1‑149 (MQ=255)
                                         tGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAacca                                                                                                        >  2:577961/1‑150 (MQ=255)
                                            ccGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCg                                                                                                     <  1:1624479/150‑1 (MQ=255)
                                               tCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTg                                                                                                  >  2:1347632/1‑150 (MQ=255)
                                                cACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCCAAC‑AACAACGGCAAAACCACCGTTGa                                                                                                 >  1:861112/1‑150 (MQ=255)
                                                 aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTgcc                                                                                                >  1:536533/1‑148 (MQ=255)
                                                  ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACg                                                                                               <  2:1817879/150‑1 (MQ=255)
                                                  ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACg                                                                                               >  1:1453077/1‑150 (MQ=255)
                                                   cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACgg                                                                                              >  2:237139/1‑150 (MQ=255)
                                                    gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGc                                                                                             >  1:898920/1‑150 (MQ=255)
                                                    gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGc                                                                                             <  2:91759/150‑1 (MQ=255)
                                                      tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGCAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCaa                                                                                           >  2:101963/1‑150 (MQ=255)
                                                       ccGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGCTGGCCCCGGCACCCAGCTTAATGGTGA‑TGACGCCCCGGCCAAC‑AACAACGGCAAAACCACCGATGACGGCAAg                                                                                          >  2:302237/1‑150 (MQ=18)
                                                        cGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAgg                                                                                         <  2:1161641/150‑1 (MQ=255)
                                                           aGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                      >  1:619973/1‑150 (MQ=255)
                                                              ccATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                                                   >  2:1673018/1‑150 (MQ=255)
                                                               cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                  >  1:774707/1‑150 (MQ=255)
                                                                aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                 <  2:1614117/150‑1 (MQ=255)
                                                                   ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                                              >  2:1598683/1‑150 (MQ=255)
                                                                    gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGTCAAAACCACCGTTGACGGCAAGGATTCCACCggga                                                                             >  1:978736/1‑148 (MQ=255)
                                                                    gTATCAAGATCGACGGTGACAATGCAATCGTAAATAACGAAGGAGATATCACCTTCTCAAAGGGTTGCTCCGTCACCCAGATTAATGGTGA‑TGACGCCACTGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                             <  2:1561018/150‑1 (MQ=17)
                                                                       tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                          >  2:1080651/1‑150 (MQ=255)
                                                                       tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                          >  1:1887426/1‑150 (MQ=255)
                                                                         cAGATCGACGGTGACAAGGCCACCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                        >  1:659081/1‑150 (MQ=255)
                                                                          aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                                       <  2:18750/150‑1 (MQ=255)
                                                                           gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                      >  2:1993180/1‑150 (MQ=255)
                                                                            aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                     >  2:1633822/1‑150 (MQ=255)
                                                                             tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAttg                                                                    <  1:1598684/150‑3 (MQ=255)
                                                                              cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGTTTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCaa                                                                   >  2:1989893/1‑150 (MQ=255)
                                                                               gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGA‑TGACGCCACg                                                                                                                               <  1:775052/90‑1 (MQ=14)
                                                                               gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGA‑TGACGCCACg                                                                                                                               >  2:775051/1‑90 (MQ=17)
                                                                                aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                                                 <  2:906772/150‑1 (MQ=255)
                                                                                aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                                                 <  1:1146518/150‑1 (MQ=255)
                                                                                aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACAAACGGTGGCACCGGAACACAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGATGACGGCAAGGATTCCACCTGTACGGAAATaaaag                                                                 >  1:1112339/1‑148 (MQ=17)
                                                                                 cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACgg                                                                >  1:853788/1‑150 (MQ=255)
                                                                                  ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGt                                                               <  2:1453076/150‑1 (MQ=255)
                                                                                    tGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaa                                                             >  1:21751/1‑150 (MQ=255)
                                                                                      aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCATATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACTGCAAAACCACCGTTGACGGCAAGGATTCCACCGTTACGGAAATCAACGGtaata                                                           <  1:1250598/150‑1 (MQ=255)
                                                                                        aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAAc                                                         <  2:1809167/150‑1 (MQ=255)
                                                                                            ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGA‑TGACGCCACGGCAAAC‑AACAACgg                                                                                                                <  1:2004497/91‑1 (MQ=255)
                                                                                            ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGA‑TGACGCCACGGCAAAC‑AACAACgg                                                                                                                >  2:2004496/1‑91 (MQ=255)
                                                                                              aTCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                                                                                                    <  2:767148/55‑1 (MQ=14)
                                                                                              aTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGa                                                                                                                                                    >  1:767149/1‑55 (MQ=255)
                                                                                                cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGA‑TGAc                                                                                                                                     >  2:731100/1‑67 (MQ=255)
                                                                                                     aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                            >  2:1331559/1‑150 (MQ=255)
                                                                                                       aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg                                          <  1:1382089/150‑1 (MQ=255)
                                                                                                        aCGAAGGCGAGAGCACCATCACAAACGGTGGCACCGGCACCCAGATTAAAGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACAACCGTTGACGGCAAGGAATCCACCGGTACGGAAATAAAAGGAAATAACGGAAAAGTTATTCAgg                                         >  1:55651/1‑150 (MQ=14)
                                                                                                             ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc                                    <  1:695707/150‑1 (MQ=255)
                                                                                                               cGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGCCGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATTACGGAAAAGTTATTCAGGACGGCGa                                  >  2:1074014/1‑150 (MQ=255)
                                                                                                                 agagCACCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACg                                                                                               >  2:1237733/1‑87 (MQ=18)
                                                                                                                 agagCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCTCCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATc                                <  2:898919/150‑1 (MQ=255)
                                                                                                                     cACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa                            <  2:1243004/150‑1 (MQ=255)
                                                                                                                       ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGt                                                                                                    >  2:1311270/1‑76 (MQ=39)
                                                                                                                       ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGt                                                                                                    <  1:1311271/76‑1 (MQ=39)
                                                                                                                       ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAGCTGGATg                          >  1:1275500/1‑150 (MQ=255)
                                                                                                                          tCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCa                       >  2:496950/1‑150 (MQ=255)
                                                                                                                           cACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAg                      >  2:1571477/1‑150 (MQ=255)
                                                                                                                            aCCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAg                                                                                          >  1:1755804/1‑81 (MQ=38)
                                                                                                                                aCGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAACCCACCGTTGACGGCAAGGAATCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcg                 <  2:1976135/150‑1 (MQ=255)
                                                                                                                                     ggcaccggcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt            <  1:570116/150‑1 (MQ=255)
                                                                                                                                       caccggcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAAACGGAAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa          >  2:676697/1‑150 (MQ=255)
                                                                                                                                        accggcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                 <  1:1672819/78‑1 (MQ=39)
                                                                                                                                        accggcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                 >  2:1672818/1‑78 (MQ=39)
                                                                                                                                        accggcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc         <  1:772462/150‑1 (MQ=255)
                                                                                                                                         ccggcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg        <  1:1080652/150‑1 (MQ=255)
                                                                                                                                          cggcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGc                                                                                                               >  1:1121944/1‑46 (MQ=255)
                                                                                                                                          cggcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGc                                                                                                               <  2:1121943/46‑1 (MQ=255)
                                                                                                                                          cggcaccCAGATAAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGGACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg       <  1:101963/150‑1 (MQ=255)
                                                                                                                                            gcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAg                                                                                          >  2:807604/1‑65 (MQ=40)
                                                                                                                                            gcaccCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAg                                                                                          <  1:807605/65‑1 (MQ=40)
                                                                                                                                               cccAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                      <  2:795272/66‑1 (MQ=39)
                                                                                                                                               cccAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                      >  1:795273/1‑66 (MQ=39)
                                                                                                                                               cccAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGAGGCGGTCACGGTATTg  >  2:1374031/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGA‑TGACGCCACGGCAAAC‑AACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTG  >  NC_000913/1466192‑1466484

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