Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,345 T→C 7.1% pseudogene (954/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,3450TC 7.1% 428.5 / 15.3 151pseudogene (954/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (50/80);  new base C (4/6);  total (62/89)
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

Only 100 of 156 total aligned reads displayed.

ACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466196‑1466494
                                                                                                                                                     |                                                                                                                                                     
aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                                       <  2:2274992/150‑1 (MQ=255)
aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                                       <  2:2271836/150‑1 (MQ=255)
aCCGGTGACAGCGCGACGGTGGATAACAAGCGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCAACACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                                       <  1:3521917/150‑1 (MQ=255)
 ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                                      <  2:3951164/150‑1 (MQ=255)
 ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                                      <  1:1384250/150‑1 (MQ=255)
 ccGGTGACAGCGCGACGGTGGAGAACAAGGGCACCATGACCGTCACCGATCCGGCGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAAGAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATGAATg                                                                                                                                                      <  2:3951172/150‑1 (MQ=255)
  cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                                                                                                                     >  2:3673821/1‑150 (MQ=255)
  cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                                                                                                                     >  2:1296855/1‑150 (MQ=255)
    gTGACAGCGCGACGGTGTATAACTAGGGCACCATGACCGTCACCGATCCGTAGTCCATGGTTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCTACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                                   <  1:1523972/150‑1 (MQ=18)
    gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                                   >  1:1457703/1‑150 (MQ=255)
     tGACAGCGCGACTGTGGATAACAAGGGAACCTTGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGTCGAGAGCTCCATCAGCAACGGTGGCACCGGCACCCAGTTTAATGGtga                                                                                                                                                  <  2:1734882/150‑1 (MQ=18)
     tGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                                                                                                                  <  2:1734841/150‑1 (MQ=255)
     tGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                                                                                                                  <  2:3131242/150‑1 (MQ=255)
      gACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                                                                                                                 <  1:2271770/150‑1 (MQ=255)
      gACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                                                                                                                 <  1:2063410/150‑1 (MQ=255)
         aGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                                                                                                              <  2:1655437/150‑1 (MQ=255)
             cgACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCa                                                                                                                                          >  1:812287/1‑150 (MQ=255)
             cgACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCa                                                                                                                                          <  2:2513487/150‑1 (MQ=255)
              gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCAc                                                                                                                                         <  2:1139321/150‑1 (MQ=255)
               aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACg                                                                                                                                        >  1:1888130/1‑150 (MQ=255)
                 ggTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                                                                                                                      <  2:2383091/150‑1 (MQ=255)
                    ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaaa                                                                                                                                   <  1:3914619/150‑1 (MQ=255)
                    ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaaa                                                                                                                                   <  1:2115289/150‑1 (MQ=255)
                    ggAAAACAAGGGCACCATGACCGTCACAGATCCGGAGTCCATGGTTATACAGAACGACTGAGACAAGGACATCGTCAAAAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCAACCAGAGAAAGGGTGATGACGCCACGGCaaa                                                                                                                                   >  1:976520/1‑150 (MQ=17)
                           aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGAc                                                                                                                                              >  1:2073079/1‑132 (MQ=255)
                           aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                            <  1:1315553/150‑1 (MQ=255)
                             gggCACCATGACCGTCACCGATCCGGTGTCCTTGGGTATCCAGATCGACGGTGACAAGGCCATCTTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACCCAGATTAATGGTGTTGACGCCACGGCAAACAACAACgg                                                                                                                          <  1:3719985/150‑1 (MQ=18)
                                 aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaa                                                                                                                      <  2:1457703/150‑1 (MQ=255)
                                 aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaa                                                                                                                      <  1:1652875/150‑1 (MQ=255)
                                   cATGACCGTCACCGATCGGGAGTCCTTGGGTATCCAGATCGACGGTGACAAGGCCTTCGTCAATAACGAAGGCGAGAGCACCATCACCAAGGGTGGCACCGGCACCCAGTTTATTGGTGATGACGCCAGGGCAAACAACAACGGCAAAac                                                                                                                    <  1:1222774/150‑1 (MQ=18)
                                   cATGACCGTCACCGATCCGGAGTCCATGGGTATCCATATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                                    <  1:2533040/150‑1 (MQ=255)
                                   cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                                    >  1:1000949/1‑150 (MQ=255)
                                   cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                                    >  1:2416314/1‑150 (MQ=255)
                                   cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                                    <  1:2198318/150‑1 (MQ=255)
                                        ccGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                                                                                                               >  2:1089069/1‑150 (MQ=255)
                                         cGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGt                                                                                                              >  1:2947978/1‑150 (MQ=255)
                                             aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                                                          <  2:1897806/150‑1 (MQ=255)
                                             aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAGTAATGGTGATGACGCCACGGCCAACAACAACGGCAAAACCACCGTTGAc                                                                                                          >  1:2615440/1‑150 (MQ=255)
                                               cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                                        >  2:2519750/1‑150 (MQ=255)
                                                gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                                                                       >  1:1690872/1‑150 (MQ=255)
                                                 aTCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCGCCGGTGACGGCa                                                                                                      >  2:1899735/1‑150 (MQ=255)
                                                  tCCGGAGTCCATGGTTATCAAGATCGACGGTGACAAGCCCATCGTCAATAACGAAGGCGAGAGCACCAACACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  1:4077719/150‑1 (MQ=255)
                                                  tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     >  2:2572770/1‑150 (MQ=255)
                                                  tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  2:118550/150‑1 (MQ=255)
                                                  tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  2:118646/150‑1 (MQ=255)
                                                  tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  1:3169300/150‑1 (MQ=255)
                                                  tCCGGAGCCCATGGGTATACAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACCCTGATGTAGGGTGATGACGCGACGGCAAACAACAACGGCAACCCCACCGTTGACGGCaa                                                                                                     <  1:1738068/150‑1 (MQ=25)
                                                     ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                                                                                                  <  1:3115030/150‑1 (MQ=255)
                                                      gAGTCCATGGGTATCCAGTTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGTGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                                 <  2:2443609/150‑1 (MQ=255)
                                                      gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGc                                                                                                                                            <  2:1357003/107‑1 (MQ=255)
                                                      gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                                 <  1:717880/150‑1 (MQ=255)
                                                       aGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACACAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                                <  2:4229708/150‑1 (MQ=255)
                                                        gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTc                                                                                               <  2:2882180/150‑1 (MQ=255)
                                                         tCCATGGGCACCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCTCCTTCGGTGGCTCCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                                                              <  1:2080180/150‑1 (MQ=18)
                                                           cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                                                          <  2:2052984/136‑1 (MQ=255)
                                                           cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                                                          >  1:2052984/1‑136 (MQ=255)
                                                            aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAAAGCTGATGAACACACGGCAAACAACAACGGCAAAACAACCGTTGACGGCAAGGATTCCAcc                                                                                           >  1:156023/1‑150 (MQ=17)
                                                            aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                           >  1:833896/1‑150 (MQ=255)
                                                                ggaTCCAGATCGACGGTGACAAGGCCATAGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCTCCGGCACCCAGATTAATGGTGATGACGCCACGGAAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                       <  2:895405/148‑1 (MQ=255)
                                                                gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                       <  2:3561776/150‑1 (MQ=255)
                                                                  aTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTaag                                                                                     >  1:1271901/1‑148 (MQ=255)
                                                                  aTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACg                                                                                     <  1:3621553/150‑1 (MQ=255)
                                                                   tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTc                                                                                               <  2:207862/139‑1 (MQ=255)
                                                                   tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                                    <  1:456866/150‑1 (MQ=255)
                                                                   tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                                    >  1:1390982/1‑150 (MQ=255)
                                                                   tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGAAAACAACAACGGCAAAACCACCGAAGACGGCAAGGATTCCACCGGTACgg                                                                                    >  1:2205801/1‑150 (MQ=255)
                                                                     cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                                                                                                  <  1:505349/134‑1 (MQ=255)
                                                                     cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                                                                                                  >  2:505349/1‑134 (MQ=255)
                                                                     cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                                  >  1:679709/1‑150 (MQ=255)
                                                                        aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                               >  2:715378/1‑150 (MQ=255)
                                                                        aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                               >  1:4113840/1‑150 (MQ=255)
                                                                        aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACAAAGATAAATCAACATGAAGCAAACCCGAAAAACAACCGCAAAACCAACGTTGACTGCAAGGATTCCACCGGTACGGAAGTc                                                                               >  1:1184223/1‑150 (MQ=17)
                                                                          cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCaa                                                                             <  2:3485657/150‑1 (MQ=255)
                                                                           gACGGTGACCAGGCCATCGTCAATACCGAAGGCGAGAGCACCCTCACCAACGGTGGCACCGGCACCCAGAGTACTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAtggct                                                                            <  2:2615440/150‑5 (MQ=17)
                                                                           gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                                                            <  1:3279178/150‑1 (MQ=255)
                                                                                    taGGGCATCGTCTATTACGAAGGCGAGAGCACCATCACCAACGGTGGCTCCGGCACCCAGATTACGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAAc                                                                   <  1:907388/149‑1 (MQ=17)
                                                                                    atggCCATCGTCTATAACGCAGGCGAGTGCACCAGCACCACCGGTGGCGACGTCACCGTGATTAAGGGTGATTACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAAc                                                                   <  2:3635901/148‑1 (MQ=17)
                                                                                     aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACg                                                                  <  1:1316767/150‑1 (MQ=255)
                                                                                     aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACg                                                                  >  1:4234110/1‑150 (MQ=255)
                                                                                        ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaa                                                               <  1:1089069/150‑1 (MQ=255)
                                                                                         cATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                                                             <  2:1372359/103‑1 (MQ=255)
                                                                                         cATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                                                             >  1:1372359/1‑103 (MQ=255)
                                                                                         cATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                                                   >  2:953808/1‑113 (MQ=255)
                                                                                         cATCGTAAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                                                   <  1:953808/113‑1 (MQ=255)
                                                                                           tCGTCAATCACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCATATTATTGGTGATGACGCCACGGCTATCAACATCGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAg                                                            <  1:3890104/150‑1 (MQ=18)
                                                                                              tCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTa                                                         >  1:3330607/1‑150 (MQ=255)
                                                                                                          gCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACCACGGCTAAACCCCCGAAGACGGCTCGGATTCCCCCGGTACGGAAATCAACGGTAATCACGGAACTGTTATTCGGGTCGGCg                                             >  1:3964236/1‑150 (MQ=17)
                                                                                                          gCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACAACCGATGACGGCTAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAGGTTGTTCAGGACGGCg                                             >  1:1853100/1‑150 (MQ=17)
                                                                                                                    cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGAGGt                                   >  2:3852257/1‑150 (MQ=255)
                                                                                                                        aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgc                               <  1:3617101/150‑1 (MQ=255)
                                                                                                                         ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcg                              >  2:3204720/1‑150 (MQ=255)
                                                                                                                             cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGTAAATGTTATTCAGGACGGCGATCTGGATGTCAgcggcgg                          <  1:1772073/150‑1 (MQ=255)
                                                                                                                              ggTGGCACCGGCACCCATATTAATGGTTATGACGCCACGGAAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGTACGGCGATCTGGATGTCAgcggcggc                         <  1:175315/150‑1 (MQ=255)
                                                                                                                                    accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc                   >  1:342791/1‑150 (MQ=255)
                                                                                                                                     ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                  <  2:2072319/150‑1 (MQ=255)
                                                                                                                                          accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt             >  2:3082902/1‑150 (MQ=255)
                                                                                                                                          accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt             >  2:2898722/1‑150 (MQ=255)
                                                                                                                                               gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGatat        <  1:3673821/150‑1 (MQ=255)
                                                                                                                                               gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGatat        <  1:1296586/150‑1 (MQ=255)
                                                                                                                                                     tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg  >  1:1494675/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466196‑1466494

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