Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,891,131 A→G 5.2% L46P (CTG→CCG)  tsaB ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,891,1310AG 5.2% 223.2 / 2.8 78L46P (CTG→CCG) tsaBtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD
Reads supporting (aligned to +/- strand):  ref base A (31/40);  new base G (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

AGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGC  >  NC_000913/1890984‑1891279
                                                                                                                                                   |                                                                                                                                                    
aGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAgg                                                                                                                                                    <  2:335354/150‑1 (MQ=255)
aGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAgg                                                                                                                                                    <  1:288628/150‑1 (MQ=255)
aGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCACTTGAGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAgg                                                                                                                                                    <  2:1290283/150‑1 (MQ=255)
   gTGGAGACGCCAATCTTCTGTAATTCCGCGCCATGCGCCAGTCCTTGCGAGATGCCTATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTAAGGGAAGTTCCGCTGGTTGTCAGTatt                                                                                                                                                 <  2:1238592/150‑2 (MQ=255)
    tGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTACCGCTGGTGGTCAGGatat                                                                                                                                                >  2:178801/1‑150 (MQ=255)
      gagaCGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATcc                                                                                                                                              >  2:1305850/1‑150 (MQ=255)
      gagaCGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATcc                                                                                                                                              <  1:490561/150‑1 (MQ=255)
       agaCGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGAAAGTTCCGCTGGTGGTCAGGATATCCt                                                                                                                                             >  1:918573/1‑150 (MQ=255)
            ccAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCAcc                                                                                                                                        >  1:1095498/1‑150 (MQ=255)
            ccAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGGTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCAcc                                                                                                                                        >  1:883473/1‑150 (MQ=255)
                     ggTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAGTGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAg                                                                                                                               >  1:898089/1‑150 (MQ=255)
                          ttCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCg                                                                                                                          <  2:634827/150‑1 (MQ=255)
                            ccGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGtt                                                                                                                        >  2:525109/1‑150 (MQ=255)
                            ccGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGGCAGGATATCCTGCACCATCGGTAAGATTCGtt                                                                                                                        <  2:1073086/150‑1 (MQ=255)
                                 ccAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGGGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTa                                                                                                                   <  2:557981/150‑1 (MQ=255)
                                    aGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGt                                                                                                                <  2:898089/150‑1 (MQ=255)
                                    aGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGt                                                                                                                <  2:1175565/150‑1 (MQ=255)
                                    aGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGt                                                                                                                <  1:730657/150‑1 (MQ=255)
                                    aGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGGCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGt                                                                                                                >  2:833088/1‑150 (MQ=255)
                                      cgcCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTc                                                                                                              <  1:311079/150‑1 (MQ=255)
                                           gCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAg                                                                                                         <  2:1072788/150‑1 (MQ=255)
                                            cccTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAgg                                                                                                        <  1:917815/150‑1 (MQ=255)
                                            cccTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAgg                                                                                                        >  1:442552/1‑150 (MQ=255)
                                            cccTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAgg                                                                                                        <  1:491737/150‑1 (MQ=255)
                                              cTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCCCGAGGGc                                                                                                      >  2:332857/1‑150 (MQ=255)
                                               ttGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCa                                                                                                     <  2:1179830/150‑1 (MQ=255)
                                                 gcgcGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGGAAGATTCGTTGAGTATGTTCACGAGGGCaaa                                                                                                   >  1:566151/1‑150 (MQ=255)
                                                   gcgATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGATCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGc                                                                                                 <  2:197467/150‑1 (MQ=255)
                                                   gcgATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGc                                                                                                 <  1:497133/150‑1 (MQ=255)
                                                   gcgATGCCAATACCAATGCGCACGCAAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATTTCAGTCAGGGAAGTTCCGCTGGTGGTAAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGc                                                                                                 <  2:460042/150‑1 (MQ=255)
                                                        gCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCaaa                                                                                            >  1:252862/1‑150 (MQ=255)
                                                         ccAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCaaaa                                                                                           >  1:1305706/1‑150 (MQ=255)
                                                         ccAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCaaaa                                                                                           <  2:1231809/150‑1 (MQ=255)
                                                          cAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAAt                                                                                          <  1:168033/150‑1 (MQ=255)
                                                           aaTACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATg                                                                                         <  1:667318/150‑1 (MQ=255)
                                                             tACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAg                                                                                       >  1:881152/1‑150 (MQ=255)
                                                               ccAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGTTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCg                                                                                     >  2:647096/1‑150 (MQ=255)
                                                                  aTGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTg                                                                                  <  2:44463/150‑1 (MQ=255)
                                                                     cgcACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACa                                                                               >  1:1226695/1‑150 (MQ=255)
                                                                      gcACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTAAATATCAGTCAGGGAAGTTCCGCTGGTGGGCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAg                                                                              <  1:1341354/150‑1 (MQ=255)
                                                                         cGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTACGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTAc                                                                           <  1:799602/150‑1 (MQ=255)
                                                                           ccAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGCGCGTGAATATCAGTCAGGGAAGTTCCGCTGGTGGTCCGGATATCCTGCACCATCGGGAAGATTCGTTGAGTATGTTGACGAGGGCAAAGCTCATAATGAGCGTTGACCGTACcg                                                                         >  2:550244/1‑150 (MQ=255)
                                                                              gTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACcgtcg                                                                      >  2:1219827/1‑150 (MQ=255)
                                                                                aaaGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGtt                                                                    >  1:598625/1‑150 (MQ=255)
                                                                                 aaGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTc                                                                   <  1:525109/150‑1 (MQ=255)
                                                                                   gCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGGTGAGTATGTACACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTccc                                                                 >  1:393921/1‑149 (MQ=255)
                                                                                    cTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCGGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCcac                                                                >  2:789202/1‑150 (MQ=255)
                                                                                    cTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCGGGATATCCTGCACCATCGGTAAGATTCGGTGAGTATGTTCACGAGGGCAAAGCTCAAAATCAGCGTTGACAGTACCGTCGTTCcac                                                                >  2:789204/1‑150 (MQ=255)
                                                                                            gacgcgCCCGTAAGCCAGAGCGTTACTATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGTGTTGACAGTACCGTCGTTCCACAGGGCGAc                                                        <  2:306594/148‑1 (MQ=255)
                                                                                             ccGCGCCCGTAAGCCATATCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGATGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACa                                                       <  2:1267521/150‑1 (MQ=255)
                                                                                                cgcCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACagag                                                    <  2:370304/150‑1 (MQ=255)
                                                                                                 gcCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCTACAGAGc                                                   <  2:275183/150‑1 (MQ=255)
                                                                                                   ccGGAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGGAAGGATATCCTGCACCAGCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAg                                                 <  2:59449/150‑1 (MQ=255)
                                                                                                           cAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGt                                         >  1:154532/1‑150 (MQ=255)
                                                                                                           cAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGt                                         <  2:234502/150‑1 (MQ=255)
                                                                                                           cAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACAATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGt                                         >  1:862825/1‑150 (MQ=255)
                                                                                                               gCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTcgcg                                     >  1:320957/1‑150 (MQ=255)
                                                                                                               gCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTcgcg                                     >  1:1177176/1‑150 (MQ=255)
                                                                                                                cGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCgg                                    <  2:1107552/150‑1 (MQ=255)
                                                                                                                   tAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTAt                                 >  2:780330/1‑150 (MQ=255)
                                                                                                                     atatCAGTCAGGGAAGTTCCGCTGGTGGTCCGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCCCAGGGCGACAGAGCAGGCCTCTGGCGCGGGATCg                               >  1:87859/1‑150 (MQ=255)
                                                                                                                     atatCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCg                               <  2:816800/150‑1 (MQ=255)
                                                                                                                          agtcagGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGc                          >  1:311985/1‑150 (MQ=255)
                                                                                                                          agtcagGGAAGTTCCGCTGGTGGTCAGGATATCCTGCAACATCGGTAAGATTCGTTGAGTATGTTCAAGAGCGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGCGCGACAGAGCAGGCCTCTGTCGGGGTATCGAAAGc                          >  1:613373/1‑150 (MQ=255)
                                                                                                                           gtctgGTAAGTTCCGCTGGTGGTCGGGATATGCTGCACCATCGGTAAGATTCGTTGCGTATGTGCACGAGGGCAAAGCTCAAACTGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGcc                         <  1:789204/150‑1 (MQ=255)
                                                                                                                           gtcagGGAAGTTCCGCTGGTGGTCGGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGcg                         <  1:789202/150‑2 (MQ=255)
                                                                                                                               gggAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTAACGTCGTTCCACAGGGCGACAGAGCAGGCGTGTGGCGAGGTATCGATAGCCAGaa                     >  1:78857/1‑150 (MQ=255)
                                                                                                                                    gTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCCGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGc                >  1:243733/1‑150 (MQ=255)
                                                                                                                                    gTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGc                >  2:1215266/1‑150 (MQ=255)
                                                                                                                                     ttCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCa               <  2:598625/150‑1 (MQ=255)
                                                                                                                                     ttCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCa               <  1:1300309/150‑1 (MQ=255)
                                                                                                                                     gtCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCa               <  2:571473/149‑1 (MQ=255)
                                                                                                                                      tCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCAt              <  1:57315/150‑1 (MQ=255)
                                                                                                                                       ccGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATg             <  2:1095498/150‑1 (MQ=255)
                                                                                                                                         gCTGGTGGTCCGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATggg           <  2:596235/150‑1 (MQ=255)
                                                                                                                                            ggtggtCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCt        >  2:37580/1‑150 (MQ=255)
                                                                                                                                                  cAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGc  >  1:1054796/1‑150 (MQ=255)
                                                                                                                                                  cAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGc  <  2:1118275/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTCGTTGAGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGC  >  NC_000913/1890984‑1891279

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