Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,507,867 G→A E363K (GAG→AAG)  rlhA → 23S rRNA 5‑hydroxycytidine C2501 synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,507,8670GA97.1% 104.9 / ‑5.9 35E363K (GAG→AAG) rlhA23S rRNA 5‑hydroxycytidine C2501 synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/15);  minor base C (1/0);  total (20/15)
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.41e-01

GAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCG  >  NC_000913/1507621‑1508113
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
gAACATTTCTTTGTTCCATCGCCGGAAACGACTTTCCACCGTGGTAGCACAGCTTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGACGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTCCCGAGCCAATGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCaaga                                                                                                                                                                                                                                                     >  2:596701/1‑250 (MQ=255)
 aaCATTTCTTTGTACCAGCGACGGAAAAGAATTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATTTTGGCGCGTGCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGaaa                                                                                                                                                                                                                                                   <  2:330171/251‑1 (MQ=255)
     tttCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCg                                                                                                                                                                                                                                               >  1:248103/1‑251 (MQ=255)
           tGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGaa                                                                                                                                                                                                                                          <  1:508453/250‑1 (MQ=255)
                                  tCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGCGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAAt                                                                                                                                                                                                                  >  2:475662/1‑251 (MQ=255)
                                       cGTGGTAGCACAGATTATTTTGTGAATGCCGGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAAt                                                                                                                                                                                                             <  1:321388/251‑1 (MQ=255)
                                            tAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGGCTGGCCCAATGAAATGCca                                                                                                                                                                                                         >  2:208308/1‑250 (MQ=255)
                                                          tggtgAATACCCGTACAGGCGATATTGGCGTGTTGGTTTTGCCAAAAGTTTTC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGGGGCGAAAGATCATCGCGAGGTTGCCTTTCCAGAGCCACTGGCAAATGGCGATGGCCTGAACGAGTTGATTAAACGGGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGGTCCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACaa                                                                                                                                                                                           <  2:338000/248‑1 (MQ=255)
                                                                   cccGTCAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTcc                                                                                                                                                                                  >  2:441808/1‑250 (MQ=255)
                                                                         aaGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAACGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGGCTGGCCCAATGAAATGCCCGCCGATTTGCACAAAATTAGTCCCCATCAc                                                                                                                                                                           >  2:37750/1‑251 (MQ=255)
                                                                                   tGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGAAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGACACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAACACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAAcc                                                                                                                                                                 >  1:307047/1‑251 (MQ=255)
                                                                                    gtggcgTTCGATTCGC‑GAAAATTAACTGGCTTGCCGGTAGGCGAAGAAGTGAAAGTGGCGAAAGATCTCCTCGATGTGGCCGTTACAGAGCCACTGGCAAATGGCGATGGCCTGAACGTTGTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCg                                                                                                                                                                <  2:83905/248‑1 (MQ=255)
                                                                                         gTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGCAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGCTTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTAAAGAAACCCGGAGCAAATCAGTACCGCGTCTGGCCCAATGGAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATc                                                                                                                                                           >  2:458344/1‑251 (MQ=255)
                                                                                         gTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATTGAGATGGCCTGAACGTGGTGATTAAACGTGAAGACGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACAGTAATc                                                                                                                                                           >  1:462467/1‑251 (MQ=255)
                                                                                             gtttCGCCGAAATTTACA‑GGCCTGCCGGTAGGCGAAGGAGTGAAGGTGGGGAAAGATCATCTCGATGTTGCCTTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACAGCACCCACTAAACCGTAATCTTGa                                                                                                                                                       <  2:210180/249‑1 (MQ=255)
                                                                                                 cGCCGAAATTTATC‑GTCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGTCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCa                                                                                                                                                    <  1:571048/250‑1 (MQ=255)
                                                                                                       aaTTTCTC‑GGCCTGCCAGTCGGCGAAGTCGTAACAGTGGCGAAAGATCATCTCGATGATGAAGTTACCGAAACACAGGCAACTGGCGATGGCCTGAACGTGTTGATTAAACGTGAAAACGAAGGTATTCGTGCCACTACGGTCCAGAAAACCGGAGACAATCAGTACCCAGTCTGGCCCATTGAAATGCCAGCAGATTTGCACAAAATTAGACCACATCCCCCACTAAACCGAAATATTGCTCATAACTgg                                                                                                                                             >  2:580904/1‑251 (MQ=255)
                                                                                                           tATC‑GGCCTGCAGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGAGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAACTCAGTACCGCTTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCCTCACTGgccg                                                                                                                                          >  2:23513/1‑248 (MQ=255)
                                                                                                                  ccTGCCGGAAGGCGACGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGAGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCAc                                                                                                                                   >  2:352734/1‑251 (MQ=255)
                                                                                                                           aGGCGAATTAGTGAAAGGGGCGAAAGATCATCTCGATTTTGCCGTCTCCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAACCTTGATCATAACTGGCAGCAGGCACTGACaaaaa                                                                                                                          <  2:457124/251‑1 (MQ=255)
                                                                                                                                  gtagtGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAg                                                                                                                   <  1:149143/251‑1 (MQ=255)
                                                                                                                                                     atcatcTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGt                                                                                                 >  2:66296/1‑250 (MQ=255)
                                                                                                                                                           gcGATGGTGCCGTTACCGAGCTACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTGCGTGCCAATACGGTCAAGAAAACCGGAGAAAAGCAGTTCCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACAGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGGGGCGGTAGACAtt                                                                                          <  2:191482/250‑1 (MQ=255)
                                                                                                                                                            cGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGAGTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACAtt                                                                                          <  1:533550/250‑1 (MQ=255)
                                                                                                                                                                     cGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGggcg                                                                                <  1:37750/251‑1 (MQ=255)
                                                                                                                                                                                              gctGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGAGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTATGAc                                                       >  1:100489/3‑251 (MQ=255)
                                                                                                                                                                                                   ccTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTc                                                   >  1:335829/1‑250 (MQ=255)
                                                                                                                                                                                                       aaCGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAGACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAg                                              <  1:121806/251‑1 (MQ=255)
                                                                                                                                                                                                        aCGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGt                                             >  1:496228/1‑251 (MQ=255)
                                                                                                                                                                                                           tgtTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGa                                           >  1:231242/1‑250 (MQ=255)
                                                                                                                                                                                                                   aaaCGTGAAGTCGTCGGGTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTCCCGCGTCGGGCCCAAGGAAATGCCAGCAGATTTGCACAAAATGCGTCCACATCACCCACTAACCCGTAATCTTGATCATAACTGGCAGCAGGCAATGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAggtta                                   <  2:14428/250‑4 (MQ=255)
                                                                                                                                                                                                                      cGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGAAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGgtgac                                 >  2:336266/1‑247 (MQ=255)
                                                                                                                                                                                                                                  ggTTTTCGTGCCAATACGGTCAAGCAACCCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGCCCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATAc                   >  2:435500/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cGTGCCAATACGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTgg              <  1:66296/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGGTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGGGTCAGCATCACGCATACGCTGGACGGGCAGTTc   >  2:306378/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ggTCAAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTcg  >  1:190507/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
GAACATTTCTTTGTTCCATCGACGGAAAAGACTTTCCACCGTGGTAGCACAGATTATTTTGTGAATGCCCGTAAAGGCGATATTGGCGCGTTCGATTCGCCGAAATTTATC‑GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCG  >  NC_000913/1507621‑1508113

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