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

GCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTC  >  NC_000913/3756421‑3756707
                                                                                                                                             |                                                                                                                                                  
gccAAGAAAGTTGTTATAGCGTGCTTTTAACCTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGc                                                                                                                                            <  2:426817/149‑1 (MQ=255)
   aaTAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGCGACTGGCGAAGATTTCGCCAg                                                                                                                                         <  1:230450/149‑1 (MQ=255)
        cgTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGATTTCGCCAGTcacg                                                                                                                                    <  1:446171/148‑4 (MQ=255)
        aGTTTTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGATTTCGCCAGTcacg                                                                                                                                    <  1:238357/149‑4 (MQ=255)
          ttGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTaa                                                                                                                                  >  2:235156/1‑149 (MQ=255)
             tcacggCGGGCTTTTACCTTGAGGGGGAAACCATACTCCCCGGGCGTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCGCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAg                                                                                                                               <  1:247088/144‑1 (MQ=255)
             tCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAGTCTGTGCGGGAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGGGTGAGGTATAACCAGGGTAGCCGTGACTGGCGAGATCTTCGCCAGCCGGTAACAg                                                                                                                               >  2:43046/1‑149 (MQ=255)
             tCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGACGGGGGCTTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGGATAACAAGGGTAGACGTGACTGGCGAAATCTTCACAAATCGGTAACAg                                                                                                                               >  1:54688/1‑149 (MQ=255)
               aTAGCGGGCTTTTAGCTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATAATTGGTCATGATCGTCTCCGTGCGGTGTGAGGTATAACAAGGGTAAACGTGACTGGCGAAATCTTCGCCAGTCGGTagaggct                                                                                                                             >  2:351052/1‑143 (MQ=255)
                tAGCGGGCTTTTAACTTGAGGGTGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTGTACACCGGCTCAACGTGACTGGCGGAATCTTCGCCAGTCGGTAACAGCtt                                                                                                                            >  2:357843/1‑149 (MQ=255)
                tAGCGGGCTTTTAACTTGAGGGGGAAACCATCCTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGCCGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCtt                                                                                                                            <  1:347125/149‑1 (MQ=255)
                      gCTTCTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGAc                                                                                                                      <  2:73512/149‑1 (MQ=255)
                          ttAACTTGAGGGGGCAACCATACTCGCCGGGCTTAAACTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTAGAAAAAGGGTAGACGTGACTGGCGAAATCTTCGCAGGTCGGTAACAGCTTCActaggacc                                                                                                                  >  2:448839/1‑142 (MQ=255)
                               tttggggggAAACCATACTCGCCGGGCTTGATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCAc                                                                                                             <  2:236680/145‑1 (MQ=255)
                               tttaGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCAc                                                                                                             <  2:376516/146‑1 (MQ=255)
                                tGAGGGGGAAACCATACTCGCCGAGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGCTAGTCGTGATGTGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACg                                                                                                            >  1:412594/1‑149 (MQ=255)
                                 gAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTCACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGa                                                                                                           <  1:538410/149‑1 (MQ=255)
                                  aGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGaa                                                                                                          <  1:380615/149‑1 (MQ=255)
                                         aaCCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACgg                                                                                                   <  1:299428/149‑1 (MQ=255)
                                                cTCGCCTGGCTTAATCTGGGCTGGAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAggg                                                                                            <  2:690435/149‑1 (MQ=255)
                                                 tCGCCGGGCTTAAGCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGt                                                                                           <  2:597786/149‑1 (MQ=255)
                                                   gCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGGTa                                                                                         >  1:704849/1‑149 (MQ=255)
                                                   gCCGGGCTTAATCTGTGCTGAAGGGGGAATATTGGTTATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCCGTCGGTAAAGGCTTTGCGACAGCTATCACGAATTTACGGGCAAGGGTTa                                                                                         >  1:68187/1‑149 (MQ=255)
                                                    ccGGGCTTAATCTGAGCTGAAGGGGGATTATTGGTCTTGATCATCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTATGGGCAAGGGTTAc                                                                                        <  1:552231/149‑1 (MQ=255)
                                                        gCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACagag                                                                                    <  2:730089/149‑1 (MQ=255)
                                                        gCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACagag                                                                                    <  1:725132/149‑1 (MQ=255)
                                                         cTTAATCTGTGCTGAAGGGGGATTATTTGGCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAgaga                                                                                   <  1:748010/149‑1 (MQ=255)
                                                          ttAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGAt                                                                                  >  1:381366/1‑149 (MQ=255)
                                                                tgtgCTGAAGGGGGATTATTGGTCGTGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCATGGGGTACAGAGATAGTTGa                                                                            >  1:579842/1‑149 (MQ=255)
                                                                tgtgCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGa                                                                            <  1:89721/149‑1 (MQ=255)
                                                                  ttcTGAAGGGGGATTATTGCTCATGACCGTCTCCTTGCGGTGTGAGGTATGACAAGGGTAGACGCGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaa                                                                          <  1:353868/147‑1 (MQ=255)
                                                                  ggctctaGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCCGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaa                                                                          <  2:571301/143‑1 (MQ=255)
                                                                    cTGAAGGGGGATTATTAGTCATGATCGTCTCCTTGCGGTGTGAGTTATAACAAAGGTAGACGTGACTTGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGaaaaa                                                                        >  2:553456/1‑149 (MQ=255)
                                                                     tGAAGTGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                       <  1:3604/149‑1 (MQ=255)
                                                                     tGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                       <  1:140617/149‑1 (MQ=255)
                                                                     tGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAg                                                                       <  2:745274/149‑1 (MQ=255)
                                                                         gggggATTATTGTTTATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGTCGTCCCGGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGc                                                                   <  1:768874/149‑1 (MQ=255)
                                                                             gATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAAc                                                               >  1:331150/1‑149 (MQ=255)
                                                                             gATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGAAGACGTGACTGGCGAAATCTTCGCCAGCCGGTAACAGCTTTCCGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTCGAAAACGCGGCTAAc                                                               >  1:317583/1‑149 (MQ=255)
                                                                                  ttGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAAttt                                                          <  1:139593/149‑1 (MQ=255)
                                                                                    ggTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGc                                                        <  1:11864/149‑1 (MQ=255)
                                                                                     gTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAACAGCGGCTAACAGTTTGcc                                                       >  1:488891/1‑149 (MQ=255)
                                                                                          gATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCg                                                  >  2:752815/1‑149 (MQ=255)
                                                                                            tCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACTGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGtt                                                >  1:63496/1‑149 (MQ=255)
                                                                                                    ttGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATAACGAATTTACGGGCCAGGGTTACAGAGATAGTTGAAAAAGCGGCTCACAATTTGCCAGCCGTTGTGGAAat                                        >  1:547984/1‑149 (MQ=255)
                                                                                                    ttGCGGTGTGAGGTATAACAAGGGGAGACGTGACTGGCGAAAGCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAGTTTACGGGCTAGGGTTACAGAGATATTTGAAAAAGCGGCTAACATTTTGCCAGCCGTTGTGGAAat                                        <  2:220689/149‑1 (MQ=255)
                                                                                                      gCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAAGTGGTTTCCAACCGTTGTGGaaagta                                      >  2:197317/1‑146 (MQ=255)
                                                                                                      gCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAatga                                      >  1:666851/1‑149 (MQ=255)
                                                                                                       cGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGACAAAGCGGCTCACAATTTGCCCGCCGTTGTGGAAatgat                                     <  2:189138/149‑1 (MQ=255)
                                                                                                         gtgtGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACCGCTGTACGACAGCTGTTACGAATTTACGGGCAAGGGCTACAGAGATAGTTGAAAAAGCGGCTCACAATTTGCCAGCCGTTGTGGAAAtggtga                                   >  2:471897/1‑149 (MQ=255)
                                                                                                         gtgtGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCACGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCCGCCGTTGTGGAAAtgatga                                   >  2:293169/1‑149 (MQ=255)
                                                                                                         gtgtGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCAGTAACAGCTTTACGACAGCGATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTAAAAAAGCGGCTAGCAATTTGCCAGCCGTTGTGGAAAtgatga                                   >  2:38318/1‑149 (MQ=255)
                                                                                                         gtgtGAGGGATAACAAGGGTAGACGTGACTGGCGAAAACTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAATAATTAGCCAGCCGTTGTGGAAAtgatga                                   <  1:511547/149‑1 (MQ=255)
                                                                                                            tgAGGTATAACAAAGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTa                                >  1:535959/1‑149 (MQ=255)
                                                                                                             gAGGTATAACAAGGGTAGACGTGACAGGAAAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTAGAGGGATAGTTGAAAAAGCGGCTAACAAGTTGCCAGCGGTTGTGGAAATGATGATTAc                               >  1:255635/1‑149 (MQ=255)
                                                                                                              aGGTATAACAAGGGTAGAGGTGACTGGCGAAATCTTCGCCAGGCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACa                              <  1:218692/149‑1 (MQ=255)
                                                                                                              aGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACa                              >  1:754602/1‑149 (MQ=255)
                                                                                                                gTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAACGAATAACAATTTGCCAGCCGTTGTGGAGATGATGAGTACaac                            >  2:257944/1‑147 (MQ=255)
                                                                                                                    aaCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTg                        >  1:626749/1‑149 (MQ=255)
                                                                                                                    aaCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTg                        >  1:773509/1‑149 (MQ=255)
                                                                                                                     aCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgc                       >  1:292701/1‑149 (MQ=255)
                                                                                                                     aCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgc                       <  1:250692/149‑1 (MQ=255)
                                                                                                                     aCAAGGGTAGACGTGACTGGCGAAAGCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgc                       <  1:247047/149‑1 (MQ=255)
                                                                                                                      cAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGGTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgcg                      >  1:83340/1‑149 (MQ=255)
                                                                                                                      cAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTgcg                      >  2:591237/1‑149 (MQ=255)
                                                                                                                         gggTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGa                   <  1:525099/149‑1 (MQ=255)
                                                                                                                         gggTAGACGGGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGa                   <  2:255313/149‑1 (MQ=255)
                                                                                                                            tAGACGTGACTGGCGAGATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATACCGACTTTACGCGCAAGGGGTTCAGAGATAGTTGACAATTCTGTGAACAGGCGGCCATCCGTTTTGGAAATCATGATTACAACGCTCcgcgcgca                >  2:552473/1‑145 (MQ=255)
                                                                                                                                cGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCg            >  1:686936/1‑149 (MQ=255)
                                                                                                                                   gACTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCca         >  2:373633/1‑149 (MQ=255)
                                                                                                                                    aCTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCcac        >  1:714116/1‑149 (MQ=255)
                                                                                                                                    aCTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCcac        >  1:549234/1‑149 (MQ=255)
                                                                                                                                    aCTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCCCGATCAATCGcccc        >  2:81488/1‑147 (MQ=255)
                                                                                                                                       ggCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAAt     >  2:496035/1‑149 (MQ=255)
                                                                                                                                          gAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATttc  >  1:753267/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTC  >  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: