Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,495,641 A→G 5.7% E544G (GAA→GGA)  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,6410AG 5.7% 214.5 / 3.0 70E544G (GAA→GGA) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base G (2/2);  total (36/34)
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.95e-01

AACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTC  >  NC_000913/3495494‑3495789
                                                                                                                                                   |                                                                                                                                                    
aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGaa                                                                                                                                                     <  6:85753/149‑1 (MQ=255)
aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGAGAAGTTTGTAAAACACCCGTCGGTTCGCAGCTGGCCTCCTGCTGGAACGAAAACATTCTGAAGCCTGaa                                                                                                                                                     >  6:281490/1‑149 (MQ=255)
 aCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAc                                                                                                                                                    <  1:7015/149‑1 (MQ=255)
    gTTACATTTGCTCAGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGAGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGGTCGATGCTGGACTCCTGCTGGAACGAATCCATTCTGACGCCGGAACATa                                                                                                                                                 >  8:389126/1‑149 (MQ=255)
        cATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTcc                                                                                                                                             >  1:117338/1‑149 (MQ=255)
         aTTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCg                                                                                                                                            >  1:122608/1‑149 (MQ=255)
         aTTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCg                                                                                                                                            >  5:153592/1‑149 (MQ=255)
         aTTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCg                                                                                                                                            >  1:148167/1‑149 (MQ=255)
           ttGATCCGCGTTGAAGGCATTAAAACCTTCAAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGc                                                                                                                                           <  8:20872/148‑1 (MQ=255)
             gATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTg                                                                                                                                         >  1:25626/1‑148 (MQ=255)
             gATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCAGAACATACTCCGCTg                                                                                                                                         >  5:315727/1‑148 (MQ=255)
              aTCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAACGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCa                                                                                                                                       <  5:281490/149‑1 (MQ=255)
              aTCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGGTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTg                                                                                                                                         >  5:236656/1‑147 (MQ=255)
                  gcgTTCAAGGCATTAAAACCTTAGAAGAACTGCTGGCGAAAGACGTCTAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCTGTTCGATGCTGGCATCCTGCTTGAACGAATACATTCTGAAGCCTGAACATACTCAGCTGCATGAt                                                                                                                                   <  8:135798/149‑1 (MQ=255)
                        aaGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTaa                                                                                                                              <  2:277415/148‑1 (MQ=255)
                            cATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGAc                                                                                                                          >  3:4421/1‑148 (MQ=255)
                             aTTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACaa                                                                                                                        >  1:97448/1‑149 (MQ=255)
                             aTTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACaa                                                                                                                        >  2:91268/1‑149 (MQ=255)
                              ttAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAAc                                                                                                                       <  2:92179/149‑1 (MQ=255)
                              ttAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCAGAACATACTCCGCTGCAGGATTCTAACGACAAc                                                                                                                       <  6:315727/149‑1 (MQ=255)
                                  aaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGAt                                                                                                                                   >  7:183882/1‑133 (MQ=255)
                                  aaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGAt                                                                                                                                   <  8:183882/133‑1 (MQ=255)
                                  aaCCTGCGAAGAACTGCTGGCGAAACACGCCAAAGGCTACGGTTGTGATGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGATAGAATACATTCTGATGCCGGAACATTCTCCGCTGCAGGATTCTAACGACAACTTcc                                                                                                                   <  2:263999/149‑1 (MQ=255)
                                     cTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCg                                                                                                                >  7:368140/1‑149 (MQ=255)
                                      ttCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATa                                                                                                                                                 <  6:304646/115‑1 (MQ=255)
                                      ttCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATa                                                                                                                                                 >  5:304646/1‑115 (MQ=255)
                                        cGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGc                                                                                                               <  2:97448/147‑1 (MQ=255)
                                        cGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCt                                                                                                              <  1:91268/148‑1 (MQ=255)
                                         gaagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTaa                                                                                                            >  8:228689/1‑149 (MQ=255)
                                         gaagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTaa                                                                                                            <  2:122608/149‑1 (MQ=255)
                                         gaagaaCTGCTCTCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCaaa                                                                                                            >  4:311013/1‑149 (MQ=255)
                                                     gCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGaaa                                                                                                  <  2:25626/147‑1 (MQ=255)
                                                      cGAAACACGGCTAAGGCTAGGGGGGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGAc                                                                                               <  4:4421/149‑1 (MQ=255)
                                                       gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGGACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACg                                                                                              >  8:344977/1‑149 (MQ=255)
                                                       gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGGACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACg                                                                                              >  8:339591/1‑149 (MQ=255)
                                                       gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACg                                                                                              >  3:69443/1‑149 (MQ=255)
                                                                    ggCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGt                                                                                 <  2:117338/149‑1 (MQ=255)
                                                                       tACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGAt                                                                              >  5:23273/1‑149 (MQ=255)
                                                                          ggTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATccc                                                                           <  6:236656/149‑1 (MQ=255)
                                                                              gtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGt                                                                       >  5:90744/1‑149 (MQ=255)
                                                                                   gTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                                                            >  6:100572/1‑123 (MQ=255)
                                                                                          aaccaaccGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCATAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGTGCGGt                                                           <  7:228689/149‑1 (MQ=255)
                                                                                          aaccaaccGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGt                                                           >  4:242879/1‑149 (MQ=255)
                                                                                           accaaccGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGt                                                           >  8:243627/1‑148 (MQ=255)
                                                                                              cgacctCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGaaa                                                       <  3:88886/144‑1 (MQ=255)
                                                                                                ccGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAAt                                                      >  3:225427/1‑148 (MQ=255)
                                                                                                  gTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATAAATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAc                                                   >  4:289987/1‑149 (MQ=255)
                                                                                                   tCGGGTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcc                                                  <  8:146107/149‑1 (MQ=255)
                                                                                                     ggTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGaaa                                                       >  4:88886/1‑142 (MQ=255)
                                                                                                     ggTTCGCTGCTGGCCTCCTGCTGGAACGAATAAATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAccc                                                 <  3:289987/148‑1 (MQ=255)
                                                                                                          gctgctGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGaa                                            <  8:368140/148‑1 (MQ=255)
                                                                                                            tgctgGCCTCCTGCTGGAACGAATACATTATGAAGCCGGAACATACTCAGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCCGTGATCCCGCCTTCTCCGGGCGGTGAAATCACCCCGGAAgcg                                         >  1:65147/1‑147 (MQ=255)
                                                                                                              ctgGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGc                                        <  2:148167/148‑1 (MQ=255)
                                                                                                               tgGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCt                                       >  2:359437/1‑148 (MQ=255)
                                                                                                                 gCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGAt                                    <  3:311013/149‑1 (MQ=255)
                                                                                                                 gCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGAAAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGTTGATCCAGCGTTCTCCGCGCGGTTAAATCACCCCGGAAGGGCTGAt                                    >  6:369220/1‑149 (MQ=255)
                                                                                                                   ctcctGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATgg                                  <  7:243627/149‑1 (MQ=255)
                                                                                                                    tcctGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGc                                 >  5:356236/1‑149 (MQ=255)
                                                                                                                        gctgGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCgg                               >  7:217804/1‑147 (MQ=255)
                                                                                                                            gAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTc                         <  6:90744/149‑1 (MQ=255)
                                                                                                                               cGAATACATTCTGAAGCCGGGACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTa                      <  7:339591/149‑1 (MQ=255)
                                                                                                                               cGAATACATTCTGAAGCCGGGACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTa                      <  7:344977/149‑1 (MQ=255)
                                                                                                                                gAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGGAGGTCGTAt                     >  8:363021/1‑149 (MQ=255)
                                                                                                                                           gAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCg                                                             <  3:4873/98‑1 (MQ=255)
                                                                                                                                           gAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCg                                                             >  4:4873/1‑98 (MQ=255)
                                                                                                                                            aaGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGa                                                                               >  8:201162/1‑79 (MQ=255)
                                                                                                                                            aaGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGa                                                                               <  7:201162/79‑1 (MQ=255)
                                                                                                                                            aaGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGaa           <  3:242879/147‑1 (MQ=255)
                                                                                                                                            aaGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATCGCGGTAGGTCGTATCGCGCGTGAAtt         >  8:350849/1‑149 (MQ=255)
                                                                                                                                               ccGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTaa      <  4:69443/149‑1 (MQ=255)
                                                                                                                                                  gAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAAtct   <  6:23273/149‑1 (MQ=255)
                                                                                                                                                  gAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAAtct   <  4:225427/149‑1 (MQ=255)
                                                                                                                                                   aaCATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAAtctc  <  6:356236/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGATGGCGGTAGGTCGTATCGCGCGTGAATTTAATCTC  >  NC_000913/3495494‑3495789

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