Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,260 G→T 7.5% intergenic (+31/‑186) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2600GT 7.5% 185.8 / 2.9 66intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/35);  new base T (3/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAG  >  NC_000913/4060116‑4060404
                                                                                                                                                |                                                                                                                                                
ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                                                              >  2:669945/1‑149 (MQ=255)
ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGCTACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                                                              <  1:517692/149‑1 (MQ=255)
 tCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg                                                                                                                                             <  1:327103/149‑1 (MQ=255)
    tCGATGATGACGAACTGGTAGAAGTGCCTCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                                                          <  1:664879/149‑1 (MQ=255)
    tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAGGGTCACCTGACGGAGAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTAATGAAAGACCTGCCAGCGATGGcc                                                                                                                                          >  2:434953/1‑148 (MQ=255)
       atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTCTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGCAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                                                                       <  2:422580/149‑1 (MQ=255)
        tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAGACCTGCCAGCGATGGCAGGtt                                                                                                                                      <  1:566656/149‑1 (MQ=255)
           tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                                                                   <  1:724912/149‑1 (MQ=255)
           tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTTTTTAGTCATAAAACACCTGCCATCGCTGGCAGGttttt                                                                                                                                   <  2:57412/149‑1 (MQ=255)
            gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                                                                  >  2:650335/1‑149 (MQ=255)
             aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                                                                 <  1:739637/149‑1 (MQ=255)
               gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                                                               >  2:163250/1‑149 (MQ=255)
               gAACCGGTAGAAGCGACTCCGACCTCTATCCGTATTCGTGAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                                                               <  2:409910/149‑1 (MQ=255)
                 aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTGTTTAGTCATAAAAAACCTGCCAGCGATGGCAAGTTTTTTTTTGc                                                                                                                             >  2:657633/1‑149 (MQ=255)
                      tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTAGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTc                                                                                                                        >  1:531277/1‑149 (MQ=255)
                        gAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                                                                      <  1:686635/149‑1 (MQ=255)
                         aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAACCGATAGTCGCCGCGCCAACCGCGCACCGAAAGACGATGAATTTCGTTTTGTAGTGATAAAATACCTGCCAGCGATGGCAGGTTTTATTTTTCTTGTCGTg                                                                                                                     >  1:704986/1‑149 (MQ=255)
                           gggccTCCGACTTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                                                                   <  2:102735/145‑1 (MQ=255)
                              ccTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcgc                                                                                                                <  2:787121/148‑1 (MQ=255)
                                tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcgcca                                                                                                              >  1:384638/1‑147 (MQ=255)
                                  cGACCTCTATCCGTATTCGTAAACGTCACCTGACTGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCacac                                                                                                            <  2:353316/149‑1 (MQ=255)
                                       tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                                                                       >  2:172022/1‑149 (MQ=255)
                                         tATCCTTATTCGTAAACGTCACCTGACGGAAAATGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTATTAGTTATAAAAACCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                    <  1:487397/149‑1 (MQ=255)
                                         tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTACGTTTTTTGGTCATAAAAAACCTGCCGGCGATGGCAGGGTTTGTTTTGCCTTTAGAGAGCGAacacgcgagcc                                                                                                     >  2:621836/1‑142 (MQ=255)
                                         tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGAAGAGTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGGGCGAGCACCGATGCg                                                                                                     >  2:616141/1‑149 (MQ=255)
                                          agtcGTAGTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                    <  2:719797/146‑1 (MQ=255)
                                            ccGTATTCGTAAACGTCACCTGACGGAAAGCGATCGTCGCCGCGCAAACCGCGCACCGAATGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                                                                  <  2:114826/149‑1 (MQ=255)
                                            ccGTATTCGTAAACGTCACCGGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                                                                  <  2:245507/149‑1 (MQ=255)
                                               tATTCGTAAACGTCACCTGACGGAAAACGGTCGTCGCCGCGCCAACCGAGCACCGAAAGACGATTAATTTCGTTGTGTAGTCATAAAAAACCTGACAGCGATGGCAGGTTTTTTTTTGCCTGTAGGGGGCGCACACCAATGAGAGTAta                                                                                               >  2:619535/1‑148 (MQ=255)
                                                       caCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                                                                       <  1:114471/148‑1 (MQ=255)
                                                       aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGGGCGCACACCAATGCGAGTATCTCCCTtat                                                                                       >  1:200816/1‑149 (MQ=255)
                                                                 gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                             <  1:650335/149‑1 (MQ=255)
                                                                 gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                             <  2:211005/149‑1 (MQ=255)
                                                                  aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                                                            >  1:343342/1‑149 (MQ=255)
                                                                      aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGCGTATCACCCTTATAAGCCAGTGGTTACg                                                                        >  2:495144/1‑149 (MQ=255)
                                                                          cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCaa                                                                    >  1:349655/1‑149 (MQ=255)
                                                                             tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                                                                 <  2:641187/149‑1 (MQ=255)
                                                                              cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTGCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                <  2:306358/149‑1 (MQ=255)
                                                                               gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                                                               <  1:163250/149‑1 (MQ=255)
                                                                               gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGGCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                                                               <  1:343533/149‑1 (MQ=255)
                                                                               gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTGTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGGGAGTATCTCCCTTATAAGTCTGTGGTTACGGCACGagag                                                               >  1:359832/1‑149 (MQ=255)
                                                                                  gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                                                            >  1:515896/1‑149 (MQ=255)
                                                                                    cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAGTAGAGAGCtt                                                          <  2:213463/149‑1 (MQ=255)
                                                                                      gtccAACCGCGCACCGAAACCCGATTATTTTCGTTTTTTAGGCACACAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                        <  2:33971/147‑1 (MQ=255)
                                                                                       gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCCCACCAATGCGAGTATCTCCCTTATACGTCTGTGGTTACGTCAATCGAGAGCTTTGt                                                       >  1:614524/1‑149 (MQ=255)
                                                                                         cAACCGTGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGGCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGCCAATAGAGAGCTTTGTca                                                     <  1:630740/149‑1 (MQ=255)
                                                                                         cAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca                                                     <  1:267219/149‑1 (MQ=255)
                                                                                          aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac                                                    <  2:334791/149‑1 (MQ=255)
                                                                                            ccGCGCACCGAAAGACGATTAATTACGTTTTTTCGTCATAACAAACCTGCCAGTGATGGCAGGTATTTTTTTGCGTGTCGTGAACGCACACCAATGCGAGTATCTCCCTTATAAGTGTGAGGTTGCGTAAACGGAGGGCCTCGTcaaca                                                  >  2:145252/1‑146 (MQ=255)
                                                                                                  aCCGAAAGACGATTAATTTCGTTTTTTAGTCATTAAAAACCTGCCATCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCCCCAACCTATGCTATTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                            <  2:330633/149‑1 (MQ=255)
                                                                                                      aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCGGTGGTTACGTCAATAGAAAGCTTTGTCACAATTGTCTGCaa                                        >  2:259107/1‑149 (MQ=255)
                                                                                                           cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                   <  1:669945/149‑1 (MQ=255)
                                                                                                                        tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAtt                      >  1:307304/1‑149 (MQ=255)
                                                                                                                        tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAtt                      >  1:5733/1‑149 (MQ=255)
                                                                                                                            tAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctt                  <  2:214610/149‑1 (MQ=255)
                                                                                                                              gTCATAAAAGACCTGCCGGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttc                <  2:375585/149‑1 (MQ=255)
                                                                                                                              gTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttc                <  1:172022/149‑1 (MQ=255)
                                                                                                                                 aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCAGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGACTGTGGGGACGTCAATAGAGAGCTCTGTCCCAGTTGTCTGCAAAGCCATACGCCGTTAATTGcgtccttt             >  2:617383/3‑149 (MQ=255)
                                                                                                                                  tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttt            <  2:343342/149‑1 (MQ=255)
                                                                                                                                  tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCAGCAAAGTCATACGCCGTTAATTGCTTTCtttt            >  1:545361/1‑149 (MQ=255)
                                                                                                                                  tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGAGAGTATCTCACTTATAAGTCTGTGGTTACGTCAGTAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGACGTTAATTGCTTTCtttt            >  1:151982/1‑149 (MQ=255)
                                                                                                                                    aaaaaCCTGCCATCGCTGGCAGGTCTTTTTTTGCCTGTCGTGAGCGCACACCAATCCGCGTAACTCCCTTATAAGCCTGTGGTTACGTCACTAGCGAGCTTTGCCCCAATTATCTGCCAAGTCATACGCCGTTAATTGCTTTCtctttt          >  2:337773/1‑149 (MQ=255)
                                                                                                                                    aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGGTTTGTCACAATTATCTGCAAAGTCATACGCCGGGAATTGCTTTCtttttt          >  1:268904/1‑149 (MQ=255)
                                                                                                                                     aacaCCTGCCAGGGATGGCCGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGACATACGCCGTTAATTGCTTTCTTTTTTg         <  2:384288/146‑1 (MQ=255)
                                                                                                                                        aCCTGCCATCGCTGGCGGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAATGCTTTCTTTTTTGgcg      >  1:500622/1‑149 (MQ=255)
                                                                                                                                            gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaag  <  1:415708/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAG  >  NC_000913/4060116‑4060404

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