Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,776,2630TG41.2% 147.0 / 12.6 80V226G (GTG→GGG) ydiFputative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base T (12/35);  new base G (33/0);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.05e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTTG  >  NC_000913/1776114‑1776412
                                                                                                                                                      |                                                                                                                                                       
cAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTg                                                                                                                                                        <  2:347852/151‑1 (MQ=255)
 aaaGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGc                                                                                                                                                       <  2:36311/151‑1 (MQ=255)
        ttGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGaaaa                                                                                                                                                 >  1:415352/1‑150 (MQ=255)
         tGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATTCAGGGGGAGgaaa                                                                                                                                                 >  2:1174411/1‑149 (MQ=255)
             ccAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTTTGGTGATGGAGGGGCAGAAAAAgggt                                                                                                                                            >  2:81570/1‑148 (MQ=255)
               aGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAg                                                                                                                                         <  1:1045611/151‑1 (MQ=255)
               aGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAg                                                                                                                                         <  2:311845/151‑1 (MQ=255)
                     tGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGGTaaaaaa                                                                                                                                      >  2:798852/1‑148 (MQ=255)
                     tGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGCATTATGGTGATTGCGGGGAAGAAAAAGGGTAAGaaaaac                                                                                                                                   >  1:409389/1‑148 (MQ=255)
                       ccTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCCGGGGCCGAAAAAGGGTAAGAAAACCAc                                                                                                                                 >  1:728473/1‑151 (MQ=255)
                       ccTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCAc                                                                                                                                 <  2:813511/151‑1 (MQ=255)
                        cTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTGTTGGGATGATGGAGGGGTAGAAAATGGGTAAAAAAACCAc                                                                                                                                 >  1:163721/1‑150 (MQ=255)
                        cTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATTGCGGGGGAGAAAAAGGGTAAAAAAGCCAc                                                                                                                                 >  1:184307/1‑150 (MQ=255)
                         ttCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGGAGAAAATGGTTAAGAAAGCCACGc                                                                                                                               >  1:861436/1‑151 (MQ=255)
                          tCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACg                                                                                                                                <  1:571895/149‑1 (MQ=255)
                           cATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGGTAAGAAAGCCAc                                                                                                                                 >  1:1057025/1‑147 (MQ=255)
                            aTTCGCGCTACCACCTGCGACAGTGAAGGCTACGCTACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGc                                                                                                                            <  1:10059/151‑1 (MQ=255)
                            aTTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGGAGGAAAAGGGTTAGAAAGGCACGCTg                                                                                                                             >  2:1153702/1‑150 (MQ=255)
                              tCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGGAGGGGCAGAAAATGGGTAAGAAAGGCACGCTGCAt                                                                                                                          >  2:316545/1‑151 (MQ=255)
                              tCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGGTTTTGATGAAGCCGGGGCAGAAACTTGTTaaaaaaagcaagc                                                                                                                               >  2:848521/1‑139 (MQ=255)
                                  gcTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTa                                                                                                                      <  2:522161/151‑1 (MQ=255)
                                    tACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTTTTTTGATGATTCAAGGGGCGAAAAAGGGTTAGAAAACCCACCTGCCACCTaa                                                                                                                     >  2:637611/1‑150 (MQ=255)
                                     accaccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTaaa                                                                                                                    <  2:578258/150‑1 (MQ=255)
                                        accTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGAAGAGGCAGGGGTAGAAAAAGGGTAAGAAAGACACGCTGCATCCCAAATCt                                                                                                                 >  1:251747/1‑150 (MQ=255)
                                        accTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTAATGTGGTGGTGCAGGGGGAGAAAATGGTTAAGAAAACCACGCTGCAGCCTAAATc                                                                                                                  >  1:309215/1‑149 (MQ=255)
                                          cTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGTTAAGAAAACCACGCTGCATCCTAAAACTGt                                                                                                               >  1:490815/1‑150 (MQ=255)
                                            gCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTcc                                                                                                             <  2:762802/150‑1 (MQ=255)
                                            gCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCg                                                                                                            <  2:230001/151‑1 (MQ=255)
                                               aCAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTAt                                                                                                        <  2:154337/151‑1 (MQ=255)
                                                   tGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGAAATGTGATGTTTCAGGGGAAGAGAATGGATAAGAAAGAGACGCTGGCTACTTAATCTGTCGGGAttt                                                                                                      >  2:358656/1‑148 (MQ=255)
                                                      aGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTTTTGTGATGATTCAGGGGCAGAAAAAGGTTAATAAAACCCCCCCTCCGCCTAAATCTGTCAGTATTGCgggg                                                                                                 >  2:520648/1‑150 (MQ=255)
                                                       ggCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGCAGGGGCAGGAAATGGGTTAGAAAAGCACGCTGCATCCCAAAACTGTCCGTATTCCggg                                                                                                  >  2:258871/1‑149 (MQ=255)
                                                        gCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTa                                                                                               <  1:597827/151‑1 (MQ=255)
                                                            cGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCtg                                                                                           <  2:334523/151‑1 (MQ=255)
                                                             gCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGGTAAGAAAGCCACCCTGCATCCTAAATCTGTCCGTATTTCGGGGTATTtgg                                                                                          >  1:1018332/1‑151 (MQ=255)
                                                             gCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGGAGAAAATGGGTTAGAAAGCCACGCTGCATCCTAAATCTTTTCGGATTCCGGGTTAttt                                                                                            >  2:708096/1‑147 (MQ=255)
                                                              ccACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGGTAAGAAAACCACGCTGCATCCTAAATCTGTCCGTATTTCGGGGTATCtgg                                                                                          >  1:215893/1‑150 (MQ=255)
                                                              ccACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTTTGGTGATGCAGGTGGAGAAAAAGGGTTAGAAAACCACGCTGCATCCTAAAACTTTTCGTATTTCGGGTTATTtgg                                                                                          >  1:499013/1‑150 (MQ=255)
                                                               cACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGTTAAGAAAGACCCGCTGCATTCTAAATCTGTTCGTATTCCGGGGTATCtggtg                                                                                        >  2:211777/1‑151 (MQ=255)
                                                                aCTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAAAGGGTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGGATTCCGGGGTATCtggggg                                                                                       >  1:379963/1‑151 (MQ=255)
                                                                   tttGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAATCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGat                                                                                     <  2:573213/150‑1 (MQ=255)
                                                                     tGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGatat                                                                                   <  2:110853/150‑1 (MQ=255)
                                                                     tGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGatat                                                                                   <  1:150349/150‑1 (MQ=255)
                                                                      gAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATAtt                                                                                  <  2:411485/150‑1 (MQ=255)
                                                                        aGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGTTTTTGTCGGGGGAGAAAAAGGTTAAGAAAATCACGCTGCATCCTAAATCTGTACGTATTCCGGGGTATCTTGGGGCTAtttt                                                                                >  2:1144968/1‑148 (MQ=255)
                                                                        aGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATtgtg                                                                               <  2:892403/151‑1 (MQ=255)
                                                                        aGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATtgtg                                                                               <  1:563137/151‑1 (MQ=255)
                                                                             aGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATTCAGGGGCAGAAAATTGTTAAGCAAGCCACGCCGGATCCTAAATCTCTCCCTATTTCCGGGTTTGtggtggggatt                                                                                  >  2:977845/1‑138 (MQ=255)
                                                                             aGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGtggtgg                                                                           <  1:668039/150‑1 (MQ=255)
                                                                              ggTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTc                                                                         <  2:433414/151‑1 (MQ=255)
                                                                                   tGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGCAGGTGGAGAAAATGGTTAAGAAAACCACGCTGCCTCCTAAATCTGTCCGTGTTTCGGGGTTTCTGGTGGAGATTTTGGTGGTCGTTc                                                                     >  1:146850/1‑150 (MQ=255)
                                                                                      aTCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTTATGATGCAGGTGGAGAAAAAGGTTAAGAAAACCACGCTGCATCCTAAATCTGTCCGGATTCCGGGGTTTCTGGTGGGTATTGTGGGGGGCGATCCgg                                                                  >  1:289160/1‑150 (MQ=255)
                                                                                       tctcGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGTTAAGAAAACCACGCCGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGGGGATATTGTGGTGGTCGATCCgg                                                                  >  1:239918/1‑149 (MQ=255)
                                                                                          cGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGa                                                                 <  1:144788/147‑1 (MQ=255)
                                                                                          cGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGAt                                                                <  2:625963/148‑1 (MQ=255)
                                                                                          cGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATc                                                               <  1:196955/149‑1 (MQ=255)
                                                                                             cGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTTGTGATGCAGGGGGAGAAAAAGGGTAAGAAAGCCACCCCGCATCCTAAAACTGTCCCTATTCCCGGTTATTTGGGGGGTATTTTGGTGGTCGCTCCCGATCaaa                                                           >  2:128884/1‑149 (MQ=255)
                                                                                             cGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGGGATGATGGAGGGGGAGAAAAAGGTTTAAAAAAGCACCCTGCATTCTAAATCTTTCCCTATTCCGGGTTATCTGGTGGATATTGTGGTGGGCGATCCGGGTCaaa                                                           >  1:1007919/1‑149 (MQ=255)
                                                                                                  tGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGGTATTGTGGGGGTCGATCCGGATCAAACCCaa                                                      >  1:499114/1‑149 (MQ=255)
                                                                                                    gTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGt                                                  <  1:182776/151‑1 (MQ=255)
                                                                                                       aTTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGTAGGGGCAGAAAATGGGTAAGAAAACCCCCCTGCCTCCTAAATCTGTCCGCATTCCGGGGTTTTTGGGGGATATTGTGGTGGTCGATCCCGGACAAAACCAACTggtt                                                >  2:896241/1‑147 (MQ=255)
                                                                                                       aTTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATg                                               <  2:821985/151‑1 (MQ=255)
                                                                                                          gCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGc                                             <  2:297028/150‑1 (MQ=255)
                                                                                                           cccAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGGCGAAAATTGTTAAGAAAAGCACGCTGCAACCTAAATCTTTCCGGATTCCCGGGTATTTGGGTGGTATTTTGGGGGTCGATTCGGATCAAAACCCAATTgattgcgg                                            >  1:184056/1‑144 (MQ=255)
                                                                                                           cccAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGGTAAGAAAACCACGCTGCATCCTAAATCTGTTCGTATTCCGGGTTTTCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCgg                                           >  1:536594/1‑151 (MQ=255)
                                                                                                            ccAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGt                                          <  1:1041048/151‑1 (MQ=255)
                                                                                                              aGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAATGGGTAAGAAAACCACGCTGCATCCTAAATCTGTCCGTATTTCGGGTTATCTGGGGGGTATTTTGGGGGTCGATCCGGATCAAAACCCACTGTATGGCgg                                           >  1:43234/1‑148 (MQ=255)
                                                                                                               ggCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCa                                       <  2:55079/151‑1 (MQ=255)
                                                                                                                      cacaATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTa                                <  1:972097/151‑1 (MQ=255)
                                                                                                                        caATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAAc                              >  1:46602/1‑151 (MQ=255)
                                                                                                                             aCGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGGATTCCGGGTTATCTGGGGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGGGCACCGGTTTACCGCtt                         >  1:321791/1‑151 (MQ=255)
                                                                                                                               ggcggTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCttt                        >  1:410678/1‑150 (MQ=255)
                                                                                                                               ggcggTATTGTGATGATGCAGGTGCAGAAAATGGCTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCttt                        >  2:339125/1‑150 (MQ=255)
                                                                                                                                   gTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTc                   >  1:419459/1‑151 (MQ=255)
                                                                                                                                          gatgatGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGa             <  1:630253/150‑1 (MQ=255)
                                                                                                                                          gatgatGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGAc            >  1:266064/1‑151 (MQ=255)
                                                                                                                                           atgatgCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGa             <  2:598583/149‑1 (MQ=255)
                                                                                                                                             gatgCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTc         <  2:1011382/151‑1 (MQ=255)
                                                                                                                                                gCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCAcc      <  2:77732/151‑1 (MQ=255)
                                                                                                                                                     gTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTtg  >  1:210365/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
CAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTTG  >  NC_000913/1776114‑1776412

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