Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 639,657 G→A 6.8% intergenic (+149/‑96) ahpC → / → ahpF alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913639,6570GA 6.8% 192.0 / 3.7 59intergenic (+149/‑96)ahpC/ahpFalkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base A (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

TTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGAC  >  NC_000913/639510‑639805
                                                                                                                                                   |                                                                                                                                                    
tttCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCgg                                                                                                                                                     >  2:550731/1‑149 (MQ=255)
   ccTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAc                                                                                                                                                  <  1:115656/149‑1 (MQ=255)
    cTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAcc                                                                                                                                                 >  4:275697/1‑149 (MQ=255)
        gTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCCGCAGCACCCGaa                                                                                                                                             <  4:282695/149‑2 (MQ=255)
        gTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCCGCAGCACCCGaa                                                                                                                                             <  3:307009/149‑2 (MQ=255)
         tCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                                                                                                                            >  4:106028/1‑149 (MQ=255)
         tCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                                                                                                                            >  1:45745/1‑149 (MQ=255)
           tttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtttt                                                                                                                                          >  2:49420/1‑149 (MQ=255)
               aCGCCATAGCGGCGTTGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTcc                                                                                                                                     >  2:243837/1‑149 (MQ=255)
                cGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTcc                                                                                                                                     <  2:45173/149‑1 (MQ=255)
                               ggCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCa                                                                                                                      >  3:446008/1‑149 (MQ=255)
                                gCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCaa                                                                                                                     <  3:86538/149‑1 (MQ=255)
                                     gCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                                <  1:23717/149‑1 (MQ=255)
                                           tCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAGCTTGCATcatg                                                                                                          <  1:43266/149‑4 (MQ=255)
                                           tCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAg                                                                                                          <  3:275697/149‑1 (MQ=255)
                                           tCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAg                                                                                                          <  3:94963/149‑1 (MQ=255)
                                             aCCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGt                                                                                                        >  4:88645/1‑149 (MQ=255)
                                              ccccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTa                                                                                                       <  4:140469/149‑1 (MQ=255)
                                                        ttacttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAgg                                                                                             >  2:75520/1‑149 (MQ=255)
                                                           cttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCg                                                                                          >  1:127800/1‑149 (MQ=255)
                                                            ttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                                >  2:250443/1‑126 (MQ=255)
                                                            ttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                                <  1:250443/126‑1 (MQ=255)
                                                               tgtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                                                                                                                            <  3:503980/95‑1 (MQ=255)
                                                               tgtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                                                                                                                            >  4:503980/1‑95 (MQ=255)
                                                               tgtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                                >  1:31096/1‑123 (MQ=255)
                                                               tgtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                                <  2:31096/123‑1 (MQ=255)
                                                                gtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGc                                                                                     >  1:181163/1‑149 (MQ=255)
                                                                 tAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCa                                                                                    <  2:127800/149‑1 (MQ=255)
                                                                    gCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatga                                                                                 >  3:187489/1‑149 (MQ=255)
                                                                     cTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgat                                                                                <  4:187489/149‑1 (MQ=255)
                                                                        ccGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgatgat                                                                             >  4:260981/1‑149 (MQ=255)
                                                                          ggggATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGt                                                                           <  3:415690/149‑1 (MQ=255)
                                                                                tCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAg                                                                     <  2:263375/149‑1 (MQ=255)
                                                                                 cacaGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGc                                                                    >  2:97845/1‑149 (MQ=255)
                                                                                   caGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGccc                                                                  >  4:202574/1‑149 (MQ=255)
                                                                                    aGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCa                                                                 <  2:181163/149‑1 (MQ=255)
                                                                                    aGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCa                                                                 <  3:88645/149‑1 (MQ=255)
                                                                                     ggCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAg                                                                >  2:131740/1‑149 (MQ=255)
                                                                                             gccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATaaa                                                        >  2:432370/1‑149 (MQ=255)
                                                                                              ccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAAc                                                       <  2:45745/149‑1 (MQ=255)
                                                                                                    tctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTttc                                                                                                                                         >  3:525360/1‑61 (MQ=255)
                                                                                                    tctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTttc                                                                                                                                         <  4:525360/61‑1 (MQ=255)
                                                                                                          gcgAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGAcaca                                           >  2:145978/1‑149 (MQ=255)
                                                                                                                tACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATg                                     >  3:449549/1‑149 (MQ=255)
                                                                                                                    tCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAata                                       >  3:104913/1‑143 (MQ=255)
                                                                                                                    tCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAata                                       <  4:104913/143‑1 (MQ=255)
                                                                                                                                cTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAgg                     >  4:472936/1‑149 (MQ=255)
                                                                                                                                cTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAgg                     <  1:243837/149‑1 (MQ=255)
                                                                                                                                 tAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaa                            <  1:390497/141‑1 (MQ=255)
                                                                                                                                 tAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaa                            >  2:390497/1‑141 (MQ=255)
                                                                                                                                  aaTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCt                   >  3:141692/1‑149 (MQ=255)
                                                                                                                                   aTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCtt                  >  1:570871/1‑149 (MQ=255)
                                                                                                                                    ttcttcGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTa                 >  3:326962/1‑149 (MQ=255)
                                                                                                                                      aatcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGAAGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTAcc               >  4:332048/3‑149 (MQ=255)
                                                                                                                                        tcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCtt             >  2:390876/1‑149 (MQ=255)
                                                                                                                                          gggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTga           >  1:221606/1‑149 (MQ=255)
                                                                                                                                                tgcGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAAtt     >  1:258866/1‑149 (MQ=255)
                                                                                                                                                  cGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGa   <  3:417157/149‑1 (MQ=255)
                                                                                                                                                  cGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGa   <  1:550731/149‑1 (MQ=255)
                                                                                                                                                  cGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGa   <  1:97845/149‑1 (MQ=255)
                                                                                                                                                   ggCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGAc  <  3:202574/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGAC  >  NC_000913/639510‑639805

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