Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 703,478 G→A 6.5% P45L (CCG→CTG)  nagB ← glucosamine‑6‑phosphate deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913703,4780GA 6.5% 189.0 / 2.9 61P45L (CCG→CTG) nagBglucosamine‑6‑phosphate deaminase
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base A (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

GAAGAAATTACGGTGCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAAT  >  NC_000913/703330‑703624
                                                                                                                                                    |                                                                                                                                                  
gaagaaATTACGGTGCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCg                                                                                                                                                    >  2:4039/1‑149 (MQ=255)
 aagaaATTACGGTGCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCgg                                                                                                                                                   >  1:209914/1‑149 (MQ=255)
 aagaaATTACGGTGCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCgg                                                                                                                                                   >  6:291791/1‑149 (MQ=255)
   gaaATTACTGTGCATAAAGCGGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGGCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCg                                                                                                                                                 <  5:35960/149‑1 (MQ=255)
    aaaTTACGGTGCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGt                                                                                                                                                >  7:106208/1‑149 (MQ=255)
       ttACGGTGCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTgcc                                                                                                                                             <  8:106208/149‑1 (MQ=255)
             tGCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGt                                                                                                                                       >  4:230763/1‑149 (MQ=255)
              gCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTc                                                                                                                                      >  5:11377/1‑149 (MQ=255)
                      cTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGcc                                                                                                                              <  2:209914/149‑1 (MQ=255)
                        gtagtaGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCg                                                                                                                                                 <  4:192561/128‑1 (MQ=255)
                        gtagtaGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCg                                                                                                                                                 >  3:192561/1‑128 (MQ=255)
                         tagtagCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCa                                                                                                                            <  5:291791/148‑1 (MQ=255)
                         tagtagCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAg                                                                                                                           >  1:100649/1‑149 (MQ=255)
                           gtagCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTa                                                                                                                         >  7:209647/1‑149 (MQ=255)
                               cTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACaaa                                                                                                                     <  8:218111/149‑1 (MQ=255)
                                       aTGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGAt                                                                                                             <  7:186998/149‑1 (MQ=255)
                                             tttCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCgg                                                                                                          <  2:184756/146‑1 (MQ=255)
                                             tttCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAg                                                                                                       >  1:221427/1‑149 (MQ=255)
                                             tttCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCAGTTCGACTTACGCTTTATAGGTGGTCCTCGGCGTGCAGCCAGTCGGCAGGCCAAGTACAAACGGACGAACGGCAg                                                                                                       >  1:232397/1‑149 (MQ=255)
                                                  gCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGt                                                                                                  <  8:294879/149‑1 (MQ=255)
                                                  gCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGt                                                                                                  <  3:230763/149‑1 (MQ=255)
                                                          aCATATTCGTCCATGTTGAAGTTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCAGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACg                                                                                           >  6:213047/1‑148 (MQ=255)
                                                          aCATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCAGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACg                                                                                           >  6:208422/1‑148 (MQ=255)
                                                              atTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCAtt                                                                                       >  2:107949/1‑148 (MQ=255)
                                                               ttCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACGTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGa                                                                                     <  5:135771/149‑1 (MQ=255)
                                                               ttCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGa                                                                                     <  7:168683/149‑1 (MQ=255)
                                                                     cATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGAt                                                                               >  2:35230/1‑149 (MQ=255)
                                                                     cATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGAt                                                                               >  5:40776/1‑149 (MQ=255)
                                                                               gTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATg                                                                     <  2:221427/149‑1 (MQ=255)
                                                                                 gACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGc                                                                   >  3:203985/1‑149 (MQ=255)
                                                                                    aaCGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGatat                                                                      >  4:63698/1‑143 (MQ=255)
                                                                                    aaCGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGatat                                                                      <  3:63698/143‑1 (MQ=255)
                                                                                     aCGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGagc                                                               <  3:120018/149‑1 (MQ=255)
                                                                                       gTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGagcag                                                             >  4:261583/1‑149 (MQ=255)
                                                                                       gTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGa                                                                 >  2:167482/1‑145 (MQ=255)
                                                                                       gTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGa                                                                 <  1:167482/145‑1 (MQ=255)
                                                                                          cTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGcc                                                           >  3:287128/1‑148 (MQ=255)
                                                                                             aaaGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCAtt                                                       >  7:153301/1‑149 (MQ=255)
                                                                                              aaGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCAtt                                                       <  2:130919/148‑1 (MQ=255)
                                                                                                 cTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGcc                                                   >  3:34650/1‑149 (MQ=255)
                                                                                                         gCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCAGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGtt                                           <  5:213047/149‑1 (MQ=255)
                                                                                                         gCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCAGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGtt                                           <  5:208422/149‑1 (MQ=255)
                                                                                                              cTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGc                                       <  3:261583/148‑1 (MQ=255)
                                                                                                              cTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCg                                      >  6:125200/1‑149 (MQ=255)
                                                                                                                     cATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGt                                 >  3:124534/1‑147 (MQ=255)
                                                                                                                      aTTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAg                              >  3:311688/1‑149 (MQ=255)
                                                                                                                      aTTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAg                              >  5:66442/1‑149 (MQ=255)
                                                                                                                       tttCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAgg                             <  2:100649/149‑1 (MQ=255)
                                                                                                                       tttCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAgg                             <  8:209647/149‑1 (MQ=255)
                                                                                                                       tttCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATAGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGTTAGTCAgg                             >  7:126398/1‑149 (MQ=255)
                                                                                                                           gACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATAGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCTACATGTTCAGCGGGAGTCAGGGGGa                         >  5:213545/1‑149 (MQ=255)
                                                                                                                             cTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATc                       <  6:11377/149‑1 (MQ=255)
                                                                                                                                 cGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGtc                   <  5:125200/149‑1 (MQ=255)
                                                                                                                                    tttATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCa                >  2:339161/1‑149 (MQ=255)
                                                                                                                                          ggtggtCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCattatt           <  6:40776/148‑1 (MQ=255)
                                                                                                                                           gtggtCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCa         >  7:36304/1‑149 (MQ=255)
                                                                                                                                            tggtCATCGGCGTGCCGCCAGTTGGCAGGCCCAGTACAAACGGACGATCGGCTGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGt                    >  8:225389/1‑137 (MQ=255)
                                                                                                                                            tggtCATCGGCGTGCCGCCAGTTGGCAGGCCCAGTACAAACGGACGATCGGCTGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGt                    <  7:225389/137‑1 (MQ=255)
                                                                                                                                            tggtCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCAc        <  4:203985/149‑1 (MQ=255)
                                                                                                                                             ggtCATCGGCGTGCCGCCAGTCGGCAGTCCCAGGACAAAAGGACGATCGGCAGTCGGTGTGCACGCATTGAGACGATTGACGATCTGGCGAGCAGCCGATTTGCCGACAGGTGCAGCGGTCGTCCGGGGGATCAGGCGCATGATTCAcc       >  7:227989/1‑149 (MQ=255)
                                                                                                                                              gtCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTAGAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCt      <  4:124534/149‑1 (MQ=255)
                                                                                                                                              gtCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCt      >  1:314465/1‑149 (MQ=255)
                                                                                                                                                  tCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAAt  <  5:68875/149‑1 (MQ=255)
                                                                                                                                                  tCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGAGCGGCAGTCGGTTTGAACGGATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGCCCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCaa   >  5:326245/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAAGAAATTACGGTGCATAAAGCTGTAGTAGCTTTCCGGATGCTCTTTCGGCAGACCGACATATTCGTCCATGTTGAAGGTGACAACGTGCTTAAAGCTGACCTGGCCTGCTTTATGCATTTCGACTAACGCTTTATAGGTGGTCATCGGCGTGCCGCCAGTCGGCAGGCCCAGTACAAACGGACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAAT  >  NC_000913/703330‑703624

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