Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 59,968 C→T V127I (GTC→ATC)  rluA ← dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091359,9680CT96.8% 1255.7 / 18.2 378V127I (GTC→ATC) rluAdual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (8/4);  new base T (194/170);  total (203/175)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 385 total aligned reads displayed.

AGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGG  >  NC_000913/59820‑60114
                                                                                                                                                    |                                                                                                                                                  
agaATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGAttt                                                                                                                                                  <  7:817790/151‑1 (MQ=255)
agaATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGAttt                                                                                                                                                  <  6:946202/151‑1 (MQ=255)
agaATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGAttt                                                                                                                                                  >  6:1638379/1‑151 (MQ=255)
agaATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTGTTAACACCACTCTTGCCGTTTTATCCGCCGCATACTCAACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGAttt                                                                                                                                                  <  5:1397639/151‑1 (MQ=255)
   aTCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCt                                                                                                                                               >  5:1202035/1‑151 (MQ=255)
      ggATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGGGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTGtt                                                                                                                                            <  6:637823/151‑1 (MQ=255)
         ttaccccccgcCAGCATTTGCACACGCTATTGAGGCGAACGGCCGGTATTCGGTTTTAACAGCACTCTTGCCGTGTTATCCGCCGGATCCTGCACCACTTCTTATTACGTCTGAGCAGGTTTAGCCGTTTCTTATCTGTCTTTCAGttta                                                                                                                                          <  4:1010237/142‑2 (MQ=255)
         tGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCg                                                                                                                                         >  6:782843/1‑151 (MQ=255)
          gACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCgg                                                                                                                                        <  7:681205/151‑1 (MQ=255)
          gACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCgg                                                                                                                                        >  7:1644544/1‑151 (MQ=255)
             ccAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCg                                                                                                                                     >  5:15354/1‑151 (MQ=255)
             ccAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCg                                                                                                                                     >  5:1311878/1‑151 (MQ=255)
             ccAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCg                                                                                                                                     >  3:604140/1‑151 (MQ=255)
             ccAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCg                                                                                                                                     <  2:1247962/151‑1 (MQ=255)
              cAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTGTTTCGGGCgg                                                                                                                                    >  8:920865/1‑151 (MQ=255)
                   ccAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTgg                                                                                                                               >  7:1253932/1‑151 (MQ=255)
                   ccAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTgg                                                                                                                               <  1:1023132/151‑1 (MQ=255)
                    cAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGc                                                                                                                              >  5:1674812/1‑151 (MQ=255)
                    cAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGc                                                                                                                              >  7:945305/1‑151 (MQ=255)
                    cAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGc                                                                                                                              <  6:874013/151‑1 (MQ=255)
                      gCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGcc                                                                                                                             >  4:1475004/1‑150 (MQ=255)
                       cATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAg                                                                                                                           <  6:442243/151‑1 (MQ=255)
                         tatGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTc                                                                                                                         >  4:782822/1‑151 (MQ=255)
                            gCACACGCAATTGATGCGAACGCCCGCTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCCTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCAAGTCGc                                                                                                                       >  8:1399252/1‑150 (MQ=255)
                              acacGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCaaa                                                                                                                    <  7:236907/151‑1 (MQ=255)
                              acacGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCaa                                                                                                                     >  1:123019/1‑150 (MQ=255)
                               cacGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAAt                                                                                                                   <  1:1502019/151‑1 (MQ=255)
                               cacGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAAt                                                                                                                   >  8:1005911/1‑151 (MQ=255)
                                acGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATc                                                                                                                  <  8:1535720/151‑1 (MQ=255)
                                acGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATc                                                                                                                  <  8:391290/151‑1 (MQ=255)
                                   cAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGc                                                                                                               <  8:1281292/151‑1 (MQ=255)
                                    aaTTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCg                                                                                                              >  3:736093/1‑151 (MQ=255)
                                     aTTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCgg                                                                                                             >  5:6046/1‑151 (MQ=255)
                                       tGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGc                                                                                                            <  5:371835/150‑1 (MQ=255)
                                        gATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCa                                                                                                           <  5:1101075/150‑1 (MQ=255)
                                           gCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGAt                                                                                                        <  8:191307/150‑1 (MQ=255)
                                           gCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGAt                                                                                                        <  1:1637990/150‑1 (MQ=255)
                                            cGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATcc                                                                                                      >  1:474487/1‑151 (MQ=255)
                                               aCGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATccacc                                                                                                   <  6:851064/151‑1 (MQ=255)
                                                cGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATccacc                                                                                                   >  8:35052/1‑150 (MQ=255)
                                                 cccccGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAg                                                                                                 <  4:864041/150‑1 (MQ=255)
                                                  cccGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATccacc                                                                                                   >  8:611798/1‑148 (MQ=255)
                                                     ggTAATCGGTTTTAAGACCACTCTTGCCGTGTTATCAGCCGCATACTCCACCACTTCATATTCGGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCGAGTCGCAAATCAGCGGCAGATCCACCAGACCt                                                                                             <  8:119193/151‑1 (MQ=255)
                                                     ggTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCt                                                                                             <  6:1375469/151‑1 (MQ=255)
                                                       tAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTc                                                                                           >  7:1097447/1‑151 (MQ=255)
                                                        aaTCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCGGTCGCAAATCAGCGGCAGATCCACCAGACCTTc                                                                                           <  4:520065/150‑1 (MQ=255)
                                                        aaTCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTc                                                                                           <  3:1264794/150‑1 (MQ=255)
                                                         aTCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTg                                                                                         <  3:666785/151‑1 (MQ=255)
                                                           cGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCg                                                                                       >  4:196072/1‑151 (MQ=255)
                                                            ggTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCgg                                                                                      >  1:1108124/1‑151 (MQ=255)
                                                             gTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCggg                                                                                     >  4:1396273/1‑151 (MQ=255)
                                                               tttAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCggggg                                                                                   >  3:294962/1‑151 (MQ=255)
                                                                ttAACACCTCTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTGATATTCCGTCTGCGCAGGTATACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCggggg                                                                                   <  6:1079921/150‑1 (MQ=255)
                                                                ttAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCggggg                                                                                   >  7:1425902/1‑150 (MQ=255)
                                                                  aaCACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGggatg                                                                                <  8:75659/151‑1 (MQ=255)
                                                                   acacCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGggatgg                                                                               <  4:284833/151‑1 (MQ=255)
                                                                              gCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAgc                                                                    >  2:807404/1‑150 (MQ=255)
                                                                                cGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGAcg                                                                  <  1:1383097/151‑1 (MQ=255)
                                                                                   gtTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCgg                                                               <  3:1138636/151‑1 (MQ=255)
                                                                                     tATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCggg                                                              <  5:896859/150‑1 (MQ=255)
                                                                                      aTCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGcc                                                            <  4:684024/151‑1 (MQ=255)
                                                                                       tCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGccc                                                           >  7:1548376/1‑150 (MQ=255)
                                                                                          gcctcATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTTGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACAt                                                        <  4:1231534/151‑1 (MQ=255)
                                                                                          gccgcATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCCGCAGATCCACCAGACCTTCTGCGGGCGATGGATGCCCCCAGACGCGGGCAACAt                                                        >  5:696556/1‑151 (MQ=255)
                                                                                             gcATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACt                                                     >  2:548617/1‑151 (MQ=255)
                                                                                             gcATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACt                                                     >  2:926038/1‑151 (MQ=255)
                                                                                                tACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCCCCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCt                                                  >  2:1282744/1‑151 (MQ=255)
                                                                                                 aCTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCt                                                  <  3:1394502/150‑1 (MQ=255)
                                                                                                    ccaccaCTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCttttt                                              <  5:260836/151‑1 (MQ=255)
                                                                                                      cacacTTCATATTGCGTCTGCGCAGTTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCGACCAGACCTTCTGCGGTGGATGGATGCCCCCAGACGGGGGCCACATACTGCTTTTTCg                                            <  1:227840/149‑1 (MQ=255)
                                                                                                      accacTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCg                                            <  2:1077831/151‑1 (MQ=255)
                                                                                                       ccacTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCgg                                           <  6:542720/151‑1 (MQ=255)
                                                                                                       ccacTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCgg                                           >  2:1539482/1‑151 (MQ=255)
                                                                                                       ccacTTCATATTCCGTCGGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCgg                                           <  4:1600867/151‑1 (MQ=255)
                                                                                                             cATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGgctcgcg                                     <  7:1302625/151‑1 (MQ=255)
                                                                                                               tatTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGgctcgcgcg                                   >  2:395534/1‑150 (MQ=255)
                                                                                                                atTCCGTCTGCGCAGGTTTACGCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGgctcgcgctc                                  <  3:465917/151‑1 (MQ=255)
                                                                                                                  tCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGgctcgcgctcgc                                >  4:1081424/1‑151 (MQ=255)
                                                                                                                  tCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGgctcgcgctcgc                                >  1:919762/1‑151 (MQ=255)
                                                                                                                   ccGTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGctcgcgctcgcg                               >  6:1578125/1‑151 (MQ=255)
                                                                                                                     gTCTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGa                             <  2:628515/151‑1 (MQ=255)
                                                                                                                       cTGCGCAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAAc                           >  7:375811/1‑151 (MQ=255)
                                                                                                                         gcgcAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAAc                           <  1:942278/149‑1 (MQ=255)
                                                                                                                         gcgcAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTg                         <  4:427175/151‑1 (MQ=255)
                                                                                                                           gcAGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGc                       <  2:296667/151‑1 (MQ=255)
                                                                                                                             atgTTTACCCGTTGCGTAACAGATTTTCTGTTTCGGGAGGTTTGTCCAGTCGCAAATAAGCGGGAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCACAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGt                     >  6:723126/3‑151 (MQ=255)
                                                                                                                             aGGTTTACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGt                     >  1:1574233/1‑151 (MQ=255)
                                                                                                                                tttACCCGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGtttt                  >  6:1316797/1‑151 (MQ=255)
                                                                                                                                tttACCCGTTTCGTAACAGATTCTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGtttt                  >  1:753977/1‑151 (MQ=255)
                                                                                                                                   ccccGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGGTTTAAc               >  5:1445019/2‑151 (MQ=255)
                                                                                                                                    cccGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAAc               <  6:242751/150‑1 (MQ=255)
                                                                                                                                    cccGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACt              >  3:162869/1‑151 (MQ=255)
                                                                                                                                     ccGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTc             <  1:64839/151‑1 (MQ=255)
                                                                                                                                      cGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTcc            <  3:21739/151‑1 (MQ=255)
                                                                                                                                      cGTTTCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTcc            >  1:434213/1‑151 (MQ=255)
                                                                                                                                          tCGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGc         <  8:1601389/150‑1 (MQ=255)
                                                                                                                                           cGTAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTc       <  3:351093/151‑1 (MQ=255)
                                                                                                                                             tAACAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCcg     >  8:538229/1‑151 (MQ=255)
                                                                                                                                              aaCAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCcgc    >  3:1502965/1‑151 (MQ=255)
                                                                                                                                              aaCAGATTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCcgc    >  1:857397/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGAATCGGATGACCCAGCGCCAGCATATGCACACGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTGTTTCGGGCGGTTTGGCCAGTCGCAAATCAGCGGCAGATCCACCAGACCTTCTGCGGGGGATGGATGCCCCCAGACGCGGGCCACATACTGCTTTTTCGGCTCGCGCTCGCGGAACTGGCGTTTTAACTCCCGCTCCGCGG  >  NC_000913/59820‑60114

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