Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,356,924 C→T P96P (CCG→CCA ycjJ ← putative amino acid/amine transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,9240CT99.1% 373.9 / ‑5.3 110P96P (CCG→CCAycjJputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (67/42);  total (68/42)
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.34e-01

Only 100 of 115 total aligned reads displayed.

GAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAG  >  NC_000913/1356835‑1357013
                                                                                         |                                                                                         
gagagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTggg                                                                                        <  2:2102158/93‑1 (MQ=255)
 agagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGt                                                                                       <  1:2545216/93‑1 (MQ=255)
 agagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGt                                                                                       <  1:1548934/93‑1 (MQ=255)
  gagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGtt                                                                                      <  1:1590862/93‑1 (MQ=255)
  gagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGGGGGtt                                                                                      >  1:3986855/1‑93 (MQ=255)
   agatagatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTa                                                                                     >  1:1551171/1‑93 (MQ=255)
    gatagatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCTGACCATAAATCCGACGTGTGGGTTat                                                                                    >  1:1119160/1‑92 (MQ=255)
    gatagatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTaa                                                                                    >  2:132524/1‑93 (MQ=255)
    gatagatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTaa                                                                                    >  1:1014388/1‑93 (MQ=255)
    gatagatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTaa                                                                                    >  1:4527474/1‑93 (MQ=255)
     atagatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCCACGTGTGGGTTAAt                                                                                   >  2:4264913/1‑93 (MQ=255)
     atagatTTTCGCAAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAAt                                                                                   <  2:478827/93‑1 (MQ=255)
      tagatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATc                                                                                  >  1:1989149/1‑93 (MQ=255)
       agatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCg                                                                                 >  2:1960589/1‑93 (MQ=255)
        gatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGa                                                                                >  2:496984/1‑93 (MQ=255)
        gatTTACGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGa                                                                                <  2:2071196/93‑1 (MQ=255)
         atTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGAc                                                                               <  1:3550706/93‑1 (MQ=255)
         atTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGAc                                                                               <  2:643057/93‑1 (MQ=255)
               gCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTgc                                                                         >  2:341879/1‑93 (MQ=255)
                  aacaaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGcgcg                                                                      >  2:3383962/1‑93 (MQ=255)
                  aacaaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGcgcg                                                                      >  1:2283097/1‑93 (MQ=255)
                  aacaaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGcgcg                                                                      >  1:114166/1‑93 (MQ=255)
                   acaaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGggcgg                                                                     >  1:1676787/1‑92 (MQ=255)
                    caaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTa                                                                    >  2:3820523/1‑93 (MQ=255)
                    caaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTa                                                                    <  1:197478/93‑1 (MQ=255)
                    caaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTa                                                                    >  1:1439448/1‑93 (MQ=255)
                     aaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAg                                                                   >  1:4284189/1‑93 (MQ=255)
                     aaGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAg                                                                   <  1:3955971/93‑1 (MQ=255)
                      aGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAgg                                                                  >  2:1831022/1‑93 (MQ=255)
                      aGACGTGAATCATCGGCAAAAAGATAGAAGCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAgg                                                                  <  2:178925/93‑1 (MQ=255)
                       gACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGATTTTTGCGCGTAGGt                                                                 >  1:882108/1‑93 (MQ=255)
                       gACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGt                                                                 >  1:2812576/1‑93 (MQ=255)
                       gACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGt                                                                 <  2:3891282/93‑1 (MQ=255)
                       gACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGt                                                                 >  1:1368455/1‑93 (MQ=255)
                         cGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTATTCGACTTTTGCGCGTAGGtat                                                               >  2:396190/1‑93 (MQ=255)
                         cGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGtat                                                               <  2:238263/93‑1 (MQ=255)
                         cGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGtat                                                               >  1:255479/1‑93 (MQ=255)
                         cGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGtat                                                               >  2:1831651/1‑93 (MQ=255)
                          gTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGtata                                                              >  1:4162858/1‑93 (MQ=255)
                          gTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGtata                                                              <  2:3652087/93‑1 (MQ=255)
                          gTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGtata                                                              <  1:3464290/93‑1 (MQ=255)
                           ttAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAg                                                             <  1:3480090/93‑1 (MQ=255)
                            tAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAgg                                                            >  2:2962362/1‑93 (MQ=255)
                            tAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAgg                                                            <  2:529376/93‑1 (MQ=255)
                             aaTCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGc                                                           >  2:4384303/1‑93 (MQ=255)
                               tcatcGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCAACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCg                                                         >  1:1419306/1‑93 (MQ=255)
                                 atcGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGaa                                                       >  2:3583247/1‑93 (MQ=255)
                                  tcGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAAc                                                      <  1:958684/93‑1 (MQ=255)
                                       aaaaaGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGcc                                                 <  1:1094924/93‑1 (MQ=255)
                                       aaaaaGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGcc                                                 <  1:1336113/93‑1 (MQ=255)
                                       aaaaaGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGcc                                                 >  2:3991764/1‑93 (MQ=255)
                                          aaGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTcc                                              >  2:2602069/1‑93 (MQ=255)
                                          aaGAGATAATCCAGCAGCAATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTcc                                              <  2:695304/93‑1 (MQ=255)
                                           agagATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCg                                             >  2:4388028/1‑93 (MQ=255)
                                            gagaTAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCgg                                            >  1:784315/1‑93 (MQ=255)
                                             agaTAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGa                                           >  2:4355718/1‑93 (MQ=255)
                                             agaTAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGa                                           <  2:4608930/93‑1 (MQ=255)
                                               aTAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGaaa                                         >  2:3339959/1‑93 (MQ=255)
                                                tAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAAc                                        <  1:1576987/93‑1 (MQ=255)
                                                 aaTCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACt                                       >  1:3776955/1‑93 (MQ=255)
                                                   tCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCTACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTgg                                     <  2:3306186/93‑1 (MQ=255)
                                                   tCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTgg                                     >  2:3185676/1‑93 (MQ=255)
                                                   tCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTgg                                     >  1:4390249/1‑93 (MQ=255)
                                                   tCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTgg                                     >  2:3185677/1‑93 (MQ=255)
                                                    ccAGCAGCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGc                                    >  1:1298251/1‑93 (MQ=255)
                                                     cagcagCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCg                                   >  2:1249571/1‑93 (MQ=255)
                                                     cagcagCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCg                                   <  1:3626581/93‑1 (MQ=255)
                                                     cagcagCGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCg                                   <  2:2074009/93‑1 (MQ=255)
                                                      agcagcGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGa                                  >  2:101968/1‑93 (MQ=255)
                                                      agcagcGATGACCAGCCGACCATAAATCCGACGAGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGa                                  >  2:162523/1‑93 (MQ=255)
                                                       gcagcGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGac                                 >  1:940367/1‑92 (MQ=255)
                                                       gcagcGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGaa                                 >  2:211285/1‑93 (MQ=255)
                                                       gcagcGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGaa                                 >  2:2892158/1‑93 (MQ=255)
                                                         cgcGATGACCAGCCGCCCATAACTCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAACCTGGCGAAcc                               <  1:4217876/92‑1 (MQ=255)
                                                         agcGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAAcc                               >  1:1982488/1‑93 (MQ=255)
                                                          gcGATGACCATCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCa                              <  2:637971/93‑1 (MQ=255)
                                                          gcGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCa                              >  2:465087/1‑93 (MQ=255)
                                                           cGATGACGAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAg                             <  2:1142065/93‑1 (MQ=255)
                                                           cGATGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAg                             >  2:1612875/1‑93 (MQ=255)
                                                             aTGACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGtt                           <  1:3811151/93‑1 (MQ=255)
                                                               gACCAGCCGACCATAAATCCGACGTGTGGGTTAATCTACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGAGGAACCAGTTTg                         >  2:3940225/1‑93 (MQ=255)
                                                               gACCAGCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTg                         >  1:241788/1‑93 (MQ=255)
                                                                    gCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTa                    <  2:216460/93‑1 (MQ=255)
                                                                    gCCGACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTa                    >  1:720679/1‑93 (MQ=255)
                                                                       gACCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCt                 >  2:3901162/1‑93 (MQ=255)
                                                                       caCCATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCt                 <  2:1288185/92‑1 (MQ=255)
                                                                          cATAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGAt              <  2:2546340/93‑1 (MQ=255)
                                                                           aTAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATa             <  1:2928949/93‑1 (MQ=255)
                                                                           aTAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATa             <  1:630063/93‑1 (MQ=255)
                                                                            tAAATCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAg            >  1:815464/1‑93 (MQ=255)
                                                                              aaTCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCg          <  2:1472691/93‑1 (MQ=255)
                                                                              aaTCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCg          >  2:4033119/1‑93 (MQ=255)
                                                                               aTCCGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCgg         <  2:4710579/93‑1 (MQ=255)
                                                                                  cGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTaa      >  1:1640794/1‑93 (MQ=255)
                                                                                  cGACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTaa      >  1:4736947/1‑93 (MQ=255)
                                                                                   gACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTaaa     >  1:1894229/1‑93 (MQ=255)
                                                                                   gACGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTaaa     >  2:2409224/1‑93 (MQ=255)
                                                                                    aCGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGCAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAAc    >  1:1237229/1‑93 (MQ=255)
                                                                                     cGTGTGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACa   <  1:3522521/93‑1 (MQ=255)
                                                                                      gtgtGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTCAACAg  >  1:2888396/1‑93 (MQ=255)
                                                                                         |                                                                                         
GAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAG  >  NC_000913/1356835‑1357013

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