Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,948,237 G→A 6.3% L132L (CTC→CTT)  nudB ← dihydroneopterin triphosphate pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,948,2370GA 6.3% 206.9 / 3.2 63L132L (CTC→CTT) nudBdihydroneopterin triphosphate pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (2/2);  total (33/30)
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.06e-01

TAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTG  >  NC_000913/1948094‑1948379
                                                                                                                                               |                                                                                                                                              
tAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAgcg                                                                                                                                            >  2:259980/1‑148 (MQ=255)
    aTCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgccgc                                                                                                                                        >  2:96714/1‑148 (MQ=255)
       cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgccgccg                                                                                                                                      <  4:91084/147‑1 (MQ=255)
          cgcAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCcaca                                                                                                                                    <  3:91905/146‑1 (MQ=255)
                    cGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATc                                                                                                                         >  3:177555/1‑147 (MQ=255)
                      aTGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAg                                                                                                                      >  1:286669/1‑148 (MQ=255)
                         tCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCa                                                                                                                   >  2:146384/1‑148 (MQ=255)
                         tCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCa                                                                                                                   <  2:283730/148‑1 (MQ=255)
                         tCTGGTGTTGGACCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCa                                                                                                                   <  2:279233/148‑1 (MQ=255)
                          cTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGcc                                                                                                                    >  7:147547/1‑146 (MQ=255)
                          cTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGcc                                                                                                                    <  8:147547/146‑1 (MQ=255)
                             gtgtTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACt                                                                                                                 <  1:259980/146‑1 (MQ=255)
                                   gCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGcc                                                                                                         <  2:282990/148‑1 (MQ=255)
                                        tttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATg                                                                                                    >  7:243546/1‑148 (MQ=255)
                                         ttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                                                                      <  6:25901/145‑1 (MQ=255)
                                         ttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                                                                      >  5:25901/1‑145 (MQ=255)
                                         ttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGt                                                                                                   >  5:36425/1‑148 (MQ=255)
                                         ttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGt                                                                                                   >  7:226774/1‑148 (MQ=255)
                                         ttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGt                                                                                                   >  8:169241/1‑148 (MQ=255)
                                            cTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCCCTTGTAAGCCAGCTGtt                                                                                                  >  5:130967/1‑146 (MQ=255)
                                              aTGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGt                                                                                              >  7:163844/1‑148 (MQ=255)
                                              aTGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGt                                                                                              >  7:249/1‑148 (MQ=255)
                                               tGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGt                                                                                              >  2:63161/1‑147 (MQ=255)
                                                gACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTg                                                                                             >  7:296826/1‑147 (MQ=255)
                                                     gCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAAGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAAc                                                                                        >  2:63596/1‑147 (MQ=255)
                                                      ccATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGa                                                                                      <  1:159943/148‑1 (MQ=255)
                                                           aaaaTATCTCCAGATAGCCCTGCCTGTTCAAGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGc                                                                                 <  1:63596/148‑1 (MQ=255)
                                                                      aGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                                                                      >  1:41701/1‑116 (MQ=255)
                                                                      aGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                                                                      <  2:41701/116‑1 (MQ=255)
                                                                      aGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGcc                                                                                <  1:284754/138‑1 (MQ=255)
                                                                      aGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGcc                                                                                >  2:284754/1‑138 (MQ=255)
                                                                       gATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGa                                                                     <  6:211154/148‑1 (MQ=255)
                                                                         tAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAg                                                                   <  7:251079/148‑1 (MQ=255)
                                                                                gcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGa                                                            >  6:264879/1‑148 (MQ=255)
                                                                                 cctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAg                                                             <  8:243546/146‑1 (MQ=255)
                                                                                  ctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGa                                                            >  7:173459/1‑146 (MQ=255)
                                                                                      tCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTAAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAAc                                                      >  1:7740/1‑148 (MQ=255)
                                                                                      tCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTAAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAAc                                                      >  1:2666/1‑148 (MQ=255)
                                                                                        aGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCa                                                    <  3:109034/148‑1 (MQ=255)
                                                                                          gCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATg                                                  <  8:226774/148‑1 (MQ=255)
                                                                                          gCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATg                                                  <  1:96714/148‑1 (MQ=255)
                                                                                           cAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGAt                                                                                     >  8:196963/1‑112 (MQ=255)
                                                                                           cAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGAt                                                                                     <  7:196963/112‑1 (MQ=255)
                                                                                               gTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTcc                                             <  8:249/148‑1 (MQ=255)
                                                                                                ttAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTAAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCg                                            <  2:7740/148‑1 (MQ=255)
                                                                                                ttAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTAAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCg                                            <  2:2666/148‑1 (MQ=255)
                                                                                                       cAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTAcg                                     >  3:242509/1‑148 (MQ=255)
                                                                                                         aaCTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGa                                                                     >  2:103133/1‑114 (MQ=255)
                                                                                                         aaCTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGa                                                                     <  1:103133/114‑1 (MQ=255)
                                                                                                           cTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTc                                 <  7:169241/148‑1 (MQ=255)
                                                                                                              ttCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACg                              >  2:75118/1‑148 (MQ=255)
                                                                                                                  aTCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCg                          <  1:146384/148‑1 (MQ=255)
                                                                                                                   tCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCgg                         <  6:129322/148‑1 (MQ=255)
                                                                                                                   tCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCgg                         <  7:63407/148‑1 (MQ=255)
                                                                                                                   tCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCgg                         <  8:163844/148‑1 (MQ=255)
                                                                                                                         gCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGcc                                                                                                         >  6:76287/1‑62 (MQ=255)
                                                                                                                         gCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGcc                                                                                                         <  5:76287/62‑1 (MQ=255)
                                                                                                                          ccGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAg                  <  7:273337/148‑1 (MQ=255)
                                                                                                                            ggTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGc                 >  8:77284/1‑147 (MQ=255)
                                                                                                                                   ccAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACg          >  3:3974/1‑147 (MQ=255)
                                                                                                                                     aGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTa        >  1:141966/1‑147 (MQ=255)
                                                                                                                                        aCTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAAtg    >  1:92793/1‑148 (MQ=255)
                                                                                                                                           tAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAAtgtg  >  6:151237/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTG  >  NC_000913/1948094‑1948379

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