Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 561,943 G→A 5.3% T749T (ACG→ACA sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,9430GA 5.3% 301.5 / 3.0 93T749T (ACG→ACAsfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (49/38);  new base A (3/2);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGA  >  NC_000913/561801‑562087
                                                                                                                                              |                                                                                                                                                
ccAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATggg                                                                                                                                             <  6:148338/148‑1 (MQ=255)
 cAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATgggg                                                                                                                                            >  8:57110/1‑148 (MQ=255)
 cAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATgggg                                                                                                                                            >  8:146204/1‑148 (MQ=255)
 cAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACAATgggg                                                                                                                                            >  3:165397/1‑148 (MQ=255)
 aaaTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACAATgggg                                                                                                                                            >  1:290618/2‑148 (MQ=255)
  aaTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGaa                                                                                                                                          >  3:184026/1‑149 (MQ=255)
   aTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAt                                                                                                                                         <  3:158059/149‑1 (MQ=255)
   aTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAt                                                                                                                                         >  2:4279/1‑149 (MQ=255)
      ttttGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTc                                                                                                                                       >  7:243808/1‑148 (MQ=255)
      ttttGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTc                                                                                                                                       >  5:345926/1‑148 (MQ=255)
       tttGATTAAAGCGCCTGGCGCAGGCGGTGAACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTCTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATAAGATGGGGAATTCCa                                                                                                                                     >  5:98289/1‑149 (MQ=255)
         tGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCAt                                                                                                                                    <  8:122946/148‑1 (MQ=255)
          gATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATc                                                                                                                                   <  7:132418/148‑1 (MQ=255)
          gATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATc                                                                                                                                   >  4:240459/1‑148 (MQ=255)
          gATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATc                                                                                                                                   <  3:92377/148‑1 (MQ=255)
          gATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATc                                                                                                                                   <  1:137811/148‑1 (MQ=255)
           aTTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAGATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGa                                                                                                                                 >  3:272445/1‑149 (MQ=255)
           aTTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACTATGGGGATTTCCATCGa                                                                                                                                 >  6:333132/1‑149 (MQ=255)
            ttAAAGGGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGCCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGAt                                                                                                                                <  5:89322/149‑1 (MQ=255)
            ttAAAGCGCCTGGCGCAGGCGGTGTTCGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGAt                                                                                                                                <  7:31751/149‑1 (MQ=255)
               aaGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGtt                                                                                                                             <  8:129372/149‑1 (MQ=255)
               aaGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGtt                                                                                                                             >  3:268485/1‑149 (MQ=255)
               aaGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGtt                                                                                                                             <  6:267439/149‑1 (MQ=255)
                   gcCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaa                                                                                                                         <  3:240459/149‑1 (MQ=255)
                         cgcAGGCTGTGTACTCATTGAAAATCAAACTTGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAata                                                                                                                   <  6:331476/149‑1 (MQ=255)
                            aGGCGGTGTACGCATTGAAGATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCTCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTa                                                                                                                <  4:272445/149‑1 (MQ=255)
                            aGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTa                                                                                                                <  2:320134/149‑1 (MQ=255)
                              gCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTag                                                                                                               >  2:333321/1‑148 (MQ=255)
                                 gtgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTagcagc                                                                                                           >  6:324401/1‑149 (MQ=255)
                                 gtgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACAATGGGGAATTCCATCGATGTTGAAAAAAATATTagcagc                                                                                                           <  4:165397/149‑1 (MQ=255)
                                         aTTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGcc                                                                                                   >  8:51894/1‑149 (MQ=255)
                                         aTTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATCCCATCGATGTTGAAAAAAATATTagcag                                                                                                            >  2:45750/1‑140 (MQ=255)
                                         aTTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATCCCATCGATGTTGAAAAAAATATTagcag                                                                                                            <  1:45750/140‑1 (MQ=255)
                                          ttGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTcc                                                                                                                                      >  7:160234/1‑113 (MQ=255)
                                          ttGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTcc                                                                                                                                      <  8:160234/113‑1 (MQ=255)
                                          ttGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTATCAGCGTTGTGcc                                                                                                   >  1:61036/1‑148 (MQ=255)
                                           tGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGaa                                                                                                                                          >  3:288799/1‑108 (MQ=255)
                                           tGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGaa                                                                                                                                          <  4:288799/108‑1 (MQ=255)
                                               aaTCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATggg                                                                                                                                             >  8:179776/1‑101 (MQ=255)
                                               aaTCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTa                                                                                                                <  3:36804/130‑1 (MQ=255)
                                               aaTCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTa                                                                                                                >  4:36804/1‑130 (MQ=255)
                                                aTCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCaa                                                                                            <  1:333321/149‑1 (MQ=255)
                                                       tGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGt                                                                                     >  7:32930/1‑149 (MQ=255)
                                                       tGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGt                                                                                     >  2:179500/1‑149 (MQ=255)
                                                         gCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTg                                                                                   >  7:89082/1‑149 (MQ=255)
                                                           aTTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGt                                                                                 <  1:4279/149‑1 (MQ=255)
                                                            ttttAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGtt                                                                                >  1:140758/1‑149 (MQ=255)
                                                                  ccGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGc                                                                           >  7:330318/1‑148 (MQ=255)
                                                                  ccGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGAAGGGGAATTCCATCGATGTAGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGc                                                                           >  4:238770/1‑148 (MQ=255)
                                                                      cTGGCGCGGCTATGCTGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACAATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAAtt                                                                      <  2:290618/149‑1 (MQ=255)
                                                                       tGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAAttt                                                                     >  4:187890/1‑149 (MQ=255)
                                                                           gcgGCTATGCGGTGATGCTGTATGCCACGGTTTATCTGTATAACAGTATCTCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAata                                                                                                                   <  3:271019/99‑1 (MQ=255)
                                                                             ggCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATAc                                                               >  2:292789/1‑149 (MQ=255)
                                                                               cTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATAccc                                                             <  7:57110/149‑1 (MQ=255)
                                                                                  tGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCa                                                          <  4:184026/149‑1 (MQ=255)
                                                                                       tGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGc                                                     <  5:336727/149‑1 (MQ=255)
                                                                                       tGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGAAGCCTTGTGCCGCCGCAAGGCGCGTTGGTTCGTGCCAAGTTTGATACCCGCATat                                                       >  6:19176/1‑146 (MQ=255)
                                                                                            cTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCg                                                <  6:137337/149‑1 (MQ=255)
                                                                                               tATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGgc                                             >  1:49148/1‑149 (MQ=255)
                                                                                                  gCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGTGGGCGCt                                          >  2:182544/1‑149 (MQ=255)
                                                                                                   ccACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGAGGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTAGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCt                                          >  6:309800/1‑148 (MQ=255)
                                                                                                   ccACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACAATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCt                                          >  8:78625/1‑148 (MQ=255)
                                                                                                    cACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCa                                        <  7:146204/149‑1 (MQ=255)
                                                                                                        gTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGATTTCCATCGATGTTGAAAAAAATATTAGCAGGGTTGTGCCGACGCAAGGCGCGTTGGTTCGGGCGAATTTTGATACCCGCATAGGCGTGCGGGCGCGCCTTAc                                    >  6:326133/1‑149 (MQ=255)
                                                                                                             tCGGTATAACCGTATCTCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCTAGGCGCGTTGGTTCGTGGCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTa                               <  5:115206/149‑1 (MQ=255)
                                                                                                                gTATAACCGTATAGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAACGCGCGTTCGTTCGTGCCAATTTTGATACCCGCATACGCGTGCTGGCGCTCATTACCGTTAcc                             >  2:288174/1‑148 (MQ=255)
                                                                                                                 tataACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGgcgc                                           >  6:199835/1‑133 (MQ=255)
                                                                                                                 tataACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGgcgc                                           <  5:199835/133‑1 (MQ=255)
                                                                                                                   taACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAgg                         >  1:180274/1‑149 (MQ=255)
                                                                                                                     aCCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGTGGGCGCTCATTACCGTTACCCAGGGc                       <  1:182544/149‑1 (MQ=255)
                                                                                                                     aCCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGc                       <  1:179500/149‑1 (MQ=255)
                                                                                                                       cGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTa                               <  6:160865/139‑1 (MQ=255)
                                                                                                                       cGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTa                               >  5:160865/1‑139 (MQ=255)
                                                                                                                        gTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGCCGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCgg                     >  1:64089/1‑148 (MQ=255)
                                                                                                                            cgcgCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAAc                >  2:95451/1‑149 (MQ=255)
                                                                                                                             gcgcTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaaaaa                                                                                                                     <  1:146457/47‑1 (MQ=255)
                                                                                                                             gcgcTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaaaaa                                                                                                                     >  2:146457/1‑47 (MQ=255)
                                                                                                                             gcgcTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAAcc               >  1:30891/1‑149 (MQ=255)
                                                                                                                             gcgcTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACAGTTACCCAGGGCGGAAAAcc               <  4:268485/149‑1 (MQ=255)
                                                                                                                               gcTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCgg             >  4:88599/1‑149 (MQ=255)
                                                                                                                               gcTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCgg             <  2:140758/149‑1 (MQ=255)
                                                                                                                                   gATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCCGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCAAAGGCGTGAGGGCGCTCATTACCGTTACCAAGGGAGGAAAAACGCTGcc         >  4:222354/1‑149 (MQ=255)
                                                                                                                                   gATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGcc         >  6:194057/1‑149 (MQ=255)
                                                                                                                                   gATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGGTGGTTCGTGCCAATTTTGATACCCGCCTAGGCGGGCGTGCGCTCAGGACCGTTACCCCGGGCGGAAAACCGGTGcc         >  2:276343/1‑149 (MQ=255)
                                                                                                                                    aTACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCg        >  2:176470/1‑149 (MQ=255)
                                                                                                                                        cAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTc                                         <  3:241180/112‑1 (MQ=255)
                                                                                                                                        cAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTc                                         >  4:241180/1‑112 (MQ=255)
                                                                                                                                        cAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTg    <  6:71524/149‑1 (MQ=255)
                                                                                                                                        cAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTg    >  5:71524/1‑149 (MQ=255)
                                                                                                                                        cAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTg    >  6:55557/1‑149 (MQ=255)
                                                                                                                                          aTACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGGGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGa  <  5:324401/149‑1 (MQ=255)
                                                                                                                                          aTACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATa                                                        <  6:230480/95‑1 (MQ=255)
                                                                                                                                          aTACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATa                                                        >  5:230480/1‑95 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGA  >  NC_000913/561801‑562087

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