Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 892,457 T→A 7.4% I164N (ATT→AAT)  rimK → ribosomal protein S6 modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913892,4570TA 7.4% 265.6 / 6.9 94I164N (ATT→AAT) rimKribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base T (35/52);  new base A (3/4);  total (38/56)
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.40e-01

GCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAAT  >  NC_000913/892308‑892605
                                                                                                                                                      |                                                                                                                                                    
gcctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGAt                                                                                                                                                      <  2:601246/150‑1 (MQ=255)
 cctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGAtt                                                                                                                                                     >  2:1530843/1‑150 (MQ=255)
     tCACGGGCATTGCGCATTCGCCGGATGATACGAGCGATTTAATCGACATGGTCGGTGGTGCGCCTCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACg                                                                                                                                                 <  1:1302930/150‑1 (MQ=255)
      cACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCCGGCGGCGGAAAGCGTGATTGACGc                                                                                                                                                >  2:34439/1‑150 (MQ=255)
      cACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGc                                                                                                                                                >  1:1343619/1‑150 (MQ=255)
        cGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCtt                                                                                                                                              <  2:549177/150‑1 (MQ=255)
          ggCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTc                                                                                                                                            <  1:1402773/150‑1 (MQ=255)
           gCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTCAATCACGCTTTcc                                                                                                                                           <  1:557680/150‑1 (MQ=255)
            cATTGCTCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTCAATCACGCTTTCcg                                                                                                                                          <  2:1400175/150‑1 (MQ=255)
              ttGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTCAATCACGCTTTCcgcc                                                                                                                                        <  1:1210364/150‑2 (MQ=255)
              ttGCGCATTCGCCGGATGATACCAGCGATTTAAGCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTGTCcgcg                                                                                                                                        <  1:916400/150‑1 (MQ=255)
                gcgccgtCGCCGGATGCTTCCAGCGATTTAATCGACATGGGCGGTGGCGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCGGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGt                                                                                                                                      <  1:1525450/144‑1 (MQ=255)
                gcgcATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGt                                                                                                                                      >  2:1663438/1‑150 (MQ=255)
                  gcATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCcgccgcc                                                                                                                                     <  2:379424/150‑5 (MQ=38)
                      tCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAc                                                                                                                                <  2:540145/150‑1 (MQ=255)
                       cGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcg                                                                                                                               >  1:1200579/1‑150 (MQ=255)
                       cGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcg                                                                                                                               <  1:1185964/150‑1 (MQ=255)
                       cGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGGGCGCCGCTGGTGGTCATGGTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcg                                                                                                                               <  1:832416/150‑1 (MQ=255)
                        gtcgtatgatACCAGCGATTTATTCGACATGGTCGGTGGTGCGCCGCTGGTGTTCATGTTGGTTGAAGGCACGCAGGGATTTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcgc                                                                                                                              <  2:326382/145‑1 (MQ=255)
                              tgatACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCatat                                                                                                                        <  2:1327160/150‑1 (MQ=255)
                                  aCCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGCCGCTTTCCGCGGTCTGCACGCGCATATTCTg                                                                                                                    >  1:1071320/1‑150 (MQ=255)
                                   ccAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAACGCGTGATTGACGCTTTCCGCGGTCTGACCGCGCATATTGTgg                                                                                                                   >  1:51581/1‑150 (MQ=255)
                                   ccAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTgg                                                                                                                   <  1:1423373/150‑1 (MQ=255)
                                   ccAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTgg                                                                                                                   <  1:394811/150‑1 (MQ=255)
                                    cAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTAGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGt                                                                                                                  >  1:236051/1‑150 (MQ=255)
                                     aGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTg                                                                                                                 <  2:1301412/150‑1 (MQ=255)
                                      gCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGc                                                                                                                <  1:1315087/150‑1 (MQ=255)
                                       cGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCa                                                                                                               >  1:1280807/1‑150 (MQ=255)
                                        gATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAg                                                                                                              >  2:561561/1‑150 (MQ=255)
                                         aTTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCACGGTGGTTGAAGGCACGCAGGGACTTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAgg                                                                                                             <  1:705019/150‑1 (MQ=255)
                                         aTTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAgg                                                                                                             <  1:577647/150‑1 (MQ=255)
                                            tactCGACTTGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGGGGTCTGAACGCGCATATTCTGGTGCAGGAat                                                                                                          <  1:592727/147‑1 (MQ=255)
                                            tAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCATGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAat                                                                                                          <  2:1631773/150‑1 (MQ=255)
                                             aaTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAata                                                                                                         >  2:1688699/1‑150 (MQ=255)
                                             aaTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAata                                                                                                         <  2:1200579/150‑1 (MQ=255)
                                                 gACATGGTCGGTGGTGCGCCGCTGGTGGTCATGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATc                                                                                                     <  1:585076/150‑1 (MQ=255)
                                                 gACAGGGGCGGTGGTGCGCCGCTGGTGGTCACGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGGGCTGGCGGAGACGCGGCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATc                                                                                                     <  1:88240/150‑1 (MQ=255)
                                                  ccATGGTCGGTGGTGCGCCGCTGGTGGTCAAGGTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCa                                                                                                    <  2:1465289/149‑1 (MQ=255)
                                                       gTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAagag                                                                                               <  1:1530843/150‑1 (MQ=255)
                                                        tCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAgg                                                                                              <  2:161584/150‑1 (MQ=255)
                                                          ggtggtGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGgcg                                                                                            >  2:696692/1‑150 (MQ=255)
                                                           gtggtgCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGgcgc                                                                                           >  1:843077/1‑150 (MQ=255)
                                                             ggtgCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCaa                                                                                         >  1:1704103/1‑150 (MQ=255)
                                                                  gccgcTTGTGGTCTAGTTGGTGGATGGCTCGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCTGGCGGCGGATAGGGTGATTGACGCTTTCGGCGGTCTGTACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCTTGGGTg                                                                                    <  1:234737/150‑1 (MQ=255)
                                                                  gccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGTCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACTCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTg                                                                                    <  1:1688699/150‑1 (MQ=255)
                                                                   ccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGGCGTGCGGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGc                                                                                   <  2:1442581/150‑1 (MQ=255)
                                                                     gcTGGAGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGGCAGGCGGCGGATAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCATGGGTGCGa                                                                                 <  2:1618480/150‑1 (MQ=255)
                                                                      cTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGat                                                                                <  1:1263311/150‑1 (MQ=255)
                                                                      cTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGaa                                                                                >  1:1679989/1‑149 (MQ=255)
                                                                            gtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATTTATCAAAGAGTCGCAAGGGTGCGATATCCGc                                                                          >  2:1129789/1‑150 (MQ=255)
                                                                             tCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCt                                                                         >  2:1229605/1‑150 (MQ=255)
                                                                               aaGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGGCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGt                                                                       <  1:1288736/150‑1 (MQ=255)
                                                                                 gttggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCt                                                                     >  2:120869/1‑150 (MQ=255)
                                                                                  ttggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTg                                                                    <  1:195342/150‑1 (MQ=255)
                                                                                    ggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGCCGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                  >  1:1225830/1‑150 (MQ=255)
                                                                                    ggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                  <  1:1530489/150‑1 (MQ=255)
                                                                                    ggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                  <  2:1067350/150‑1 (MQ=255)
                                                                                       tgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttgt                                                               <  2:1343619/150‑1 (MQ=255)
                                                                                         aaGGCCCGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCGGGTGCAGGAATATATCAAAGCGGCGCAAGGGTGCGCTATCCGCTGTCGGGttgttg                                                             >  1:501517/1‑150 (MQ=255)
                                                                                         aaGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGttgttg                                                             >  2:710589/1‑150 (MQ=255)
                                                                                           ggCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGc                                                           <  2:954166/150‑1 (MQ=255)
                                                                                             cACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTGTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGa                                                         >  2:175822/1‑150 (MQ=255)
                                                                                             cACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGa                                                         >  2:490595/1‑150 (MQ=255)
                                                                                              aCGCAGGGAATTGGCGTCGTGCGGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGAt                                                        <  2:843077/150‑1 (MQ=255)
                                                                                                gCAGGGAATTCGCGTCGCGCTGGCGGAGACGCGGCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGa                                                      <  2:231557/150‑1 (MQ=255)
                                                                                                 cAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGaa                                                     <  2:999427/150‑1 (MQ=255)
                                                                                                   gggAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGt                                                   <  1:710589/150‑1 (MQ=255)
                                                                                                    ggAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTg                                                  <  1:561561/150‑1 (MQ=255)
                                                                                                        ttGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCg                                              <  1:696692/150‑1 (MQ=255)
                                                                                                              tcgtGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGa                                        <  2:1280807/150‑1 (MQ=255)
                                                                                                                gtGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGTACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGAtt                                      <  1:534286/150‑1 (MQ=255)
                                                                                                                 tGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGATGTGGTCGCTGc                                          <  1:233181/145‑1 (MQ=255)
                                                                                                                 tGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGATGTGGTCGCTGc                                          >  2:233181/1‑145 (MQ=255)
                                                                                                                 tGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTg                                     >  2:1143924/1‑150 (MQ=255)
                                                                                                                 tGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTg                                     >  2:1149599/1‑150 (MQ=255)
                                                                                                                     ggCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGACCGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAAcg                                 >  1:252678/1‑150 (MQ=255)
                                                                                                                     ggCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAAcg                                 >  2:623548/1‑150 (MQ=255)
                                                                                                                           gaCGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACggcgggc                           <  1:1335836/150‑1 (MQ=255)
                                                                                                                              gcgTCAGGCGGCGGAAAGCGTGATTGACGCGTTCCGCGGGCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGaa                        <  1:129593/150‑1 (MQ=255)
                                                                                                                                   cAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAagag                    <  1:120869/150‑1 (MQ=255)
                                                                                                                                    aGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAgg                   >  1:88149/1‑150 (MQ=255)
                                                                                                                                     ggcggcGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAggg                  <  1:1663438/150‑1 (MQ=255)
                                                                                                                                     ggcggcGGAAAGCGGGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAggg                  <  1:522355/150‑1 (MQ=255)
                                                                                                                                     ggcggaGGAAAGCGTGATTGACGCTTTCCGCTGTCTGAACGCGGATTTTCTGGTGCAGGAATATATCAATGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGGTGCGATTGAACGGCGGGCGAAAGAggg                  <  1:673042/150‑1 (MQ=255)
                                                                                                                                      gcggcgGAAAGGGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGGGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGc                 <  2:259469/150‑1 (MQ=255)
                                                                                                                                       ccgcggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCg                >  2:1108043/3‑150 (MQ=255)
                                                                                                                                          cggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAtt             >  1:1331192/1‑150 (MQ=255)
                                                                                                                                           ggAAAGCGTGATTGACGCTTTCCGCGGTCTGACCGCGCATATTCTGGTGCAGGAATATATAAAAGATTAACAAGGGTGAGATATACGGTGTCTGGTTGGTGGCGAGGAAGTGGTCGCTGAGATTCAACCGCGCGCGAAAGCggtagagtt            >  2:297051/1‑143 (MQ=255)
                                                                                                                                           ggAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAttt            >  2:187242/1‑150 (MQ=255)
                                                                                                                                           ggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAttt            >  2:681248/1‑150 (MQ=255)
                                                                                                                                            gAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAtttt           >  1:145738/1‑150 (MQ=255)
                                                                                                                                             aaaGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTc          >  1:813800/1‑150 (MQ=255)
                                                                                                                                                gCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGtt       <  1:969966/150‑1 (MQ=255)
                                                                                                                                                 cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTc      >  1:47374/1‑150 (MQ=255)
                                                                                                                                                 cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTc      >  1:583581/1‑150 (MQ=255)
                                                                                                                                                     aTTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAAt  >  1:737566/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAAT  >  NC_000913/892308‑892605

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: