Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,671,4150CT43.4% 97.7 / 43.6 76K160K (AAG→AAAyhjBputative DNA‑binding transcriptional response regulator
Reads supporting (aligned to +/- strand):  ref base C (40/3);  new base T (0/33);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.99e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGCC  >  NC_000913/3671266‑3671559
                                                                                                                                                      |                                                                                                                                                
tCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTcc                                                                                                                                                  <  1:855540/150‑1 (MQ=255)
  ctgctgGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCtt                                                                                                                                                >  2:744043/1‑150 (MQ=255)
       ggctgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTAttt                                                                                                                                            <  2:529178/149‑1 (MQ=255)
        gctgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGAt                                                                                                                                         <  2:1404906/151‑1 (MQ=255)
        gctgTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGAt                                                                                                                                         >  2:66592/1‑151 (MQ=255)
             aaaaTTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTcgc                                                                                                                                     <  1:113518/150‑1 (MQ=255)
                aTTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTcgcccgc                                                                                                                                 >  2:366700/1‑151 (MQ=255)
                   aTCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCTTTATTTGATTCGCccgccg                                                                                                                               <  2:1240470/150‑1 (MQ=255)
                     ccTCTTCAGGAGGAGATATTTAACATCATTGCCGCCGGGGTGCGTTTTTTCCCTGCCAGACGGCGATACGGGGTTTCCAGATGCGTTTTTCCCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCTTTTTTTGATTCTCCcgccgc                                                                                                                              <  2:595336/149‑1 (MQ=255)
                      ctctTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTa                                                                                                                            >  2:992145/1‑150 (MQ=255)
                            aggaggAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGt                                                                                                                     >  2:816476/1‑151 (MQ=255)
                                ggagATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTaa                                                                                                                  >  1:759482/1‑150 (MQ=255)
                                  agaTATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAAt                                                                                                               >  2:1195924/1‑151 (MQ=255)
                                      atTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTcc                                                                                                           >  1:1227235/1‑151 (MQ=255)
                                      atTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTcc                                                                                                           >  2:308927/1‑151 (MQ=255)
                                       tttAACATCTTGCCCGCTGGGGTGCTTTTTTTCCCTCCCAGAGGGAGATACAGGGATTCCAGATGCGTTTTTCCCTTCCCGGTACTGATTTTCAGCGCTCTGCCGTTCTCTTTTTTTGTTCCGCCCGCCGCTAACATGGTTAAAATCTcc                                                                                                           <  1:851458/150‑1 (MQ=255)
                                       tttAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTccc                                                                                                          >  2:1287205/1‑151 (MQ=255)
                                         tAACATCATTGCCGCCTGGGAGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGc                                                                                                        >  2:227102/1‑151 (MQ=255)
                                         tAACATCATTGCCGCCGGGGTGGTTTTTTTCCCTCCCAGACGGCGATACGGGGATTCCAGATGCGTTTTTCCCGTTCCGGTACTGATATTCGGCGCTCTGCCGATCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCg                                                                                                         <  1:230036/150‑1 (MQ=255)
                                          aaCATCATTGCCGCCGGGGTGCGTTTTTCCACTCCCAGCCGGCGATACAGGGATCCCAGATGCGTTTTTCCCTTTCCGGTCCTGATATCCGGCGCTCTGCCGTTCTCTTTTTTTGTTCCGCCCGCCGCTAACATGGTTAAAATCTCCCGc                                                                                                        <  2:529197/150‑1 (MQ=255)
                                            catcatTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTgg                                                                                                     >  2:983606/1‑151 (MQ=255)
                                              tcatTGCCGCCTGGGTGCGTTTTTTCACTCCCAGACGGCGATACAGGGATTCCAGATGCTTTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCTTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggc                                                                                                    <  1:608923/150‑1 (MQ=255)
                                              tcatTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcg                                                                                                   >  2:686713/1‑151 (MQ=255)
                                               catTGCCGCGTGGGTGCTTTTTTTCACTCCCAGACGGCGATACAGGTATTCCAGATGCGTTTTTCCCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCTTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggc                                                                                                    <  1:1100421/149‑1 (MQ=255)
                                               catTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgg                                                                                                  >  1:1021972/1‑151 (MQ=255)
                                                 ttGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgggc                                                                                                >  1:340254/1‑151 (MQ=255)
                                                  tGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgggcg                                                                                               >  1:493178/1‑151 (MQ=255)
                                                  tGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgggc                                                                                                >  2:1072308/1‑150 (MQ=255)
                                                    ccgctgggttgctttttttCCCTCCCAGAGGGCGATACAGGGATTCCAGTTGCTTTTTTCCCGTTCGGTTATTTATATCCAGCGCTCTGCCGTTCTCTTTATTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTggcgggcg                                                                                               <  1:642066/135‑1 (MQ=255)
                                                         tGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACg                                                                                        >  1:44788/1‑151 (MQ=255)
                                                          ggttgtttttttcCACTCCCAGCCGGCGATACAGGTTTCCCAGATGCGTTTTTCCCGTTCCGGTACTGATTTTCAGCTCTCTGCCGTTCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAAc                                                                                         <  2:333104/141‑1 (MQ=255)
                                                          gggTGCTTTTTTTCCCTCCCAGACGGCGATACAGGGATCCCAGATGCGTTTTTCCCGTTCCGGTACTGATATTCAGCGCTCTCCCGATCTCTTTATTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGa                                                                                       <  1:398574/151‑1 (MQ=255)
                                                          gggTGCGTTTTTCCCCTCCCAGAGGGCGATACGGGGATCCCAGTTGCTTTTTTCCCTTCCCGTTACTGTTATCCAGCGCTCTCCCGTTCTCTTTTTTTTTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCtt                                                                                            <  2:839332/146‑1 (MQ=255)
                                                          gggTCCGTTTTTTCCCTTCCAGACGGCGATACGGGGATTCCAGATGCTTTTTTCCCGTCCCGGTCCTGATATTCAGCGCTCTGCCGATCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACg                                                                                        <  1:1251402/150‑1 (MQ=255)
                                                              gCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTg                                                                                   >  1:1199397/1‑151 (MQ=255)
                                                               cGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTg                                                                                   >  1:989644/1‑150 (MQ=255)
                                                                gATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTg                                                                                   >  2:881152/1‑149 (MQ=255)
                                                                gATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTg                                                                                   >  1:1082482/1‑149 (MQ=255)
                                                                  ttttcCCCTCCCAGACGGCAATACAGGGTTTCCAGATGCGTTTTTCCCGTCCCGGTACTGATATCCAGCGCTCTGCCGTTCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTga                                                                                  <  1:1066077/148‑1 (MQ=255)
                                                                   ttttCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTg                                                                                   >  2:1367247/1‑146 (MQ=255)
                                                                          ttCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATg                                                                       >  2:1414993/1‑151 (MQ=255)
                                                                            cgaggcGGAGATACAGGTATCCAAGTTTCGTTTTCACCGTTCCGGTATTGATATACCGCGCTCTGCCGTTCTTTTTTTTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAAt                                                                        <  2:544787/143‑1 (MQ=255)
                                                                            ccAGCCGGCAATACAGGGATCCAAGATGCGTTTTTCCCGTCCCGTTACTAATATCCAGCGCTCTGCCGTTCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTaa                                                                         <  1:1011726/147‑1 (MQ=255)
                                                                              aGACGGCGATACGGGGATTCCAGATGCGCTTTTCCCGTTCCGTTACTGATATTCAGCGCTCTCCCGATCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGt                                                                      <  2:1082382/148‑1 (MQ=255)
                                                                              aGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCt                                                                   >  1:1437182/1‑151 (MQ=255)
                                                                               gACGGCGATACGGGGATCCAAGTTGCGCTTTTCCCGTCCCGTTACTTATTTTCAGCGCTCTGCCTTTCTTTTTTTTTTTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTc                                                                     <  2:577192/148‑1 (MQ=255)
                                                                                  ggggATACGGGGATTCCAGATGCTTTTTTCCCGTCCCGTTATTGATATCCAGCCCTTTGCCGTTCTCTTTTTTTGTTTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCtttt                                                                <  1:895193/147‑1 (MQ=255)
                                                                                  ggCGATAGGGGAATTCCAGGGGCTTTTTTCCCTTCCCGGTAGTAATATCCAGGGCTCGCCCGGTCTCTTTTTTTGTTCCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTaa                                                                         <  2:1115829/141‑1 (MQ=255)
                                                                                  ggCGATACGGGGTTTCCAGTTGCTTTTTTCCCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCttt                                                                 <  1:66602/149‑1 (MQ=255)
                                                                                  ggCAATACAGGGATTCCAGATGCGTTTTTCCCGTCCCGGTATTGATATCCAGCCCTCTGCCGTTCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTcc                                                                    <  1:1356626/146‑1 (MQ=255)
                                                                                   gCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTcc                                                              >  1:97707/1‑151 (MQ=255)
                                                                                    cGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCg                                                             >  2:966583/1‑151 (MQ=255)
                                                                                     aaTACAGGGATCCCAGATGCTTTTTTCCCGTCCCGGTACTGATTTTCGGCCCTCGCCCGTTCTCTTTTTTTGATCCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCgg                                                            <  2:989550/150‑1 (MQ=255)
                                                                                      aTACGGGAATCCCAGATGCTTTTTTCCCGTTCCGTTACTATTATTCATCGCTCTGCCTTTCTCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTcc                                                              <  1:1214829/148‑1 (MQ=255)
                                                                                      aTACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGc                                                           >  1:986978/1‑151 (MQ=255)
                                                                                        acggtgttcccAGATCCGTTTTTCCCTTTCGGGTACATAGATACGGCGCTCGGCCGTTCTCTTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgt                                                         <  2:363963/142‑1 (MQ=255)
                                                                                        acgggtaTTCCAGTTCCTTTTTTCCCTTCCCGGTACTAATATCCAGCGCTCTGCCTTTCTCTTTTTTTTATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCgg                                                            <  2:256420/142‑1 (MQ=255)
                                                                                        aCAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgt                                                         >  2:869518/1‑151 (MQ=255)
                                                                                            ggtattcaagtggcgttttTCCCTTCCGGTTATTTATATACGGCGCTCTCCCTATCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGcgtcg                                                       <  1:771647/135‑1 (MQ=255)
                                                                                               tcccagactcttttttCCCGTTCCGGCCCTTTTCTTCCGCCCTCGCCCGATCTTTTTTTTTTATCCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTcc                                                    <  2:1189066/137‑1 (MQ=255)
                                                                                                 ccAGATGCGTTTTTCCCGTCCCGTTACTGATATTCAGCGCTCTCCCGATCTCTTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCa                                                <  1:24467/151‑1 (MQ=255)
                                                                                                 ccAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGcc                                                 >  1:116743/1‑150 (MQ=255)
                                                                                                   aGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGt                                              >  1:216436/1‑151 (MQ=255)
                                                                                                    ggtgcctttttCCCTTGCCGGTCCTGATATCCGCCGCTCTCCCGTTTTCTTTTTTTTTCTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAg                                               <  2:494481/142‑1 (MQ=255)
                                                                                                           ttttCCCGTTCCGGTATTGATATCCGGCGCTCTCCCTTTCTTTTTTTTTGTTCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCa                                         <  2:236022/148‑1 (MQ=255)
                                                                                                           tttccccttcccGGTACTAATTTCCGGCGCTCTGCCTTTCCCTTTTTTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCa                                         <  1:960671/138‑1 (MQ=255)
                                                                                                               aCCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGa                                   >  1:1044566/1‑150 (MQ=255)
                                                                                                                 cGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAAc                                >  2:481480/1‑151 (MQ=255)
                                                                                                                   ttCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAAcat                              >  1:1307319/1‑151 (MQ=255)
                                                                                                                   ttCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAAcat                              >  1:530210/1‑151 (MQ=255)
                                                                                                                       ggTACTGATATTCGGCGCTCTGCCGATCTCTTTTTTTGATCCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAAcatc                             <  1:578194/148‑1 (MQ=255)
                                                                                                                           aCTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATcccc                       >  1:654309/1‑151 (MQ=255)
                                                                                                                                       cgcTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACt            >  2:666996/1‑150 (MQ=255)
                                                                                                                                         ctctGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTAtt         >  1:1088998/1‑151 (MQ=255)
                                                                                                                                         ctctGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTAtt         >  2:1412141/1‑151 (MQ=255)
                                                                                                                                              ccGTTCTCTTTTTTTGTTCCCCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCa              <  1:408731/141‑1 (MQ=255)
                                                                                                                                                 aTCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGcc  >  2:810841/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TCCTGCTGGCTGTAAAATTATCCTCTTCAGGAGGAGATATTTAACATCATTGCCGCCTGGGTGCGATTTTTCACTTCCAGACGGCGATACAGGGATTCCAGATGCGCTTTTACCGTTCCGGTACTGATATTCAGCGCTCTGCCGATCTCCTTATTTGATTCGCCCGCCGCTAACATGGTTAAAATCTCCCGCTGGCGGGCGCTTAACGATTTGAGATCTTTAATGTCCTTTTCCGGCGTCGTCCGCCAGTCTCCAGGCAGAAACATCATCCCCATCGCCGCACTATTTACCGCC  >  NC_000913/3671266‑3671559

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