Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913862,5590GT38.7% 103.6 / 41.1 75H16Q (CAC→CAA) ybiWputative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (14/32);  new base T (29/0);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 1.76e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCAT  >  NC_000913/862416‑862709
                                                                                                                                               |                                                                                                                                                       
gTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCtt                                                                                                                                                  <  1:859927/151‑1 (MQ=255)
 tGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGAGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCTCAATAAACACCAGCGCATTTTTTTGCGCttt                                                                                                                                                 >  1:444816/1‑151 (MQ=255)
  gCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑a                                                                                                                                               <  2:239307/151‑1 (MQ=255)
     aTTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTCCAATATTCACCAGCGCCTTTTTTTGCGCcttt                                                                                                                                                >  1:280691/1‑147 (MQ=255)
        cGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCgg                                                                                                                                         <  2:481025/151‑1 (MQ=255)
        cGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCgg                                                                                                                                         <  2:1091595/151‑1 (MQ=255)
          ccAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTc                                                                                                                                       <  1:1193229/151‑1 (MQ=255)
           cAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTTCAAAATGCCCCAGCGCCATTTTTTGCGCTTT‑TATTCGGGCg                                                                                                                                      >  2:1052268/1‑151 (MQ=255)
            aGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTTCAATATGCACCAGCGCATTTTTTTGGGCTTT‑TATGCgg                                                                                                                                         >  2:947365/1‑147 (MQ=255)
            aGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAAAATTGACCCGCGCAGTTTTTTTCGCGTT‑TATgt                                                                                                                                           >  1:1252081/1‑144 (MQ=255)
             ggTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCt                                                                                                                                    <  2:563527/151‑1 (MQ=255)
              gTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTg                                                                                                                                   <  1:842312/151‑1 (MQ=255)
               tGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACCAAATGCACCCGCGCCTTTTTTTTCGCTTT‑AATTGGGGCGCTg                                                                                                                                   >  1:493682/1‑150 (MQ=255)
                gATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAg                                                                                                                                 <  1:1374803/151‑1 (MQ=255)
                  tGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAAAATGCACCAGCGCATTTTTTTGTGCTTT‑TAAGCGGGTGGTGAGCg                                                                                                                               >  2:368418/1‑151 (MQ=255)
                     gcCAGTGCCCGCGCGCGACGGACTGAGATTGGCTTATCGAGATGTTGTGGATACATCTCGGTATAGTGCTGCGCACGCTGGGTACAGACTGGCGGTTTCACCATGTGCGCCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCgta                                                                                                                             <  1:123448/150‑2 (MQ=255)
                     gcCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAATAATCCCCCGCGCCTTTTTGTGCGCTTT‑TATTCGGTCGCTGAGAgtg                                                                                                                             >  2:988429/1‑150 (MQ=255)
                             cAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTc                                                                                                                    <  1:1180756/151‑1 (MQ=255)
                                cgcgcgACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAAAATGCACCAGCGCCTTTTTTTGGGCTTT‑TATGCGGGCGCTGGGCGGGTCCAGttt                                                                                                                     >  2:1144110/1‑147 (MQ=255)
                                      aCGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAATATGCACCAGCGCATTTTTTTGCGCTTT‑TATTCGGGCGCTGGGCGTGTCCCGTTTCAGTGTgg                                                                                                             >  2:1259879/1‑149 (MQ=255)
                                      aCGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAATATGCACCAGCGCATTTTTTTGCGCTTT‑TATGCGGGCGCCGAGCGTGTGCAGGTTCAGGGTgg                                                                                                             >  1:745756/1‑149 (MQ=255)
                                        gTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAATATGCACCAGCGCACTTTTTTGCGCTTT‑AATTCGGTCGCTGGGCGGGTCCCGTTTCCGTGTGGTCAt                                                                                                         >  2:964922/1‑151 (MQ=255)
                                          aCCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATg                                                                                                        <  1:1033453/150‑1 (MQ=255)
                                           ccGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACCATATGCACCAGCGCAATTTTTTGGGCTTT‑AATATGGTGGCCGAGCGTGTTAAGGTTCAGTGCGGTCCTgg                                                                                                       >  2:714538/1‑150 (MQ=255)
                                           ccGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACCAAATGCACCCGCGCCATTTTTTGCGCTTT‑TATGCGGGGGCCGAGGGTGGGCCGCTTTAGTGTgggcgttg                                                                                                       >  1:288496/1‑144 (MQ=255)
                                           ccGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAATATGCACCAGCGCATTTTTTTGGGCTTT‑AATGCGGGCGCTGGGCGGGTCCAGTTTCCGTGTgggcct                                                                                                         >  2:1407727/1‑144 (MQ=255)
                                                aTCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAAAAAACAACAGAGGAATTTTTTTTGCTTT‑AAATCGGGGGCTGAGCGGGTCCCCTTTTAATGTGGGCATGggtgtt                                                                                                  >  1:1019014/1‑150 (MQ=255)
                                                   ggCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGGTTTAAAATGTACACCAGCGCATTTTTGTGCGCGTT‑AATGTGGGCGCTGAGCGGGTGCAGGTTCAGTGGGGTCGTGgttgtt                                                                                                  >  2:868730/1‑147 (MQ=255)
                                                    gCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCt                                                                                             <  1:54429/151‑1 (MQ=255)
                                                    gCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGCTTCACAATTTTCACCACTGCATTTTTTGGGGCTTTAAATTCGG‑CGCCGAGCGCGTCCACATTTAGGGTgggtccgggt                                                                                                     >  2:1192457/1‑135 (MQ=255)
                                                       tATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCg                                                                                           <  2:926180/150‑1 (MQ=255)
                                                           gagaTGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCCCAATATTCACCAGTGCCGTTTTTTGCGGCTT‑TATGCGGGCGGTTGGCGTGTTCCCGTTCACGGGGGTCATGGTTGTTTTTCTTGTac                                                                                        >  1:668332/1‑148 (MQ=255)
                                                           gagaTGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGGCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTaa                                                                                        <  1:512622/149‑1 (MQ=255)
                                                            agaTGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGt                                                                                     <  2:99737/151‑1 (MQ=255)
                                                             gaTGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTc                                                                                    <  2:1113177/151‑1 (MQ=255)
                                                              aTGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGGTTCTCAAAATGCACCCGCGCCTTTTTTTGGGCTTT‑TATGCGGGCGCCGAGCGTGTCCCGGTTCTGTGTGGGCATGGTTGTTTTCCTCGTAAGGGcg                                                                                   >  1:1114493/1‑151 (MQ=255)
                                                               tgttgtTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCCTTTTTTTGCGCTTT‑TATGCGGTCGCTGAGCGTGTCCCGTTTTTGTGTGGGCATGGTTGTTATCCTCGTAAgg                                                                                      >  2:957235/1‑147 (MQ=255)
                                                                 ttgttgATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAAAATGCACCCACGCATTTTTTTGCGCTTT‑AATGGGGGCGCTGAGCGTGTCCAGTTTTAGTGTGGGCATGGTTGTTATCCTCGTAAGGTCGCgg                                                                                >  2:87547/1‑151 (MQ=255)
                                                                  tgttgATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCgg                                                                                <  1:659789/150‑1 (MQ=255)
                                                                   gttgATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAATATTTCCCAGCCCCTTTTTTTTGGCTTT‑TATGCAGGCGCCGGGAGGGGCCCGTTTCAGTGGGGGCCTGGTGGTTTTTCTCCTAAAGTc                                                                                    >  1:1365107/1‑145 (MQ=255)
                                                                   gttgATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGtt                                                                              <  2:1108958/151‑1 (MQ=255)
                                                                       aTACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGGTTCAAAATATGCAACAGCGCATTTTTTTGTGCTTT‑TATTGGGGGGCTTCGGGTGTGCAGTCTTAGTGGGGGGATGGTTGGTATTCTCGCAAGGTGGCgggt                                                                              >  1:1271419/1‑145 (MQ=255)
                                                                          cATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCTCAATAAGACCCCGCGCACTTTTGTGGGCTTT‑AATGCGGGGCCTGAGCGTGTCCAGATTTAGTGGGGGTCTGGGTGGTTTCCTCGTCAGGGCGCGGGTAACCCtt                                                                       >  1:48605/1‑151 (MQ=255)
                                                                          cATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCtt                                                                       <  2:1191014/151‑1 (MQ=255)
                                                                             ctcGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTCCAATATGCACCAGCGCATTTTTTTGCGCTTT‑TATGCGGGCGCTGAGCGTGTCCCGTTTTTGTGTGGGCATTGTTGTTATCCTCGTAAGGGCGCGGGTTACCCTTTCt                                                                    >  2:1047203/1‑151 (MQ=255)
                                                                              tcGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGGTTCTCAATATGCACCCACGCCTTTTTTTTGGCTTT‑TTTGCCGGCGCTGCGCGGGTCCCGTGTCCCGGTGGTGTCGGGTTTTTTCCCTGTAAGGGCGCGGGGTACCCttt                                                                      >  2:1190775/1‑148 (MQ=255)
                                                                               cGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCt                                                                    <  1:1112861/149‑1 (MQ=255)
                                                                                      gTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTACAATATGCACCAACGCATTTTTTTGCGGCTT‑TATGTGGGCGGGGAGCGTGTCCCAGTTCCGTGTGGGCATGGTTGTTATCCTCTTGAGGGCGCGGTTAATCCTCTCTtgcgggcac                                                           >  1:243141/1‑150 (MQ=255)
                                                                                       tgctgcGCGCGCTCGGTACAGACTGGCGGTTTTACAATATGCCCCAGCGCCTTTTTTTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTTAGTGTGGGCATGGTTGTTATCCTCGTTAAGTGGCGGGTTACCCTTTTTGGCAGGCAt                                                           >  2:760682/1‑150 (MQ=255)
                                                                                          tgcGCGCGCTCGGTACAGACTGGCGGTTTCCCAAAATGCACCCGCGCATTTTTGTGCGCTTT‑TATGCGGGCGCTGCGCGTGGCCAGGTTTAGTGTGGGCCTGGGTGTTTTTCTCGTAAGGTCGCGGGTAAACCCTTCTTGCAGGGATAct                                                        >  1:1008440/1‑150 (MQ=255)
                                                                                           gcgcgcgcTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATActg                                                       <  1:166150/150‑1 (MQ=255)
                                                                                           gcgcgcgcTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATActg                                                       <  2:1089032/150‑1 (MQ=255)
                                                                                           gcgcgcgcTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCACGCATActg                                                       <  2:1000786/150‑1 (MQ=255)
                                                                                            cgcgcgcTCGGTACAGACTGGCGGTTTCACTATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATActgct                                                     <  2:782392/151‑1 (MQ=255)
                                                                                                cgcTCGGTACAGACTGGCGGTTTCACCAAATTCAACACGGCCCTTTTTTGCGCCTT‑CAAATGGGGGGTGTGCGAGTGCAGGTTTTGTGTGGGCCTGGTTGTTTTTCTCCGAAGGGCGCGGGTAACCCCTTCTGGGAGGCATActcattg                                                   >  1:1414597/1‑144 (MQ=255)
                                                                                                  cTCGGTACAGACTGGCGGTTTTACAAAATGCACCCACGCATTTTTTTTCGCTTT‑AATGGGGTCGCTGAGCGTGTCCCGTTTCAGTGTGGTCCTGGTTGTTTTTCTCGTAAAGGCGCGGGTAAACCCTTTTTGCAGGCATACTGCTggg                                                  >  1:767338/1‑148 (MQ=255)
                                                                                                   tCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTTTGCGCTTT‑AATGAGGTCGCTGAGCGGGTCCAGTTTCAGTGGGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGGaa                                               >  2:1123990/1‑150 (MQ=255)
                                                                                                    cGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCaaa                                              <  1:643548/150‑1 (MQ=255)
                                                                                                              tGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGtt                                   >  2:269102/1‑151 (MQ=255)
                                                                                                               ggCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTc                                  >  2:1347498/1‑151 (MQ=255)
                                                                                                                 cGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTg                                <  2:41901/151‑1 (MQ=255)
                                                                                                                     ttCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGgcgc                            >  1:138293/1‑151 (MQ=255)
                                                                                                                        acaATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATc                         <  1:1255710/151‑1 (MQ=255)
                                                                                                                         caATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATc                         >  1:119152/1‑150 (MQ=255)
                                                                                                                          aaTATGCACCAGCGCTTTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGCGTCCAGTTTCTGTGTGGTCGTGGTTGTTATCCTCGCTAGGTCTCGGTTAACCCTTTCTGGCAGGCTTACTGTTGCGCAAAGTCGAGCAGGTCTGgcgc                            >  1:478503/1‑146 (MQ=255)
                                                                                                                          aaTATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCaa                       <  2:1264099/151‑1 (MQ=255)
                                                                                                                           atatGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAg                      <  1:427402/151‑1 (MQ=255)
                                                                                                                            tatGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGc                     >  2:359277/1‑151 (MQ=255)
                                                                                                                                  ccAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGtttt               <  2:1305892/151‑1 (MQ=255)
                                                                                                                                   cAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGttttt              >  2:1108779/1‑151 (MQ=255)
                                                                                                                                     gcgcATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTcc            >  1:1388100/1‑151 (MQ=255)
                                                                                                                                         aTTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCgggg        >  2:21848/1‑151 (MQ=255)
                                                                                                                                          tttttGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGc       >  1:338059/1‑151 (MQ=255)
                                                                                                                                           ttttGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGc       <  2:340088/150‑1 (MQ=255)
                                                                                                                                            tttGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCATTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCg      <  2:584397/150‑1 (MQ=255)
                                                                                                                                               gtgCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCAt  <  2:555993/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
GTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTCACAATATGCACCAGCGCATTTTTGTGCGCTTT‑AATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCAT  >  NC_000913/862416‑862709

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