Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,318,079 A→T 5.4% intergenic (‑73/+134) rimP ← / ← metY ribosome maturation factor for 30S subunits/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,318,0790AT 5.4% 371.9 / 8.8 130intergenic (‑73/+134)rimP/metYribosome maturation factor for 30S subunits/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (67/56);  new base T (4/3);  total (71/59)
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.70e-01

Only 100 of 131 total aligned reads displayed.

TGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCC  >  NC_000913/3317931‑3318226
                                                                                                                                                    |                                                                                                                                                   
tGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACCAGCCCACCCCCAACACATAAAAAAAGGGCCTAAAGCGCAGTTATTCTGTAGTCAGATAACAAAAAACACCGATaa                                                                                                                                                    >  1:324894/1‑150 (MQ=255)
     aCCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCgg                                                                                                                                               >  2:1666617/1‑150 (MQ=255)
     aCCAGTTCAAAACCCAGGGCCTAAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAAAGTGTACAAGCCAACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCgg                                                                                                                                               >  1:1053395/1‑150 (MQ=255)
      ccAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTATGTAGTCAGATAACAAAAAACCCCGATAAATCggg                                                                                                                                              >  2:1536114/1‑150 (MQ=255)
        aGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACACGCCCCCCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCggggg                                                                                                                                            >  2:898320/1‑149 (MQ=255)
            cAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCCCCCCCAAGCCATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAAACCCGATAAATCGGGGCTTta                                                                                                                                        >  2:1282290/1‑150 (MQ=255)
            cAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATc                                                                                                                                                 <  2:113966/141‑1 (MQ=255)
            cAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCAAGTTATTCTGTAGTCAAATAACAAAAAACCCCGATAAAACGGGGCTTta                                                                                                                                        >  2:1709486/1‑150 (MQ=255)
             aaaaCCAAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGGGtttttt                                                                                                                                       <  2:1007320/150‑3 (MQ=255)
              aaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtata                                                                                                                                      >  1:1308293/1‑150 (MQ=255)
               aaCCCATGGCCTCAACTGTCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtatat                                                                                                                                     <  1:1331737/150‑1 (MQ=255)
                aCCCAGTGCCTCAACTGGGGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtatata                                                                                                                                    <  2:1456446/150‑1 (MQ=255)
                aCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtatata                                                                                                                                    <  1:1430621/150‑1 (MQ=255)
                aCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtatata                                                                                                                                    >  2:356398/1‑150 (MQ=255)
                aCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtatata                                                                                                                                    >  1:1078247/1‑150 (MQ=255)
                 cccAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATaa                                                                                                                                   >  1:159869/1‑150 (MQ=255)
                         cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAccc                                                                                                                           >  1:1707603/1‑150 (MQ=255)
                            aaCTTGCCCATTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCtat                                                                                                                        <  1:203241/150‑1 (MQ=255)
                            aaCTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCtat                                                                                                                        >  2:710166/1‑150 (MQ=255)
                            aaCTGGCGCAGTAATCATCTCTGATAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAAATGACCCCtat                                                                                                                        >  1:1245242/1‑150 (MQ=255)
                             aCTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCtata                                                                                                                       <  1:802902/150‑1 (MQ=255)
                             aCTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCtata                                                                                                                       >  1:1133310/1‑150 (MQ=255)
                              cTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATaa                                                                                                                      >  1:1514824/1‑150 (MQ=255)
                              cTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATaa                                                                                                                      >  1:1161340/1‑150 (MQ=255)
                                 gcgcAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGACTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAATTCGGGGCTTTATATAACTGAACCCTATAACCg                                                                                                                   <  1:903628/150‑1 (MQ=255)
                                   gcAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAACACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCa                                                                                                                 <  1:313943/150‑1 (MQ=255)
                                     aGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCaag                                                                                                               <  1:867108/150‑2 (MQ=255)
                                         atcatcTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCACCTGCg                                                                                                           >  2:664360/1‑150 (MQ=255)
                                         atcatcTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCg                                                                                                           <  1:1560561/150‑1 (MQ=255)
                                                ctGTTAATTTTTGCTCTAATGTGGACAAGCCCCCCCCCAAGACATAAAAAACGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCTATAAATCGGGGCGTTATATAACTGAACCCTATAACCGCAACTGCGGTCTggc                                                                                                    >  2:73576/1‑149 (MQ=255)
                                                 tGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAg                                                                                                   <  2:434316/150‑1 (MQ=255)
                                                 tGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAg                                                                                                   >  2:358511/1‑150 (MQ=255)
                                                     aaTTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTTTTGTGTAGTCAGATAACAAAAAACCAAGATAAATCGGTGCTTTATATAAATTAACAATAATACAGCAACTGCCGTCTGTAGCaat                                                                                               >  1:1285992/1‑148 (MQ=255)
                                                     aaTTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACt                                                                                               >  1:1499557/1‑150 (MQ=255)
                                                         tttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGCTAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTcc                                                                                           <  1:1341913/150‑1 (MQ=255)
                                                           tGCTCAAATGTCGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCATTATATAACTGAACCCTATAACCGCAACAGCGGTCTGGAGCACTTTCCAg                                                                                         <  1:193369/150‑1 (MQ=255)
                                                             ctctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                       >  1:1499370/1‑150 (MQ=255)
                                                              tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCTGTCTGGAGCACTTTCCAGAAg                                                                                      <  1:1536114/150‑1 (MQ=255)
                                                              tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAg                                                                                      <  2:643678/150‑1 (MQ=255)
                                                                 aaTGTGGACAAGCCCACCCACAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAt                                                                                   <  2:1596086/150‑1 (MQ=255)
                                                                  aTGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtt                                                                                  >  1:393621/1‑150 (MQ=255)
                                                                    gtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttt                                                                                <  1:574107/150‑1 (MQ=255)
                                                                      ggACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttttt                                                                              >  2:1732531/1‑150 (MQ=255)
                                                                      ggACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATAAAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGCAGGAtttttt                                                                              >  2:1059539/1‑150 (MQ=255)
                                                                       gACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGACCCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTc                                                                             >  2:635756/1‑150 (MQ=255)
                                                                       gACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCTGAAGGATTTTTTc                                                                             <  1:374801/150‑1 (MQ=255)
                                                                        aCAGGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCa                                                                            >  2:135795/1‑150 (MQ=255)
                                                                        aCAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCa                                                                            <  1:1080897/150‑1 (MQ=255)
                                                                        aCAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCa                                                                            <  2:1732794/150‑1 (MQ=255)
                                                                         cAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaa                                                                           <  1:1650791/150‑1 (MQ=255)
                                                                           aGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaa                                                                           >  1:1409370/1‑148 (MQ=255)
                                                                           aGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaa                                                                           <  2:1409370/148‑1 (MQ=255)
                                                                           aGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAAt                                                                         >  1:758373/1‑150 (MQ=255)
                                                                           aGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAAt                                                                         >  1:891140/1‑150 (MQ=255)
                                                                            gCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATc                                                                        >  1:755571/1‑150 (MQ=255)
                                                                             cccACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATcc                                                                       >  2:490629/1‑150 (MQ=255)
                                                                             cccACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATcc                                                                       >  2:66799/1‑150 (MQ=255)
                                                                             cccACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACAGCAACTGCGGTCTGGAGAACTTTCCAGAAGGATTTTTTCAAATcc                                                                       >  2:160467/1‑150 (MQ=255)
                                                                               cACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                     >  2:1677950/1‑150 (MQ=255)
                                                                                 cccccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGTATTTTTTCAAATCCCACt                                                                   >  1:1123014/1‑150 (MQ=255)
                                                                                    ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACg                                                                <  1:1534515/150‑1 (MQ=255)
                                                                                          cATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGcc                                                          >  2:1350189/1‑150 (MQ=255)
                                                                                           aTAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCg                                                         <  2:563854/150‑1 (MQ=255)
                                                                                            tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGa                                                        <  2:1724771/150‑1 (MQ=255)
                                                                                             aaaaaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAACCTATAACCGCAACTGCGGTCTTGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                       >  2:85066/1‑150 (MQ=255)
                                                                                                aaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCACACTACGAAGGCGGAAGTc                                                    >  2:154033/1‑150 (MQ=255)
                                                                                                 aaatgtccTAAAGCCCAGTTTTTCTGTAGTAAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCACACTACGAAGGCCGAAGTCt                                                   <  1:1226287/144‑1 (MQ=255)
                                                                                                    gggCCTAAAGCAGAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGTGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTca                                                <  2:387746/150‑1 (MQ=255)
                                                                                                     ggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAACCCCCGATAAATCGGGGCTTTATATAACGGAACCCTATAACCGCAACTGCGGTCTGGAGCAATTTCCAGAAGGATTATTTCAAATCCCACTACGAAGGCCTAAGTCTTcac                                               >  1:299566/1‑150 (MQ=255)
                                                                                                     ggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTcac                                               <  1:1195482/150‑1 (MQ=255)
                                                                                                     ggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGGCTTcac                                               >  2:386495/1‑150 (MQ=255)
                                                                                                     ggCCTAAAGCCCAGTTATTCAGTAGTCAGATAACAAAAAACCCCGAAAAATCGGGGCTTTATATAACTGACCCCTATAACCGCAACTGCGGTCTGGAGCACTATCCAGATGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTcac                                               <  1:90489/150‑1 (MQ=255)
                                                                                                         tAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGta                                           >  2:1391790/1‑150 (MQ=255)
                                                                                                             gCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAtt                                       >  1:506230/1‑150 (MQ=255)
                                                                                                                cAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGa                                    >  1:1145458/1‑150 (MQ=255)
                                                                                                                 aGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaa                                   <  2:159869/150‑1 (MQ=255)
                                                                                                                     aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTCCGAAGGCCGAAGACATAACAGTATATTTTAAAAAg                               >  2:279451/1‑150 (MQ=255)
                                                                                                                      ttCTGTAGTCATATAACAAAAAACCCCGATAAATCGTGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGTCGAAGTCTTCACAGTATATTTGAAAAAgg                              <  2:498247/150‑1 (MQ=255)
                                                                                                                      ttCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAgg                              >  2:686603/1‑150 (MQ=255)
                                                                                                                          gTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGActc                          <  2:1101054/150‑1 (MQ=255)
                                                                                                                          gTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGActc                          <  2:1660433/150‑1 (MQ=255)
                                                                                                                            aGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTTGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACTAAGGCCGACGTCGTCACAGTATATTTGAAAAATGACTCTa                        >  2:837690/1‑150 (MQ=255)
                                                                                                                               cAGATAACAAAAAACCCCGATAATTCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGACGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAgg                     <  2:1053395/150‑1 (MQ=255)
                                                                                                                               cAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAgg                     >  1:1636816/1‑150 (MQ=255)
                                                                                                                               cAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAgg                     >  2:52741/1‑150 (MQ=255)
                                                                                                                                aGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACAACGAAGGCCGAAGTCTTAACAGTATTTTTGAAAAAGGCCTCTAAggg                    >  2:1569246/1‑150 (MQ=255)
                                                                                                                                 gATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGa                   <  2:758373/150‑1 (MQ=255)
                                                                                                                                 gATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGa                   >  1:1121625/1‑150 (MQ=255)
                                                                                                                                   tAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGaaa                 >  1:553065/1‑150 (MQ=255)
                                                                                                                                    tatataaagCCCCGATTTATCGGGGTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAg                >  1:19918/6‑150 (MQ=255)
                                                                                                                                     aCAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAATGGACTCTAAGGGAAAGt               >  1:426318/1‑150 (MQ=255)
                                                                                                                                     aCAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGt               >  1:1010100/1‑150 (MQ=255)
                                                                                                                                      cAAAAAACCCCGATAAATCGTGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTg              <  1:900611/150‑1 (MQ=255)
                                                                                                                                      cAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTg              <  1:358511/150‑1 (MQ=255)
                                                                                                                                       ataaagCCCCGATTTATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGAc                            >  1:1447579/3‑135 (MQ=255)
                                                                                                                                       ataaagCCCCGATTTATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGAc                            <  2:1447579/133‑1 (MQ=255)
                                                                                                                                       aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTgg             <  1:710166/150‑1 (MQ=255)
                                                                                                                                       aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTgg             >  1:1425262/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCACGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTgg             <  2:1245242/150‑1 (MQ=255)
                                                                                                                                                  aTAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCGCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGcc  <  2:393621/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCC  >  NC_000913/3317931‑3318226

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