Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,495,131 T→C 5.7% V374A (GTG→GCG)  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,1310TC 5.7% 211.1 / 3.0 71V374A (GTG→GCG) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base T (28/38);  new base C (2/2);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

CTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAA  >  NC_000913/3494986‑3495279
                                                                                                                                                 |                                                                                                                                                    
ctgGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACtggtgg                                                                                                                                                    <  4:269730/148‑1 (MQ=255)
ctgGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACtggtgg                                                                                                                                                    <  3:259580/148‑1 (MQ=255)
 tgGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTCCCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGc                                                                                                                                                  >  1:354722/1‑149 (MQ=255)
 tgGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGc                                                                                                                                                  >  5:126761/1‑149 (MQ=255)
   ggCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGTTGCGGTACTGGTGGGc                                                                                                                                                  <  2:112625/147‑1 (MQ=255)
   ggCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGc                                                                                                                                                  >  1:112625/1‑147 (MQ=255)
      gTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATAc                                                                                                                                             <  7:65121/149‑1 (MQ=255)
      gTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGAGGGCGATAc                                                                                                                                             >  4:147261/1‑149 (MQ=255)
       tAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAACCCCTGCTCGGTGCGGTACTGGTGGGCGATAcc                                                                                                                                            >  3:216388/1‑149 (MQ=255)
        aGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATAc                                                                                                                                             >  8:254141/1‑147 (MQ=255)
         gACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAg                                                                                                                                          <  4:224109/149‑1 (MQ=255)
            gTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGa                                                                                                                                       <  4:242282/149‑1 (MQ=255)
              aGGCGGTATTGGTCATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTATACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTTTCAGCGAAGAAAACAAAACCCTGCTCGATGCGGTACTTGTGGGCGATACCAGCGAc                                                                                                                                      >  2:233286/1‑148 (MQ=255)
                gCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTAc                                                                                                                                   <  4:325527/149‑1 (MQ=255)
                 cGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACg                                                                                                                                  <  6:295871/149‑1 (MQ=255)
                 cGGTATTCGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCGACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGACGACAACAAAACCCTCCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTAcc                                                                                                                                  <  6:259826/149‑2 (MQ=255)
                  ggTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACgg                                                                                                                                 <  4:149181/149‑1 (MQ=255)
                  ggTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACgg                                                                                                                                 <  1:50155/149‑1 (MQ=255)
                      ttGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACTCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAAc                                                                                                                             >  4:218820/1‑149 (MQ=255)
                           gATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACctgct                                                                                                                        <  6:128111/149‑1 (MQ=255)
                                gcACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAACCGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAAc                                                                                                                   >  7:372258/1‑149 (MQ=255)
                                gcACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAACACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAAc                                                                                                                   >  8:372199/1‑149 (MQ=255)
                                  aCGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTTGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTg                                                                                                                 <  5:71514/149‑1 (MQ=255)
                                    ggTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGt                                                                                                               <  5:178310/149‑1 (MQ=255)
                                     gTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTTAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTg                                                                                                              <  6:111848/149‑1 (MQ=255)
                                         cACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGa                                                                                                          <  8:253970/149‑1 (MQ=255)
                                          aCGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGaa                                                                                                         <  5:337884/149‑1 (MQ=255)
                                          aCGCCTGGCGCACGTAGCTACGTTAACCTCGACGAAAGTAAAGAGCTCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCACCGACTACGGTAACCTGATCCAACTGGTGCTGaa                                                                                                         >  2:18424/1‑149 (MQ=255)
                                            gCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACg                                                                                                       <  2:306445/149‑1 (MQ=255)
                                            gCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACg                                                                                                       <  6:57193/149‑1 (MQ=255)
                                              cTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGc                                                                                                      >  3:284874/1‑148 (MQ=255)
                                              cTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGc                                                                                                      >  2:59551/1‑148 (MQ=255)
                                                 gcgcACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGTTTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAAGTGGTGCTGAACGCTATc                                                                                                  <  2:149870/149‑1 (MQ=255)
                                                 gcgcACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTAt                                                                                                   <  3:146891/148‑1 (MQ=255)
                                                     aCGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAAc                                                                                              >  3:83526/1‑149 (MQ=255)
                                                     aCGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGACAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGaa                                                                                               >  4:245660/1‑148 (MQ=255)
                                                           cTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCgg                                                                                        <  5:43978/149‑1 (MQ=255)
                                                                    ccTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCgg                                                                               <  1:341704/149‑1 (MQ=255)
                                                                      tCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTACCCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGAt                                                                             >  4:291921/1‑149 (MQ=255)
                                                                      tCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAAc                                                                                                        <  6:62833/122‑1 (MQ=255)
                                                                      tCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAAc                                                                                                        >  5:62833/1‑122 (MQ=255)
                                                                          cgaAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTccc                                                                         <  3:245660/149‑1 (MQ=255)
                                                                            aaaGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGc                                                                                                                                         >  3:63472/1‑83 (MQ=255)
                                                                            aaaGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGc                                                                                                                                         <  4:63472/83‑1 (MQ=255)
                                                                             aaGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGa                                                                      <  2:249350/149‑1 (MQ=255)
                                                                              aGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGAt                                                                     <  4:284874/149‑1 (MQ=255)
                                                                               gTAAAGAGATCTACAAACGCCTGATTTTCAGCGTAGACAACAAAAACCTGCTCTTTGCTGTACTGGTTGTCGATACAATCGACTACGGTAACCTGCTGCAACTTGTTCTGAACGCTATCGAACTGCCGGAATACCCGGATTCCCTGATc                                                                    <  1:233286/149‑1 (MQ=255)
                                                                                 aaaGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGCGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTTAt                                                                     <  7:197766/146‑1 (MQ=255)
                                                                                 aaaGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGCGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGAt                                                                     >  8:197762/1‑146 (MQ=255)
                                                                                 aaaGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGCGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGAt                                                                     <  7:197762/146‑1 (MQ=255)
                                                                                 aaaGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGCGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTAAt                                                                     >  8:197766/1‑146 (MQ=255)
                                                                                       aTCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCg                                                                                         <  7:79472/120‑1 (MQ=255)
                                                                                       aTCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCg                                                                                         >  8:79472/1‑120 (MQ=255)
                                                                                        tCTACAAACGCCTGATTGTCTGCGAAGACAACAAAACACTGCTCGGTGCGGTACTGGTGGGCGATACCATCGACTACGTAAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGc                                                            <  8:336933/148‑1 (MQ=255)
                                                                                          tACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGcaca                                                         >  1:275859/1‑149 (MQ=255)
                                                                                            cAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTTCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGca                                                           <  2:354722/145‑1 (MQ=255)
                                                                                            cAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGcac                                                          <  7:254141/146‑1 (MQ=255)
                                                                                               aCGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGGGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCgg                                                    <  7:345173/149‑1 (MQ=255)
                                                                                                      aTTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGc                                               <  3:218820/147‑1 (MQ=255)
                                                                                                       ttGTCAGCGAAGACAACAAAACCCTGCTCGGTGCTGTACTGGTGTGCGATACCAGCGAGTACGGTAACGTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAATCGCGGATTCCCTGATCCTGCCAGCACACTCGTGTAGCGGc                                            <  1:18424/149‑1 (MQ=255)
                                                                                                        tGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGc                                            <  6:126761/148‑1 (MQ=255)
                                                                                                                aaGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTc                                    >  5:120565/1‑148 (MQ=255)
                                                                                                                aaGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCt                                   <  4:83526/149‑1 (MQ=255)
                                                                                                                 aGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTa                                  >  7:145372/1‑149 (MQ=255)
                                                                                                                 aGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTa                                  >  8:260212/1‑149 (MQ=255)
                                                                                                                 aGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGAGTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTa                                  >  6:250899/1‑149 (MQ=255)
                                                                                                                     aacaaAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACctgc                                                                                                                         <  4:309575/58‑1 (MQ=255)
                                                                                                                     aacaaAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACctgc                                                                                                                         >  3:309575/1‑58 (MQ=255)
                                                                                                                     aacaaAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAAccc                                                                                 >  2:73642/1‑98 (MQ=255)
                                                                                                                     aacaaAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAAccc                                                                                 <  1:73642/98‑1 (MQ=255)
                                                                                                                     aacaaAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCgg                              >  5:177495/1‑149 (MQ=255)
                                                                                                                       caaAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGgtg                            >  4:311303/1‑149 (MQ=255)
                                                                                                                             ccTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGAt                       >  4:29209/1‑148 (MQ=255)
                                                                                                                                                gtggGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAgcgcgc    >  4:259952/1‑148 (MQ=255)
                                                                                                                                                 tggGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCaa  <  7:372199/149‑1 (MQ=255)
                                                                                                                                                 tggGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCaa  <  3:147261/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAA  >  NC_000913/3494986‑3495279

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