Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,519 T→G 6.7% intergenic (+44/‑204) ygaM → / → nrdH DUF883 domain‑containing protein YgaM/glutaredoxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,5190TG 6.7% 148.7 / 3.2 60intergenic (+44/‑204)ygaM/nrdHDUF883 domain‑containing protein YgaM/glutaredoxin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base G (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

CTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAAC  >  NC_000913/2800384‑2800661
                                                                                                                                       |                                                                                                                                              
cTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTGTTTGACCCCGCAGCCTGTCTTCggg                                                                                                                                         >  1:203300/1‑143 (MQ=255)
 ttttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGGCTTTGACCCCGAAGCCTGTCTTCgggg                                                                                                                                        >  1:388444/1‑143 (MQ=255)
  tttGTTCGTGCGAGACCCTGGTACAGCGTGCGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                                                                       <  1:620755/143‑1 (MQ=255)
   ttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGCAGCCTGTCTTCGGGGtt                                                                                                                                      >  1:212877/1‑143 (MQ=255)
      ttCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCt                                                                                                                                   >  1:277112/1‑143 (MQ=255)
      tgcGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTAACCCCGAAGACAGGCTTCGGGGtcaaa                                                                                                                                   <  1:200460/141‑5 (MQ=255)
       tCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCCGGGTTTCtt                                                                                                                                  >  1:461432/1‑143 (MQ=255)
         gTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtttt                                                                                                                                >  1:776524/1‑143 (MQ=255)
          tGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATACGCCAATAGTAATGCAAAAATGATACTCCATACGCGCCTTGGCCTCCGAAACGTCTCTTCGGGGTTTCtcttt                                                                                                                               >  1:638973/1‑143 (MQ=255)
               ggcccTTGCGTAGCGTGCGTACAGTAGCTCCGGTAGGTATTCTTCTCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATCATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTg                                                                                                                          <  2:131912/141‑1 (MQ=255)
                aCCCTGGTGTAGCGAGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                                                                         <  1:572087/143‑1 (MQ=255)
                 ccgcggTGGAGCGGGGGGACAGTAGCTGCGGTAGGTATTTGTATCGGGGCGCTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                        <  2:607527/139‑1 (MQ=255)
                 cccTGGTGGAGCGGGGGTACGGGAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                        <  2:602918/143‑1 (MQ=255)
                    tGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGaa                                                                                                                     <  1:612968/143‑1 (MQ=255)
                      gtgtAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGCTGAATc                                                                                                                   >  2:54997/1‑143 (MQ=255)
                        gtAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCCGTCTTCGGGGTTTCTTTTTGACTGGTGGATcac                                                                                                                 >  1:742939/1‑143 (MQ=255)
                           gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTGCTTTTTGCCTGGTGAATCACaaa                                                                                                              >  1:781578/1‑143 (MQ=255)
                               gggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATc                                                                                                          >  2:492699/1‑143 (MQ=255)
                                ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCAGTGTGAAGCATGCGCAAAGAGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCAGAGGCCTGTCTTCGGGGGTTCTTTTTGCCTGGTGAAGCACGAAAATcg                                                                                                         >  2:578313/1‑142 (MQ=255)
                                     agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTa                                                                                                    <  1:164450/143‑1 (MQ=255)
                                        agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAccc                                                                                                 >  2:676007/1‑143 (MQ=255)
                                           tGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGt                                                                                              >  2:370252/1‑143 (MQ=255)
                                              ggtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGAGTCTTTGACCCCGAAGACTGTCTTAGGGGTTTCTCTTTCACTGGTGAATCACAAAAAACCCCCTACCCCATCAc                                                                                           >  2:129647/1‑143 (MQ=255)
                                                  ggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTc                                                                                       <  2:178696/143‑1 (MQ=255)
                                                   gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                                                                      >  2:587577/1‑143 (MQ=255)
                                                   gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGACTGGTGAATCACAAAAATCACCCTACCCCGTCACGCTCa                                                                                      >  2:251462/1‑143 (MQ=255)
                                                   gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTATTTTTGCAAGGTGAATCCCAAAAATCCCCCTACCCCGTCACGCTCa                                                                                      <  1:808351/143‑1 (MQ=255)
                                                       ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTATTTGACAACGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATAACAACAATCACCCTACCCCGTCACGCTCAGATc                                                                                  >  1:368819/1‑143 (MQ=255)
                                                          tttgcgcgcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                                                               <  2:365888/138‑1 (MQ=255)
                                                            tCGGCGCACTGTTAAGCATGCGCCAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAggg                                                                             <  2:28892/143‑1 (MQ=255)
                                                            tCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAACCTGTCATCGGCTATTCTGTTTGCCTGGTGCATCCCCAAAGTCCACCTCCCCGGTCACGCTCCGAACCAggg                                                                             >  2:282230/1‑143 (MQ=255)
                                                               gcgcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTaa                                                                          <  1:281605/143‑1 (MQ=255)
                                                                  cACTGTTAAGCGTGCGATAATCGTAGTGCAAAACTGATCATAAATACGCGTCTTTGACCCCGCAGCCTGTCTCCGGGGTTGCTTTTTGCCTGGTGAATCACAAAAGTCCCCCTACCCCGTCACGCTCATATCCAGGGTAAttt                                                                       <  2:227597/143‑1 (MQ=255)
                                                                  cACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAAttt                                                                       >  1:34893/1‑143 (MQ=255)
                                                                   cgtGTTAAGTATGCGGAAATCGTAGTGCAAAAAGGATAATAAGTACGCGTCTTTGTCCCCGAAGCCTGTCTTCGGTGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTc                                                                      <  2:800193/141‑1 (MQ=255)
                                                                   aCTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCACGTCGTCCGGGTTTCTTCCTGCCTGCTGAATCACAAAAATCCCCCTACCCCGTCCCCCCCACATCCAGGCTAAtctc                                                                      >  2:744793/1‑143 (MQ=255)
                                                                            cATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAGGCCTGTCTTCCGGGGTACTTTTTGCCTGGTGCACACCACATAACCCCCTACCCCGTCACCCGCATATCCAGCGTAATTCCCACAACTag                                                             >  2:130282/1‑142 (MQ=255)
                                                                             aTGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAGTCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAtt                                                            >  1:265329/1‑143 (MQ=255)
                                                                                cgcAAATCGTAGTGAAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc                                                         <  2:203300/143‑1 (MQ=255)
                                                                                   aaaTCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtat                                                      >  1:53150/1‑143 (MQ=255)
                                                                                            tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAACATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgt                                             <  1:730144/143‑1 (MQ=255)
                                                                                            tGCAAAAATGATAATAAATACGCGTCGTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgt                                             <  2:280533/143‑1 (MQ=255)
                                                                                                aaaaTGATTATAAATACGCGTCTTTGTCCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                         <  2:720906/143‑1 (MQ=255)
                                                                                                aaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                         <  1:480585/143‑1 (MQ=255)
                                                                                                aaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                         <  1:523956/143‑1 (MQ=255)
                                                                                                              taCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCAAAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCaa                           <  1:185129/143‑1 (MQ=255)
                                                                                                                 gcgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAGTTTCGACCACTATTTGCTCTATATTGTGTGGTTGAATCTTTTTTCAACTa                        >  1:51526/1‑143 (MQ=255)
                                                                                                                   gTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACa                      <  2:18648/143‑1 (MQ=255)
                                                                                                                    tCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACAt                     <  1:605774/143‑1 (MQ=255)
                                                                                                                      tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATACCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATAGATTGTGTGGTTGAGTCTTTTTCCACCCACATCt                   >  2:105512/1‑143 (MQ=255)
                                                                                                                        tGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                 <  1:611569/143‑1 (MQ=255)
                                                                                                                           ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGAATAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTCCATCTAGTAt              >  1:91477/1‑143 (MQ=255)
                                                                                                                            cccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtc             >  1:418076/1‑143 (MQ=255)
                                                                                                                            cccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGATGCACAGCGAGGGGTTGAATCTTTTTTCAACTACATCTAGTAtc             >  2:354375/1‑143 (MQ=255)
                                                                                                                             ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCGCTAGTTGCTAGATATTGTGTGGTTGAATCTTTTTTCAACTGCATCTAGTAtct            >  1:258716/1‑143 (MQ=255)
                                                                                                                             ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtct            <  2:520372/143‑1 (MQ=255)
                                                                                                                             ccGAAGCCTGTCTTCGGGGTGTCTTTTTGCCTGGTGAAGCACAAAAATCCCCCTACCCCGTACCGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATGGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtct            <  2:604779/143‑1 (MQ=255)
                                                                                                                              cGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtctc           <  2:260293/143‑1 (MQ=255)
                                                                                                                                aaGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTg         <  1:492699/143‑1 (MQ=255)
                                                                                                                                     tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCa    >  2:211445/1‑143 (MQ=255)
                                                                                                                                       tCTTCGGGGTTTCTTTTTTCTTGGTCAATTTCAATATTCCCCCTTCCCCGTCACGCTAAAATCCAGGGTAATTTCGACCACTATTTGCGATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAAc  <  1:602049/143‑1 (MQ=255)
                                                                                                                                       tCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTAACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAAc  >  1:74467/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                              
CTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAAC  >  NC_000913/2800384‑2800661

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