Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,756,561:1 +G 8.0% intergenic (‑51/+114) aldB ← / ← yiaY aldehyde dehydrogenase B/L‑threonine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,756,5611.G 8.0% 300.8 / 22.4 75intergenic (‑51/+114)aldB/yiaYaldehyde dehydrogenase B/L‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (36/33);  new base G (4/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-02

GCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTC  >  NC_000913/3756421‑3756707
                                                                                                                                             |                                                                                                                                                  
gccAAGAAAGTTGTTATAGCGTGCTTTTAACCTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGc                                                                                                                                            <  2:426817/149‑1 (MQ=255)
   aaTAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGCGACTGGCGAAGAT‑TTCGCCAg                                                                                                                                         <  1:230450/149‑1 (MQ=255)
        cgTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGAT‑TTCGCCAGTcacg                                                                                                                                    <  1:446171/148‑4 (MQ=255)
        aGTTTTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑GATTTCGCCAGTcacg                                                                                                                                    <  1:238357/149‑4 (MQ=255)
          ttGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTaa                                                                                                                                  >  2:235156/1‑149 (MQ=255)
             tcacggCGGGCTTTTACCTTGAGGGGGAAACCATACTCCCCGGGCGTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCGCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAg                                                                                                                               <  1:247088/144‑1 (MQ=255)
             tCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAGTCTGTGCGGGAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGGGTGAGGTATAACCAGGGTAGCCGTGACTGGCGAG‑ATCTTCGCCAGCCGGTAACAg                                                                                                                               >  2:43046/1‑149 (MQ=255)
             tCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGACGGGGGCTTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGGATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCACAAATCGGTAACAg                                                                                                                               >  1:54688/1‑149 (MQ=255)
               aTAGCGGGCTTTTAGCTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATAATTGGTCATGATCGTCTCCGTGCGGTGTGAGGTATAACAAGGGTAAACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTagaggct                                                                                                                             >  2:351052/1‑143 (MQ=255)
                tAGCGGGCTTTTAACTTGAGGGTGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTGTACACCGGCTCAACGTGACTGGCGGA‑ATCTTCGCCAGTCGGTAACAGCtt                                                                                                                            >  2:357843/1‑149 (MQ=255)
                tAGCGGGCTTTTAACTTGAGGGGGAAACCATCCTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGCCGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCtt                                                                                                                            <  1:347125/149‑1 (MQ=255)
                      gCTTCTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGAc                                                                                                                      <  2:73512/149‑1 (MQ=255)
                          ttAACTTGAGGGGGCAACCATACTCGCCGGGCTTAAACTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTAGAAAAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCAGGTCGGTAACAGCTTCActaggacc                                                                                                                  >  2:448839/1‑142 (MQ=255)
                               tttggggggAAACCATACTCGCCGGGCTTGATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCAc                                                                                                             <  2:236680/145‑1 (MQ=255)
                               tttaGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCAc                                                                                                             <  2:376516/146‑1 (MQ=255)
                                tGAGGGGGAAACCATACTCGCCGAGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGCTAGTCGTGATGTGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACg                                                                                                            >  1:412594/1‑149 (MQ=255)
                                 gAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTCACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                                                                           <  1:538410/149‑1 (MQ=255)
                                  aGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGaa                                                                                                          <  1:380615/149‑1 (MQ=255)
                                         aaCCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACgg                                                                                                   <  1:299428/149‑1 (MQ=255)
                                                cTCGCCTGGCTTAATCTGGGCTGGAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAggg                                                                                            <  2:690435/149‑1 (MQ=255)
                                                 tCGCCGGGCTTAAGCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGt                                                                                           <  2:597786/149‑1 (MQ=255)
                                                   gCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGGTa                                                                                         >  1:704849/1‑149 (MQ=255)
                                                   gCCGGGCTTAATCTGTGCTGAAGGGGGAATATTGGTTATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCCGTCGGTAAAGGCTTTGCGACAGCTATCACGAATTTACGGGCAAGGGTTa                                                                                         >  1:68187/1‑149 (MQ=255)
                                                    ccGGGCTTAATCTGAGCTGAAGGGGGATTATTGGTCTTGATCATCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTATGGGCAAGGGTTAc                                                                                        <  1:552231/149‑1 (MQ=255)
                                                        gCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACagag                                                                                    <  2:730089/149‑1 (MQ=255)
                                                        gCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACagag                                                                                    <  1:725132/149‑1 (MQ=255)
                                                         cTTAATCTGTGCTGAAGGGGGATTATTTGGCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAgaga                                                                                   <  1:748010/149‑1 (MQ=255)
                                                          ttAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGAt                                                                                  >  1:381366/1‑149 (MQ=255)
                                                                tgtgCTGAAGGGGGATTATTGGTCGTGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCATGGGGTACAGAGATAGTTGa                                                                            >  1:579842/1‑149 (MQ=255)
                                                                tgtgCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGa                                                                            <  1:89721/149‑1 (MQ=255)
                                                                  ttcTGAAGGGGGATTATTGCTCATGACCGTCTCCTTGCGGTGTGAGGTATGACAAGGGTAGACGCGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaa                                                                          <  1:353868/147‑1 (MQ=255)
                                                                  ggctctaGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCCGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaa                                                                          <  2:571301/143‑1 (MQ=255)
                                                                    cTGAAGGGGGATTATTAGTCATGATCGTCTCCTTGCGGTGTGAGTTATAACAAAGGTAGACGTGACTTGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaaaa                                                                        >  2:553456/1‑149 (MQ=255)
                                                                     tGAAGTGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                       <  1:3604/149‑1 (MQ=255)
                                                                     tGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                       <  1:140617/149‑1 (MQ=255)
                                                                     tGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                       <  2:745274/149‑1 (MQ=255)
                                                                         gggggATTATTGTTTATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGTCGTCCCGGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGc                                                                   <  1:768874/149‑1 (MQ=255)
                                                                             gATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAAc                                                               >  1:331150/1‑149 (MQ=255)
                                                                             gATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGAAGACGTGACTGGCGAA‑ATCTTCGCCAGCCGGTAACAGCTTTCCGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTCGAAAACGCGGCTAAc                                                               >  1:317583/1‑149 (MQ=255)
                                                                                  ttGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAAttt                                                          <  1:139593/149‑1 (MQ=255)
                                                                                    ggTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGc                                                        <  1:11864/149‑1 (MQ=255)
                                                                                     gTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAACAGCGGCTAACAGTTTGcc                                                       >  1:488891/1‑149 (MQ=255)
                                                                                          gATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCg                                                  >  2:752815/1‑149 (MQ=255)
                                                                                            tCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACTGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGtt                                                >  1:63496/1‑149 (MQ=255)
                                                                                                    ttGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATAACGAATTTACGGGCCAGGGTTACAGAGATAGTTGAAAAAGCGGCTCACAATTTGCCAGCCGTTGTGGAAat                                        >  1:547984/1‑149 (MQ=255)
                                                                                                    ttGCGGTGTGAGGTATAACAAGGGGAGACGTGACTGGCGAA‑AGCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAGTTTACGGGCTAGGGTTACAGAGATATTTGAAAAAGCGGCTAACATTTTGCCAGCCGTTGTGGAAat                                        <  2:220689/149‑1 (MQ=255)
                                                                                                      gCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAAGTGGTTTCCAACCGTTGTGGaaagta                                      >  2:197317/1‑146 (MQ=255)
                                                                                                      gCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAatga                                      >  1:666851/1‑149 (MQ=255)
                                                                                                       cGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGACAAAGCGGCTCACAATTTGCCCGCCGTTGTGGAAatgat                                     <  2:189138/149‑1 (MQ=255)
                                                                                                         gtgtGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACCGCTGTACGACAGCTGTTACGAATTTACGGGCAAGGGCTACAGAGATAGTTGAAAAAGCGGCTCACAATTTGCCAGCCGTTGTGGAAAtggtga                                   >  2:471897/1‑149 (MQ=255)
                                                                                                         gtgtGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCACGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCCGCCGTTGTGGAAAtgatga                                   >  2:293169/1‑149 (MQ=255)
                                                                                                         gtgtGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCAGTAACAGCTTTACGACAGCGATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTAAAAAAGCGGCTAGCAATTTGCCAGCCGTTGTGGAAAtgatga                                   >  2:38318/1‑149 (MQ=255)
                                                                                                         gtgtGAGGGATAACAAGGGTAGACGTGACTGGCGAA‑AACTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAATAATTAGCCAGCCGTTGTGGAAAtgatga                                   <  1:511547/149‑1 (MQ=255)
                                                                                                            tgAGGTATAACAAAGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTa                                >  1:535959/1‑149 (MQ=255)
                                                                                                             gAGGTATAACAAGGGTAGACGTGACAGGAAAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTAGAGGGATAGTTGAAAAAGCGGCTAACAAGTTGCCAGCGGTTGTGGAAATGATGATTAc                               >  1:255635/1‑149 (MQ=255)
                                                                                                              aGGTATAACAAGGGTAGAGGTGACTGGCGAA‑ATCTTCGCCAGGCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACa                              <  1:218692/149‑1 (MQ=255)
                                                                                                              aGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACa                              >  1:754602/1‑149 (MQ=255)
                                                                                                                gTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAACGAATAACAATTTGCCAGCCGTTGTGGAGATGATGAGTACaac                            >  2:257944/1‑147 (MQ=255)
                                                                                                                    aaCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTg                        >  1:626749/1‑149 (MQ=255)
                                                                                                                    aaCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTg                        >  1:773509/1‑149 (MQ=255)
                                                                                                                     aCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgc                       >  1:292701/1‑149 (MQ=255)
                                                                                                                     aCAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgc                       <  1:250692/149‑1 (MQ=255)
                                                                                                                     aCAAGGGTAGACGTGACTGGCGAA‑AGCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgc                       <  1:247047/149‑1 (MQ=255)
                                                                                                                      cAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGGTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgcg                      >  1:83340/1‑149 (MQ=255)
                                                                                                                      cAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgcg                      >  2:591237/1‑149 (MQ=255)
                                                                                                                         gggTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGa                   <  1:525099/149‑1 (MQ=255)
                                                                                                                         gggTAGACGGGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGa                   <  2:255313/149‑1 (MQ=255)
                                                                                                                            tAGACGTGACTGGCGAG‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATACCGACTTTACGCGCAAGGGGTTCAGAGATAGTTGACAATTCTGTGAACAGGCGGCCATCCGTTTTGGAAATCATGATTACAACGCTCcgcgcgca                >  2:552473/1‑145 (MQ=255)
                                                                                                                                cGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCg            >  1:686936/1‑149 (MQ=255)
                                                                                                                                   gACTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCca         >  2:373633/1‑149 (MQ=255)
                                                                                                                                    aCTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCcac        >  1:714116/1‑149 (MQ=255)
                                                                                                                                    aCTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCcac        >  1:549234/1‑149 (MQ=255)
                                                                                                                                    aCTGGCGAAGAT‑TTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCCCGATCAATCGcccc        >  2:81488/1‑147 (MQ=255)
                                                                                                                                       ggCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAAt     >  2:496035/1‑149 (MQ=255)
                                                                                                                                          gAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATttc  >  1:753267/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAA‑ATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTC  >  NC_000913/3756421‑3756707

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