Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,449,388 C→T E166K (GAA→AAA)  rpsC ← 30S ribosomal subunit protein S3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,449,3880CT98.7% 302.4 / ‑4.1 79E166K (GAA→AAA) rpsC30S ribosomal subunit protein S3
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/47);  total (31/48)
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.84e-01

CAACAGCAGCCATACCACCCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTCGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTGATGCTGTCAGC  >  NC_000913/3449241‑3449535
                                                                                                                                                   |                                                                                                                                                   
cAACAGCAGCCATACCACCCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTgg                                                                                                                                                   >  1:829081/1‑150 (MQ=255)
  aCAGCAGCCATACCACCCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTAc                                                                                                                                                >  1:1183201/1‑151 (MQ=255)
    agcagcCATACCACCCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGt                                                                                                                                              >  1:704431/1‑151 (MQ=255)
    agcagcCATACCACCCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGt                                                                                                                                              <  2:961602/151‑1 (MQ=255)
     gcagcCATACCACCCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTg                                                                                                                                             <  2:168504/151‑1 (MQ=255)
             accaccCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTcc                                                                                                                                     <  1:923639/151‑1 (MQ=255)
                accCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCcgc                                                                                                                                   >  1:922698/1‑150 (MQ=255)
                     ggATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCgccgcc                                                                                                                             >  2:413371/1‑151 (MQ=255)
                      gATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGccc                                                                                                                            >  2:1086180/1‑151 (MQ=255)
                         ctcGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGccc                                                                                                                            >  2:485029/1‑148 (MQ=255)
                         ctcGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGccc                                                                                                                            >  2:788624/1‑148 (MQ=255)
                           cGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCa                                                                                                                           <  2:724768/147‑1 (MQ=255)
                           cGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAg                                                                                                                          <  1:785101/148‑1 (MQ=255)
                           cGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAg                                                                                                                          <  2:1003680/148‑1 (MQ=255)
                                ttGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCg                                                                                                                   <  2:321864/150‑1 (MQ=255)
                                  gAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCt                                                                                                                 <  2:805745/150‑1 (MQ=255)
                                   aaGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTaa                                                                                                               <  2:858799/151‑1 (MQ=255)
                                   aaGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTa                                                                                                                >  1:1196494/1‑150 (MQ=255)
                                    aGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAAc                                                                                                              >  1:1068503/1‑151 (MQ=255)
                                     gATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACt                                                                                                             >  2:1121513/1‑151 (MQ=255)
                                     gATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACt                                                                                                             <  1:777590/151‑1 (MQ=255)
                                        ccACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTc                                                                                                           <  1:1038965/150‑1 (MQ=255)
                                         cacaCTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCa                                                                                                          >  1:585626/1‑150 (MQ=255)
                                            acTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACtt                                                                                                      <  1:168124/151‑1 (MQ=255)
                                                tAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTaa                                                                                                   <  2:578675/150‑1 (MQ=255)
                                                 aaCGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAAt                                                                                                  <  1:560883/150‑1 (MQ=255)
                                                    gCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCt                                                                                              <  2:1134985/151‑1 (MQ=255)
                                                    gCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCt                                                                                              >  2:392624/1‑151 (MQ=255)
                                                      cGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCttt                                                                                            <  1:201715/151‑1 (MQ=255)
                                                      cGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCttt                                                                                            <  1:271344/151‑1 (MQ=255)
                                                         ttccacCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGTCCTTCGCGGTACCATTTGGTACGTGCGATTTCCGAGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAg                                                                                          <  2:735807/147‑1 (MQ=255)
                                                         ttACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAgc                                                                                         >  2:492865/1‑151 (MQ=255)
                                                          tACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAgc                                                                                         <  1:1081490/150‑1 (MQ=255)
                                                              ccGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGccc                                                                                     <  2:203816/150‑1 (MQ=255)
                                                              ccGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCa                                                                                    <  1:1167738/151‑1 (MQ=255)
                                                                gTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGccc                                                                                     >  1:633775/1‑148 (MQ=255)
                                                                gTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGCACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGccc                                                                                     >  1:538324/1‑148 (MQ=255)
                                                                    gtggtgTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGc                                                                               >  1:629358/1‑150 (MQ=255)
                                                                      ggtgTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCAtt                                                                            <  1:676553/151‑1 (MQ=255)
                                                                         gtgCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCg                                                                         <  2:842556/151‑1 (MQ=255)
                                                                           gcgcTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGtt                                                                       <  1:824441/151‑1 (MQ=255)
                                                                            cgcTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTc                                                                      <  1:939161/151‑1 (MQ=255)
                                                                              cTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTg                                                                    >  1:712721/1‑151 (MQ=255)
                                                                                tCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGt                                                                   <  1:843118/150‑1 (MQ=255)
                                                                                 cAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTAc                                                                 <  2:932454/151‑1 (MQ=255)
                                                                                  agagGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGAACCATTTGGTACGGGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTAc                                                                 >  1:539949/1‑150 (MQ=255)
                                                                                    agGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGc                                                              >  1:782577/1‑151 (MQ=255)
                                                                                        gttgtAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACg                                                           <  1:893274/150‑1 (MQ=255)
                                                                                        gttgtAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACg                                                           <  2:294960/150‑1 (MQ=255)
                                                                                             aGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCAt                                                     >  2:991897/1‑151 (MQ=255)
                                                                                              gTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATa                                                    <  2:737638/151‑1 (MQ=255)
                                                                                                cGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGc                                                  <  2:970624/151‑1 (MQ=255)
                                                                                                       aGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACgg                                           <  1:468463/151‑1 (MQ=255)
                                                                                                        gCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACgg                                           >  2:1042524/1‑150 (MQ=255)
                                                                                                        gCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGa                                          <  2:828939/151‑1 (MQ=255)
                                                                                                          aCGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAAc                                        <  2:1204408/151‑1 (MQ=255)
                                                                                                            gCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACAt                                      <  2:505667/151‑1 (MQ=255)
                                                                                                            gCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACAt                                      >  1:803863/1‑151 (MQ=255)
                                                                                                             cAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATa                                     <  1:950183/151‑1 (MQ=255)
                                                                                                                 gtgtGCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAAcgc                                 <  1:1181344/151‑1 (MQ=255)
                                                                                                                  tgtgCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAAcgcg                                >  1:273385/1‑151 (MQ=255)
                                                                                                                   gtgCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGTACGACGGAACATAACGCGa                               <  2:633658/151‑1 (MQ=255)
                                                                                                                   gtgCAGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGa                               <  2:126127/151‑1 (MQ=255)
                                                                                                                   gtgCAGCGGTACGCGACCTTCGCGGTACCATTCGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAAcgcg                                <  1:589909/150‑1 (MQ=255)
                                                                                                                       aGCGGTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGtt                           >  1:153466/1‑151 (MQ=255)
                                                                                                                           gTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTcc                         >  1:211530/1‑149 (MQ=255)
                                                                                                                           gTACGCGACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTcc                         >  2:381964/1‑149 (MQ=255)
                                                                                                                               gcgACCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTg                    >  2:21923/1‑150 (MQ=255)
                                                                                                                                  aCCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTgaga                 >  1:166124/1‑150 (MQ=255)
                                                                                                                                  aCCTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGaa                <  2:920929/151‑1 (MQ=255)
                                                                                                                                   ccTTCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAg               <  2:385711/151‑1 (MQ=255)
                                                                                                                                      tCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGt              <  2:1218346/149‑1 (MQ=255)
                                                                                                                                      tCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTg             <  1:231271/150‑1 (MQ=255)
                                                                                                                                      tCGCGGTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTg             <  2:843912/150‑1 (MQ=255)
                                                                                                                                          ggTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTGATGCt        <  2:1219526/151‑1 (MQ=255)
                                                                                                                                           gTACCATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTGATGCTg       <  1:824862/151‑1 (MQ=255)
                                                                                                                                               cATTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTGATGCTGTCAg   <  1:1032629/151‑1 (MQ=255)
                                                                                                                                                aTTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTGATGCTGTCAGc  >  1:307441/1‑151 (MQ=255)
                                                                                                                                                aTTTGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTGATGCTGTCAGc  >  2:13817/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAACAGCAGCCATACCACCCAGGATCTCGCCTTTGAAGATCCACACTTTAACGCCGATTACACCGTAAGTGGTGTGCGCTTCAGAGGTGTTGTAGTCGATGTCAGCACGCAGAGTGTGCAGCGGTACGCGACCTTCGCGGTACCATTCGGTACGTGCGATTTCCGCGCCGCCCAGACGGCCGCTAACTTCAACTTTAATACCTTTAGCGCCCAGACGCATTGCGTTCTGTACAGCACGCTTCATAGCACGACGGAACATAACGCGACGTTCCAGCTGAGAAGTGATGCTGTCAGC  >  NC_000913/3449241‑3449535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑