Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,7330TG41.9% 122.0 / 17.5 75G20G (GGT→GGGyegDHsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base T (11/32);  new base G (31/0);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 7.06e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCG  >  NC_000913/2147584‑2147883
                                                                                                                                                      |                                                                                                                                                        
aCGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTa                                                                                                                                                         <  1:982469/151‑1 (MQ=255)
  gCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTaaa                                                                                                                                                       <  1:352532/151‑1 (MQ=255)
  gCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTaaa                                                                                                                                                       <  1:888259/151‑1 (MQ=255)
    tGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGaaaac                                                                                                                                                     >  1:10373/1‑149 (MQ=255)
            gctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGc                                                                                                                                             <  1:811432/151‑1 (MQ=255)
                tCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGTGACGGGAAAACGCATTTTCtta                                                                                                                                          >  1:1004973/1‑148 (MQ=255)
                tCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGTCGGTAAACCGCATTTTCt                                                                                                                                            >  2:475641/1‑148 (MQ=255)
                 cATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTaaaa                                                                                                                                        <  1:213660/151‑1 (MQ=255)
                    atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAAt                                                                                                                                      <  1:642661/150‑1 (MQ=255)
                      tCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGGGTGGCGGGAAAACGCACTTGCTTAAAAttg                                                                                                                                    >  2:535192/1‑148 (MQ=255)
                      tCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATgg                                                                                                                                    <  2:823821/150‑1 (MQ=255)
                       ccGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCTGCGGGGCGGGAAAAAGCATTTTTTaaaaa                                                                                                                                       >  1:694647/1‑146 (MQ=255)
                          gTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCTTCGTGGCGGGTAAACGCCCTTTTTTAAAATGGgaaa                                                                                                                                >  2:782952/1‑150 (MQ=255)
                            tGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGTAAAACCCATTTTTTAAAAATAGaaaa                                                                                                                                >  1:753078/1‑148 (MQ=255)
                               cgcgACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAg                                                                                                                          <  1:594873/151‑1 (MQ=255)
                                gcgACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGc                                                                                                                         <  1:861835/151‑1 (MQ=255)
                                  gACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGTGAAGGGAAAACGCCTTTTCTAAAAAAGGAAAACGAAAGCa                                                                                                                        >  2:906673/1‑150 (MQ=255)
                                   aCGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGTGAGGGGAAAACGCATTTGTTTAAAAATGAAAACg                                                                                                                              >  2:90571/1‑143 (MQ=255)
                                       cAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGGGTGGCGGTAAAACGCACTTTTTAAAAATGGAAAAAAACAACACgc                                                                                                                     >  1:542761/1‑148 (MQ=255)
                                       cAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACgctgc                                                                                                                  <  2:52605/151‑1 (MQ=255)
                                       cAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGGAAACCGCATTTTCTAAAAAAGGAAAACGACAGCACgctgc                                                                                                                  >  1:847290/1‑151 (MQ=255)
                                         aCATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCAGCGTGTGGGGTAAACGCATTTGCTTAAAAAAGGAAACGACaaccaggtgctg                                                                                                                >  2:638473/1‑140 (MQ=255)
                                           aTTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCCTTTGCTAAAAAAGGAAAAAGACAGCACGCCGCTGcc                                                                                                              >  1:324887/1‑151 (MQ=255)
                                             ttCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCt                                                                                                             <  1:1050115/150‑1 (MQ=255)
                                               cGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCa                                                                                                          <  1:559863/151‑1 (MQ=255)
                                                gCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCTTCGGGAGGGGTAAACCGATTTGCTAAAAAAGGGAAAAGa                                                                                                                             >  2:151789/1‑131 (MQ=255)
                                                gCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCaa                                                                                                         <  1:101202/151‑1 (MQ=255)
                                                 cTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGTGGGGGGAAAAAGCCTTTTCTAAAAAAGGAAAAAAACAGAACGCTCCTGCCTTTaa                                                                                                         >  2:610778/1‑150 (MQ=255)
                                                  tAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAACGACAGCACGCCGCTGCCTTCaa                                                                                                         >  1:603363/1‑149 (MQ=255)
                                                    aaGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCTGCGTGGCGGTAAAACGCCTTTTATCAAAATTGAAAAACACAGGACGCCGCGGcccttaa                                                                                                         >  1:449659/1‑142 (MQ=255)
                                                     aGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGTGAAACGCCTTTTCTTAAAATGGAAAAAGACAACACGCTGCCGCCTTCCATGCtt                                                                                                   >  2:1057790/1‑151 (MQ=255)
                                                       tCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGACCTGGGGTAAGGGTAACCCCATTTTTTAAAAAATAAAAAAGACACCACGCTGATGCCTTCAATGCttt                                                                                                  >  2:973731/1‑150 (MQ=255)
                                                        cACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTg                                                                                                 <  1:929908/150‑1 (MQ=255)
                                                             cccTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTATGCGGGAGGGGAAACCGCAATTTCTAAAAATGGAAAACGATAGCACGCTGCTGCCTTCAATGCTTTGCGCGcc                                                                                           >  1:159338/1‑151 (MQ=255)
                                                                ttcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGGGGGACGGGAAAACGCATTTTCTTAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCaa                                                                                         >  1:18177/1‑150 (MQ=255)
                                                                    tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAAcg                                                                                       <  2:226791/148‑1 (MQ=255)
                                                                     cACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCTGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAAAGACAGCACCCGGCTGCCTTCTATGCTTTTCGCGCCAAcgcggg                                                                                   >  2:238510/1‑149 (MQ=255)
                                                                       ccGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTCTAAAAATGGAAAACAACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCCAcgcgcg                                                                                   >  2:159550/1‑147 (MQ=255)
                                                                          gCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCTTGTGTGGAGGGAAACCGCCCTTGCTACAAAAGGAAAACAACAGCACGCTGCCGCCCTCAATTCTTTTCGCGCCGACGCGGGAAg                                                                                >  2:597072/1‑149 (MQ=255)
                                                                            aTGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTTCTAAAAAAGGAAAAAGACAGCACGCCGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCgg                                                                             >  1:102520/1‑150 (MQ=255)
                                                                               gggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGAGGGGAAAACGCATTTTTTTAAAATGGGAAACGACAGCACGCCGCTGCCTTCAATGTTTTGTGCGCCAACGCGTGAAGCgg                                                                             >  1:788832/1‑147 (MQ=255)
                                                                                ggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGc                                                                        <  2:159474/151‑1 (MQ=255)
                                                                                    tattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGaa                                                                    <  1:535448/150‑1 (MQ=255)
                                                                                     attTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAAt                                                                   <  2:431961/150‑1 (MQ=255)
                                                                                      tttCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTTGTGGGGGGTAACCGCCTTTGCTAAAAATAGAAAAAGACAGCACGCTGCTGGCCTCAATGCTTTGCGCGCCCACGCGTGAAGCGGGCAGCGAATgg                                                                 >  2:586954/1‑151 (MQ=255)
                                                                                      tttCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGGGACGGGAAACCGCATTTGCTTAAAATGGAAAACGAAAGCACGCTGCCGCCCTCAATGCTTTGTGCGCCAACGCGGGAAAGGGGAAGCGAATgg                                                                 >  2:428304/1‑151 (MQ=255)
                                                                                             tttATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCATTTTCTTAAAAAGGAAAACGACAGCACCCCCCTGCCCTCAATGCTTTGCGCGCCAACGCGGGAAGCGGGAAGCGAATGGCTGTAcc                                                          >  1:171861/1‑151 (MQ=255)
                                                                                              ttATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAACACGCCGCTGCCCTCAATGCTTTTCGCGCCAAAGCGGGAAGCGGGAAGCGAATGGCTGTAcc                                                          >  1:680099/1‑150 (MQ=255)
                                                                                                aTTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGc                                                        <  1:475886/150‑1 (MQ=255)
                                                                                                 ttGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCATTTGTGGCGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAGCACGCCGCTTCCTTCAATTCTTTGCGCGCCAAAGCGTTAAACGGGAAAGGGATGGGTGTACcccc                                                       >  2:743114/1‑150 (MQ=255)
                                                                                                  tGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCca                                                      <  2:1044437/150‑1 (MQ=255)
                                                                                                        tGATTACGGTACAGCAAACTGTTCAGTGGCGGTTATGCGCGACGGGAAACCCCATTTTCTTCAAATGGGAAACGACAGCCCGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGCGAAACGGGAAGCGACAGGCCGTACCGCCATCctggt                                               >  2:608128/1‑149 (MQ=255)
                                                                                                          aTTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGGCGGTAAAACGCCTTTGCTAAAAAAGGAAAACGAAAGCACGGCGCTGCCTTCAAAGGTTTGCGCGCCAACGGGGGAAAGGGGAAGCGGATGGGTGTTCCGCcctcc                                                   >  2:825207/1‑144 (MQ=255)
                                                                                                          aTTACGGTACAGCAAACTGTTCAGTGGCGGTACAGACTGACGGAGAACCCCATTTGCTAAAAAAGGAAAAAAACAGCACCCTGCCGGCTTTAATAGCTTTGGCGCCAACGAGGGGAACGGAGAGGGAAGGGCTGTACCGCcatcat                                                  >  1:994772/1‑146 (MQ=255)
                                                                                                           ttACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTGCTAAAAATGGAAAAGAAAACAACCCGCCTCCCTTCAAGCTTTTCGCGGCCAAGCGGGGAACGGGAAACGGAATGGCGGTACCCCcaccttatgtt                                             >  1:169514/1‑141 (MQ=255)
                                                                                                               ggTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCgg                                        <  1:498313/151‑1 (MQ=255)
                                                                                                                   cAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTGCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGGAAGCGAATGGCTGTACCGCCCTCATGATGTTTCgg                                        >  1:562806/1‑147 (MQ=255)
                                                                                                                   cAGCAAACTGTTCAGTGGCGGTCACGCGTGGGGGGAAAACGCACTTTCTAAAAATGGGAAACGACAGCACGCCGCTGCCTTCCATGCTTTGCGCGCCAACGCGTGAAGCGGGAAGGGAATGGCTGTACCGCCATCATGATGTTCCgg                                        >  2:13700/1‑147 (MQ=255)
                                                                                                                     gCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAATGGAAAACGAAAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAACGCGGGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGtttc                                          >  1:257898/1‑141 (MQ=255)
                                                                                                                      cAAACTGTTCAGTGGCGGTCATTCGTGACGGGAAACCCCCTTTGCTAAAAAATGAAAACGCCAGCACCCCGCCGCCCTCAATGCCTTTCGCGCCCACGCGCGAAACGGGAAGCGGAAGGCTGTACCCCCCTCCTGATGTTcc                                          >  2:551782/1‑142 (MQ=255)
                                                                                                                      cAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAgacga                                 <  1:272772/151‑1 (MQ=255)
                                                                                                                       aaaCTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAgacgac                                <  1:973165/151‑1 (MQ=255)
                                                                                                                        aaCTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTCTAAAAATAGAAAAAGGAAGCACGCTGCTGGCCTCCAAGCTTTGCGCGCCAACCAGGGAAGCGGGAAACGGAAGGGTGTACCGCCCTCATGATGTTTCCGCAga                                    >  2:199769/1‑146 (MQ=255)
                                                                                                                            gTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGAt                             <  2:861811/149‑1 (MQ=255)
                                                                                                                            gTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGa                           <  1:280740/151‑1 (MQ=255)
                                                                                                                             ttCAGTGGCGGGCATGCGTGGCGGGAAAACGCCTTTTCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAAACCGTGAAGCGGGAAGGGAATGGCTGTACCCCCACCATGATGTTTCgg                                        >  2:631233/1‑137 (MQ=255)
                                                                                                                                  tcgCGGTCATGCTTGACCGTAAACCGCATTTGCTAAAAATGGAAAACCACAGCACGGTGCTCGCTTCAATGCTTTGCGCTCCCCCGCGTGAAGCGGTAAGCGAATGGCTGCACCGCCATCATGATGTGCCGGCCGACGACGATGAAAc                        <  2:135794/146‑1 (MQ=255)
                                                                                                                                      ggTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGgcg                 >  2:314410/1‑151 (MQ=255)
                                                                                                                                       gTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGgcgc                <  2:1019947/151‑1 (MQ=255)
                                                                                                                                         cATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTg              <  1:819899/151‑1 (MQ=255)
                                                                                                                                         cATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTg              <  1:715973/151‑1 (MQ=255)
                                                                                                                                          aTGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGt             >  2:388594/1‑151 (MQ=255)
                                                                                                                                                gACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTc       <  1:683571/151‑1 (MQ=255)
                                                                                                                                                 aCGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCg      <  2:760277/151‑1 (MQ=255)
                                                                                                                                                    ggTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCggg    <  1:485477/151‑1 (MQ=255)
                                                                                                                                                      tAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCg  >  2:953753/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
ACGCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCG  >  NC_000913/2147584‑2147883

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: