Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,053 (T)9→8 100% intergenic (+47/‑80) rlmB → / → yjfI 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,0450T.85.6% 324.7 / 60.9 104intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (11/4);  new base . (43/46);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 9.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 126 total aligned reads displayed.

CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATG  >  NC_000913/4409946‑4410140
                                                                                                   |                                                                                               
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCAGTGCGCCAGCGAAGCTCATTTCTCAGAATTATGAAAAGCCATCCAGACTTGGATGGt                                                                                                 >  1:1730639/1‑100 (MQ=255)
   gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAACTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttt                                                                                              >  1:2896141/1‑100 (MQ=255)
    aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGctttt                                                                                             <  2:2637762/100‑1 (MQ=38)
     aCGTCTCGGTTGCGACCGTAATTTGCTTATGTGAAGCGGTGCGCCAGCGCAGCTAATGTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGcttt                                                                                              <  1:1112856/98‑1 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           >  1:75577/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           >  2:1344322/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGCttttttt                                                                                           <  2:1525967/100‑1 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTGAGAATTATGAAAGGCCATCCAGATATGGATGGtttttttt                                                                                          >  2:1166183/1‑100 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                          >  2:1340846/1‑100 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                          <  1:2890864/100‑4 (MQ=37)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                          <  2:444598/100‑4 (MQ=37)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                          >  1:2844971/1‑99 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcata                                                                                         <  2:154793/100‑5 (MQ=37)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       >  1:1448541/1‑97 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  2:1479007/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  2:1862267/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                          <  1:2853799/97‑4 (MQ=37)
          tCGGTTCCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataat                                                                                       <  2:1914406/100‑7 (MQ=37)
           cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGTATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                      <  1:611246/100‑5 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGATTTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCta                                                                                    <  1:2782806/100‑1 (MQ=255)
               tGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAg                                                                                 >  2:373474/1‑100 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGc                                                                                <  1:782240/100‑1 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGc                                                                                <  1:782258/100‑1 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAATTGGATGG‑TTTTTTTTGTCTATAGCt                                                                               >  2:1411060/1‑100 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGt                                                                            >  2:13881/1‑100 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGt                                                                            <  2:592160/100‑1 (MQ=255)
                      ggAGTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTg                                                                              <  1:896197/96‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTaa                                                                          <  2:2262845/100‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTaa                                                                          >  1:1979514/1‑100 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑CTTTTTTTGTCTATAGCTGGTaa                                                                          <  1:2092774/100‑1 (MQ=255)
                         aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGa                                                                        >  2:2384788/1‑100 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                      <  2:289473/100‑1 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                      <  1:1579534/100‑1 (MQ=255)
                             gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                   >  1:1900623/1‑100 (MQ=255)
                               ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                                 <  1:2795377/100‑1 (MQ=255)
                               ttATTTGAAGCGCTGCGCCAGCCCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCCTTAGACGG‑TTCCTTTTGTCTATAGCTGCTAAGATAAGTAc                                                                 >  1:657893/1‑100 (MQ=255)
                                tATTTGAAGCTGTGCGCCAACGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTAGATGG‑TTTTTTTTGTCTATAGCTGGTAATATAATTAc                                                                 <  1:2596168/99‑1 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                >  2:565354/1‑100 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCTTTGGATGGTTTTTTTTTGTCTTTTGCTGGTAAGATAATTACg                                                                >  2:969432/1‑100 (MQ=255)
                                  tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                              >  2:1611185/1‑100 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                            >  2:2611892/1‑100 (MQ=255)
                                         cGGTGCGCCAGCGCAGCTAATATCTCAGAATCATGAACAGACATCCAGATTTGGATGG‑TCTTTTTTGTCTAAAGCAGCTAAGATAATTACGTATTGCa                                                         >  1:1286234/1‑98 (MQ=255)
                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                      >  2:2272792/1‑100 (MQ=255)
                                              cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                   >  2:1156936/1‑100 (MQ=255)
                                               gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                 <  2:747728/100‑1 (MQ=255)
                                                   gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                             <  2:2305118/100‑1 (MQ=255)
                                                   gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                             <  2:1348873/100‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTtta                                            <  1:969430/99‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                          <  2:867112/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                         >  2:2332715/1‑100 (MQ=255)
                                                        gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                        >  2:1240312/1‑100 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                      >  2:1976284/1‑100 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                      >  2:2059533/1‑100 (MQ=255)
                                                           aaTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAatgt                                     >  1:1758308/1‑98 (MQ=255)
                                                             tttCTCAGAATTATGAAAAGCCATACAGAATTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTCTTAATAtt                                    >  2:2878603/1‑100 (MQ=255)
                                                              ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                  <  2:248259/100‑1 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTATGTATTTCAAATATTCATTTTATAAATAttt                                   >  1:2558912/1‑97 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                <  2:1872101/100‑1 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                >  1:1410314/1‑99 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                               <  2:1565168/100‑1 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                             >  2:883720/1‑100 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                             >  2:424444/1‑100 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTAAGTATTGCAAATATTCATTTTATTAATATTTAAACtt                             >  1:2074977/1‑100 (MQ=255)
                                                                      gacaaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          >  2:326122/7‑100 (MQ=37)
                                                                      aTTATGAAAAGCCATCCAGATTTGTAT‑ATTTTTTTTTGTCTATAGTTGCTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          <  1:2766746/100‑1 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         >  1:417005/1‑100 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATACATTTTATTAATATTTAAACTTAATc                         <  2:1776608/100‑1 (MQ=255)
                                                                        tATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                        <  1:102353/100‑1 (MQ=255)
                                                                         aTGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                       <  1:743074/100‑1 (MQ=255)
                                                                          tGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTTTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                      >  1:418227/1‑100 (MQ=255)
                                                                          aaaaaaACCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                      >  2:2843438/3‑100 (MQ=255)
                                                                          aaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                      >  1:1998618/3‑100 (MQ=38)
                                                                          aaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTAAAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                      >  1:401326/3‑100 (MQ=37)
                                                                            aaaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                    <  2:1349821/100‑1 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGGTTTTTTTTT‑TTTTTTTTTTGTAAGATAATTACGTATTGCAAATATTCATTTTTTTAATATTTAAACTTTATCTTTATg                   >  1:475914/1‑100 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTTTATAGCTGGTAAAATAATTACGTATTATAAATATTTATTTTATTAATATTTATACTTACTatatat                    >  1:1097040/1‑99 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  1:2868175/100‑1 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  2:1375101/100‑1 (MQ=255)
                                                                               aGCCATCCAAATTTGGATGG‑TTTTTTTTGTCTATAGTTTGTAAGATAATTAAGTTTTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                  <  1:603156/99‑1 (MQ=255)
                                                                                gCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAATTAAACTTAATCTATATTGaa                >  1:2154150/1‑100 (MQ=255)
                                                                                accATCCAAATCTGGATGG‑TTTTTTTTGTCTATCGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                >  1:1746518/2‑100 (MQ=255)
                                                                                accATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                  >  2:1055912/2‑98 (MQ=255)
                                                                                accATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                >  2:2668727/2‑100 (MQ=255)
                                                                                 ccATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat               <  2:1740296/100‑1 (MQ=255)
                                                                                 ccATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat               >  2:2367324/1‑100 (MQ=255)
                                                                                 ccATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat               <  2:465408/100‑1 (MQ=255)
                                                                                  cATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaata              <  2:2271902/100‑1 (MQ=255)
                                                                                   aTCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataa             <  1:1344320/100‑1 (MQ=255)
                                                                                    tCCAGATTTGGATGGTTTTTTTTT‑TTTTTTTTTTGTAAGATAATTACGTATTGCAATTATTCATTTTATTAATATTTAAACTTAATCTTTTTGGaat               >  1:2139527/1‑97 (MQ=255)
                                                                                      cAGATTTGGATGG‑TTTTTTTTGTCTATAGTTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           <  1:635475/99‑1 (MQ=255)
                                                                                       aGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata         <  1:1001245/100‑1 (MQ=255)
                                                                                        gATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat        >  2:1913288/1‑100 (MQ=255)
                                                                                         aTTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg       <  1:2367313/100‑1 (MQ=255)
                                                                                         aTTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg       <  2:659032/100‑1 (MQ=255)
                                                                                          tttGGATGGTTTTTTTTT‑TTTTTTTCTGGTACGCTAATTACGTACTGCAAATATTCTTTTTATTAATATTTACCCCTACTCTATCTGTAATACTATATg       >  1:219137/1‑99 (MQ=255)
                                                                                             ggATGG‑TTTTTTTTGTCTATAGCTTGTAATATAATTAAGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   <  1:1842281/100‑1 (MQ=255)
                                                                                              gATGGTTG‑TTTTTGTTTATAGCTGGTAAGATAATTAAGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAAGATATGACATg  <  1:2186934/100‑1 (MQ=255)
                                                                                              gATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg  <  1:2186909/100‑1 (MQ=255)
                                                                                              gATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg  <  2:1194077/100‑1 (MQ=255)
                                                                                              gATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg  <  1:266181/100‑1 (MQ=255)
                                                                                                   |                                                                                               
CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATG  >  NC_000913/4409946‑4410140

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑