Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,336 C→G 9.0% 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,3360CG 9.0% 722.9 / 45.0 224pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (103/98);  new base G (9/12);  total (114/110)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

ATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466187‑1466485
                                                                                                                                                     |                                                                                                                                                      
aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                                                                                                        <  2:859483/150‑1 (MQ=255)
aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                                                                                                        <  2:1219214/150‑1 (MQ=255)
  tGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa                                                                                                                                                      <  1:4589735/150‑1 (MQ=255)
         aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCGTGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                               <  2:1054817/150‑1 (MQ=255)
         aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                               <  1:4370437/150‑1 (MQ=255)
         aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                               <  2:3643765/150‑1 (MQ=255)
            ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGt                                                                                                                                            <  2:2271033/150‑1 (MQ=255)
             gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGAGACCAGGCCATAGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATAAAAGGtg                                                                                                                                           >  1:3701122/1‑150 (MQ=14)
              tGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                                                                                                          <  2:4776793/150‑1 (MQ=255)
                 cAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgatga                                                                                                                                       <  1:3583369/150‑1 (MQ=255)
                  aGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                                                                                                      <  1:3850793/150‑1 (MQ=255)
                          ggTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                                                                                                              <  1:4595658/150‑1 (MQ=255)
                             ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaaa                                                                                                                           <  1:3209436/150‑1 (MQ=255)
                              gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                          >  2:3810650/1‑150 (MQ=255)
                               aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                                                                                                         >  1:2612574/1‑150 (MQ=255)
                                tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaa                                                                                                                        <  2:1094508/150‑1 (MQ=255)
                                  aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                                                                                                      <  1:1123628/150‑1 (MQ=255)
                                  aCAAGGGCAACATGACCTTCACAGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGAACCATAACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                                                                                                      <  2:664302/150‑1 (MQ=18)
                                   cAAGGGCACCATGACCGTCACCGATCCGGTGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                     <  1:1793900/150‑1 (MQ=255)
                                   cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                     <  1:2888536/150‑1 (MQ=255)
                                   cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                     <  1:1916045/150‑1 (MQ=255)
                                    aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                    <  1:2828144/150‑1 (MQ=255)
                                     aGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                                                                                                   <  2:5114088/150‑1 (MQ=255)
                                       ggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                                                                                                                 >  1:2593028/1‑150 (MQ=255)
                                         cACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa                                                                                                               <  2:3294451/150‑1 (MQ=255)
                                          aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                       <  2:2429014/141‑1 (MQ=255)
                                          aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                       >  1:2429014/1‑141 (MQ=255)
                                          aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaa                                                                                                              >  1:5562104/1‑150 (MQ=255)
                                            cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                            >  1:5030757/1‑150 (MQ=255)
                                            cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCAGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                            >  2:2651195/1‑150 (MQ=255)
                                                 ccGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑Cg                                                                                                      >  1:603019/1‑150 (MQ=255)
                                                    tCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                                                                                                             <  1:4207927/109‑1 (MQ=255)
                                                      aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgatg                                                                                                                                        <  2:323358/112‑1 (MQ=255)
                                                      aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGAc                                                                                                 >  2:1589924/1‑150 (MQ=255)
                                                      aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGAc                                                                                                 <  2:5583375/150‑1 (MQ=255)
                                                      aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGAc                                                                                                 >  1:3439927/1‑150 (MQ=255)
                                                        cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccac                                                                                                         >  1:3547774/1‑141 (MQ=255)
                                                        cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACgg                                                                                               >  1:3683110/1‑150 (MQ=255)
                                                        cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACgg                                                                                               >  1:4033026/1‑150 (MQ=255)
                                                        cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACgg                                                                                               >  1:4792462/1‑150 (MQ=255)
                                                        cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGCCAAACAACAACCGCAAAACCAC‑CGTTCACgg                                                                                               >  1:3352715/1‑150 (MQ=255)
                                                              ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGa                                                                                         >  2:1931773/1‑150 (MQ=255)
                                                               gAGTCGATGGGTATCCAGATCGACGGTGTCCAGGCCCTCGTCTCTAACGTAGGCGAGAGGACCCTCACGAACGGTGGCCCCGGCACCGAGATTACTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGAt                                                                                        <  2:3154401/150‑1 (MQ=17)
                                                               gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGAt                                                                                        <  2:2037086/150‑1 (MQ=255)
                                                                 gTCCATGGGTATCCAGAGCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTc                                                                                      <  1:1062810/150‑1 (MQ=255)
                                                                   ccATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCa                                                                                    >  1:1561154/1‑150 (MQ=255)
                                                                     aTGGGTATCCAGATCGACGGTGACAAGGCCATTGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCAcc                                                                                  >  2:1977623/1‑150 (MQ=255)
                                                                     aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCAcc                                                                                  >  1:1084619/1‑150 (MQ=255)
                                                                     aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCACACAACCACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCAcc                                                                                  >  2:4184712/1‑150 (MQ=255)
                                                                     aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCAcc                                                                                  <  2:3485921/150‑1 (MQ=255)
                                                                      tGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGGGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCg                                                                                 <  2:580891/150‑1 (MQ=255)
                                                                      tGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCg                                                                                 <  1:5250042/150‑1 (MQ=255)
                                                                       gggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                                                                                                    <  2:2609610/83‑1 (MQ=255)
                                                                         gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTa                                                                              >  2:1314650/1‑150 (MQ=255)
                                                                          tATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTAc                                                                             <  1:5042926/150‑1 (MQ=255)
                                                                           aTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACg                                                                            <  1:2476961/150‑1 (MQ=255)
                                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  2:4720783/1‑150 (MQ=255)
                                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  1:1439955/1‑150 (MQ=255)
                                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  2:618515/1‑150 (MQ=255)
                                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  1:542483/1‑150 (MQ=255)
                                                                              cAGATCGACGGTGACAAGGACATCGTCAATTACGAAGGCGAGAGCACCATCACCAACGGTGGCTCCGGCACCCTGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                         <  1:559327/150‑1 (MQ=255)
                                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTaac                                                                                                                                               <  2:1376400/80‑2 (MQ=255)
                                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                                        >  1:4820275/1‑150 (MQ=255)
                                                                                 aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGTCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                      <  2:2900220/150‑1 (MQ=255)
                                                                                 aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                      >  1:4292735/1‑150 (MQ=255)
                                                                                    gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAAttgct                                                                   <  2:4033026/150‑5 (MQ=255)
                                                                                    gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAAttgct                                                                   <  2:3683110/150‑5 (MQ=255)
                                                                                    gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAAttgct                                                                   <  2:1813289/150‑5 (MQ=255)
                                                                                      cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAAttgctgg                                                                 <  2:5626527/150‑7 (MQ=11)
                                                                                           aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata                                                            >  2:5324546/1‑150 (MQ=255)
                                                                                              aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACg                                                         <  1:1466943/150‑1 (MQ=255)
                                                                                                   aTCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTGGCGCCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                    <  2:4811108/150‑1 (MQ=255)
                                                                                                   aTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACTCCACGGCAAAAAACAACTGCAAAAACAC‑CTTTGACTGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                    >  1:5297369/1‑150 (MQ=17)
                                                                                                    tCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                       <  2:4661718/130‑1 (MQ=255)
                                                                                                    tCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                       >  1:4661718/1‑130 (MQ=255)
                                                                                                    tCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                       >  1:216890/1‑130 (MQ=255)
                                                                                                      gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACTGCGCAGCACACAACAACTG‑AACTCCACACGCTCACGGAGAGCAGTCAAGCGGTCCGGAGATCTACGGTAATAACggtcgagtg                                                 >  1:1495738/1‑143 (MQ=9)
                                                                                                      gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAAAAACAACTGAAAAACCAA‑CGATCACGGAAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtt                                                 >  1:2280558/1‑150 (MQ=17)
                                                                                                       tCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACgg                                                                           >  1:5805309/1‑123 (MQ=255)
                                                                                                           tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCa                                            >  1:90540/1‑150 (MQ=255)
                                                                                                            aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg                                           <  1:938601/150‑1 (MQ=255)
                                                                                                            aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg                                           >  1:4721855/1‑150 (MQ=255)
                                                                                                            aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg                                           <  1:1353811/150‑1 (MQ=255)
                                                                                                              cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGa                                         <  2:2593028/150‑1 (MQ=255)
                                                                                                                     gagaGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt                                  <  2:3131620/150‑1 (MQ=255)
                                                                                                                        agCACCATCACCAACGGTGGCACCGGGACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCTCCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg                               <  1:5421159/150‑1 (MQ=255)
                                                                                                                            ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATAAACGGAAATAACGGGAATGTTATTCAGGACGGCGATCTGGATg                           >  1:5673540/1‑150 (MQ=255)
                                                                                                                              aTCACCATCGGTGGCAGCGGCATCGAGATTACTGGTGATGACGCCACTGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTc                         <  2:5030572/150‑1 (MQ=17)
                                                                                                                                 aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                                        >  1:2347867/1‑100 (MQ=37)
                                                                                                                                      cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcgg                 <  2:5417458/150‑1 (MQ=255)
                                                                                                                                         tGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGcggcggcgg              >  2:4534602/1‑150 (MQ=255)
                                                                                                                                         tGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGcggcggcgg              >  2:1420320/1‑150 (MQ=255)
                                                                                                                                           gcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCa                                                                                             >  1:4027581/1‑69 (MQ=40)
                                                                                                                                           gcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCa                                                                                             <  2:1559653/69‑1 (MQ=40)
                                                                                                                                                  caccGAGATTACGGGTGATGACGCCACGGAAAACAACAAAGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt     <  2:434000/150‑1 (MQ=14)
                                                                                                                                                  caccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt     <  1:3810650/150‑1 (MQ=255)
                                                                                                                                                   accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCaa                                                                                            <  1:3057066/62‑1 (MQ=40)
                                                                                                                                                   accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt    <  2:4792462/150‑1 (MQ=255)
                                                                                                                                                     ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGAGGTCAGCGGCGGCGGTCACGGTATTGa  >  2:2167420/1‑150 (MQ=255)
                                                                                                                                                     ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAAAAAAGGAAAAGGTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATAGa  >  2:5592365/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAC‑CGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466187‑1466485

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: ‑