Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,495,860 G→A 6.0% W617* (TGG→TAG)  nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,495,8600GA 6.0% 221.9 / 2.8 68W617* (TGG→TAG) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (35/33)
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.98e-01

ATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTA  >  NC_000913/3495712‑3496004
                                                                                                                                                    |                                                                                                                                                
aTCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTg                                                                                                                                                  <  3:404985/149‑1 (MQ=255)
        ttCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGc                                                                                                                                          <  1:596757/149‑1 (MQ=255)
             cGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGAtt                                                                                                                                     <  2:478418/149‑1 (MQ=255)
                gCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGaa                                                                                                                                  <  1:308239/149‑1 (MQ=255)
                      aaaTCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATATGGCGTCAGCTGATTGAAGCCGGc                                                                                                                            >  2:388717/1‑149 (MQ=255)
                        aTCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCtt                                                                                                                          >  2:22608/1‑149 (MQ=255)
                            ccccGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGaa                                                                                                                      >  4:169714/1‑149 (MQ=255)
                             cccGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGaaa                                                                                                                     >  3:70739/1‑149 (MQ=255)
                             cccGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGaaa                                                                                                                     >  3:543541/1‑149 (MQ=255)
                              ccGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAAc                                                                                                                    >  4:325537/1‑149 (MQ=255)
                              ccGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAAc                                                                                                                    <  3:140726/149‑1 (MQ=255)
                                  aaGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCTTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGt                                                                                                                <  3:458225/149‑1 (MQ=255)
                                   aGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTc                                                                                                               <  1:49186/149‑1 (MQ=255)
                                   aGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCGTCGAAACCGGTc                                                                                                               <  3:43141/149‑1 (MQ=255)
                                        tGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGcc                                                                                                          >  4:573984/1‑149 (MQ=255)
                                           tGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTAt                                                                                                       <  3:84471/149‑1 (MQ=255)
                                           tGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTAt                                                                                                       <  1:296453/149‑1 (MQ=255)
                                                 tAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGaaa                                                                                                 <  1:336258/149‑1 (MQ=255)
                                                 tAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGaaa                                                                                                 <  1:321841/149‑1 (MQ=255)
                                                              gcgTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATgg                                                                                    <  3:23214/149‑1 (MQ=255)
                                                                gTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTc                                                                                                                                             >  1:149301/1‑90 (MQ=255)
                                                                gTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTc                                                                                                                                             <  2:149301/90‑1 (MQ=255)
                                                                  gAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGaa                                                                                >  4:435149/1‑149 (MQ=255)
                                                                   aaTTTAGTCTCTACACCAAGATCACGGGCTCCCAGCGTCTGGCAATGTGTGGCGCACGGAACGCCGATCTGCCGGAGATCTCGCGTCAGCTGATTGAATCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCTTATGGCGaaa                                                                               >  4:175650/1‑149 (MQ=255)
                                                                    aTTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGaaaa                                                                              <  3:169714/149‑1 (MQ=255)
                                                                     tttAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGaaa                                                                                                 >  4:73429/1‑129 (MQ=255)
                                                                     tttAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGaaa                                                                                                 <  3:73429/129‑1 (MQ=255)
                                                                     tttAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAAc                                                                             >  3:258794/1‑149 (MQ=255)
                                                                     tttAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAAc                                                                             >  3:199215/1‑149 (MQ=255)
                                                                     tttAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTAGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAAc                                                                             >  4:175066/1‑149 (MQ=255)
                                                                     tttAATCGCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTAGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAAc                                                                             >  2:412440/1‑149 (MQ=255)
                                                                           ctctACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTAGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGt                                                                       <  1:412440/149‑1 (MQ=255)
                                                                           ctctACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTAGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGt                                                                       <  3:175066/149‑1 (MQ=255)
                                                                            tctACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGTGAAAACCTGCGTg                                                                      >  3:529409/1‑149 (MQ=255)
                                                                            tctACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTg                                                                      >  1:540730/1‑149 (MQ=255)
                                                                            tctACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTg                                                                      >  4:247749/1‑149 (MQ=255)
                                                                                cacCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTa                                                                  <  2:474966/149‑1 (MQ=255)
                                                                                 acCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAg                                                                 >  1:265159/1‑149 (MQ=255)
                                                                                 acCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAg                                                                 >  3:221017/1‑149 (MQ=255)
                                                                                 acCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAg                                                                 >  3:290247/1‑149 (MQ=255)
                                                                                  ccAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGc                                                                >  2:137151/1‑149 (MQ=255)
                                                                                   cAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCa                                                               >  1:353972/1‑149 (MQ=255)
                                                                                    aaGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGc                                                                                          >  2:447682/1‑121 (MQ=255)
                                                                                    aaGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGc                                                                                          <  1:447682/121‑1 (MQ=255)
                                                                                      gATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCt                                                            >  1:81564/1‑149 (MQ=255)
                                                                                       aTCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTg                                                           <  4:378323/149‑1 (MQ=255)
                                                                                          aCTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGaaa                                                                                                 <  1:591766/108‑1 (MQ=255)
                                                                                          aCTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGaaa                                                                                                 >  2:591766/1‑108 (MQ=255)
                                                                                            tGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGcc                                                      >  3:307723/1‑149 (MQ=255)
                                                                                            tGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGcc                                                      <  1:325337/149‑1 (MQ=255)
                                                                                            tGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGcc                                                      <  1:315727/149‑1 (MQ=255)
                                                                                                tCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGaaa                                                                                                 <  1:459587/102‑1 (MQ=255)
                                                                                                tCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGaaa                                                                                                 >  2:459587/1‑102 (MQ=255)
                                                                                                tCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTa                                                  >  4:533703/1‑149 (MQ=255)
                                                                                                 cccAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTAc                                                 <  3:573984/149‑1 (MQ=255)
                                                                                                       gTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGTAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGtt                                           <  2:540730/149‑1 (MQ=255)
                                                                                                       gTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGtt                                           <  3:325537/149‑1 (MQ=255)
                                                                                                       gTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGtt                                           <  4:290247/149‑1 (MQ=255)
                                                                                                                tGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGc                                  <  4:166748/149‑1 (MQ=255)
                                                                                                                 gTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGc                                        >  2:355119/1‑142 (MQ=255)
                                                                                                                 gTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGc                                        <  1:355119/142‑1 (MQ=255)
                                                                                                                            aGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGc                      >  3:512658/1‑149 (MQ=255)
                                                                                                                                    gATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACt              >  3:136217/1‑149 (MQ=255)
                                                                                                                                    gATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACt              <  1:492509/149‑1 (MQ=255)
                                                                                                                                    gATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACt              >  1:198773/1‑149 (MQ=255)
                                                                                                                                     aTCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTg             >  1:104417/1‑149 (MQ=255)
                                                                                                                                          ccGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGaaaa        <  3:273084/149‑1 (MQ=255)
                                                                                                                                               gaTCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCt   >  4:326847/1‑149 (MQ=255)
                                                                                                                                                aTCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTa  <  3:241241/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTCTACACCAAGATCACTGGCTCCCAGCGTCTGGCGATGTTTGGCGCACAGAAAGACGATCTGCCGGAGATCTGGCGTCAGCTGATTGAAGCCGGCTTCGAAACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTA  >  NC_000913/3495712‑3496004

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