Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,337,7230AC40.8% 112.7 / 19.6 73G566G (GGT→GGGgyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (37/5);  new base C (0/29);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.06e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACG  >  NC_000913/2337575‑2337872
                                                                                                                                                     |                                                                                                                                                     
cgcgcgcgCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAcc                                                                                                                                                     >  2:829086/1‑151 (MQ=255)
 gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCt                                                                                                                                                    <  1:1094413/151‑1 (MQ=255)
 gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCt                                                                                                                                                    >  1:489588/1‑151 (MQ=255)
  cgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCtt                                                                                                                                                   <  1:230525/151‑1 (MQ=255)
    cgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTc                                                                                                                                                 <  1:105618/151‑1 (MQ=255)
                         cAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCg                                                                                                                            >  1:830294/1‑151 (MQ=255)
                           aCTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGt                                                                                                                           >  1:723836/1‑150 (MQ=255)
                            cTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTAt                                                                                                                         >  1:1067283/1‑151 (MQ=255)
                               aTAAACTTTCATCGAATAGACGCGCCCACGGCTGGAGAAGCAAAGAATATGGTCGTGAGTGTTCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATCCGTGCGGCAGTTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                       <  2:793051/150‑1 (MQ=255)
                                tAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                       >  1:477326/1‑149 (MQ=255)
                                tAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGACGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTAtt                                                                                                                        >  1:673796/1‑148 (MQ=255)
                                 aaaCTTTCATCAAATAGACGCGCCCACGGCGGGAGAAGCACAAAATATGTTCGGGATTTTCCCCCCCCAGCAGTCGTTCAAAAAAGTTTTTTTCTTTATTGCGTGCGGCAGTTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAg                                                                                                                     <  1:698044/150‑1 (MQ=255)
                                      tgcgccGAATAGACGCCACGGCGGGTGGAGAAGCACAGAATCTGGGCGGGAGTTTTCGCCACCAGGAGCAGGACGATAAAATTTGCTTTTTTCAAACGTGGGGCAGATTTACCTTCCGCGCCACGACGCTGCGCTTCGTATTCAGaa                                                                                                                   <  1:103547/142‑1 (MQ=255)
                                       tCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCgg                                                                                                              >  2:935328/1‑151 (MQ=255)
                                        cATCAAATAGACGCGCCCACGGCTGGAGAAGCACAAAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTTTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCg                                                                                                               <  2:489609/149‑1 (MQ=255)
                                            gAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGAt                                                                                                         >  2:1007939/1‑151 (MQ=255)
                                             aaTAGACGCGACCCCGGCGGGAGAAGCAAAAAATATGGTCGTGAGTGTTCCCCACCAGCAGTCGGTCAATAAAGTCTTCTTTTTTAATACGTGCGGCAAATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                            <  2:69622/147‑1 (MQ=255)
                                             aaTAGACGCCACCCGGGCGGGAGAAGCACAAAATGGGTTGGTGAGTTTCCGCCCCCAGCATTCGTTCAAAAAATTCTTTTTTTTTAATACGTGCGGCAGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATa                                                                                                        <  1:888375/151‑1 (MQ=255)
                                                aaccgggcccACGGCGGGAAAAGCACAGAATATGGGCTTGAGTGTTCGCCACAAGGAGTCGCACCAAAAGGCCTTCTTTTTTATACGGTGCGGTATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGc                                                                                                                <  1:900622/132‑1 (MQ=255)
                                                aGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACtt                                                                                                     >  1:595638/1‑151 (MQ=255)
                                                aGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACt                                                                                                      >  2:660106/1‑150 (MQ=255)
                                                  ccgcccccACGGTTGAAGAGACAAAGAATGGGTTCGTGAGTGTCCCCCACCAGCAGTGGGTCAATAAAGTTTTTTTTTTAAATGCGTGGGGAATTTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTaa                                                                                                   <  2:360566/145‑1 (MQ=255)
                                                  aCGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGTTAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTa                                                                                                    >  2:518269/1‑150 (MQ=255)
                                                  aCGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTa                                                                                                    >  2:944555/1‑150 (MQ=255)
                                                        ccACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGc                                                                                             >  1:804741/1‑151 (MQ=255)
                                                         cACGGGGGGAGAAGCAAAGAATATGTTCGTGAGTTTCCGCCCCAGGCAGTCGGTGAAAAACGTCTTCTTCTTAAATGCGTGCGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                                                                                  <  1:309610/145‑1 (MQ=255)
                                                         cACGGCTGGAGAAGCAAAAAATATGGTCGGGAGTGTCCCCCACCAGCAGTCGGTCGAAAAAGTCTTCTTCTTTAATCCGGGCGGAAGTTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACg                                                                                                 <  1:118410/146‑1 (MQ=255)
                                                          aCGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGccc                                                                                           >  2:588073/1‑151 (MQ=255)
                                                            ggCGGGAAAAGCAAAAAATATGGTCGTGGGTGTCCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATCCGTGCGGCATTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCt                                                                                          <  1:35151/150‑1 (MQ=255)
                                                             gCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTgg                                                                                        >  2:679943/1‑151 (MQ=255)
                                                              cTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                       >  2:619595/1‑151 (MQ=255)
                                                               tGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                      >  1:1190566/1‑151 (MQ=255)
                                                                ggAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                      >  1:319023/1‑150 (MQ=255)
                                                                ggAAAAGCACAAAATTGGGTCGGGAGTGTTCGCCACCAGCAGTCGGTCAATAAAGTCTTCTTTTTAAATCCGTGCGGCAAATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                      <  2:1022193/150‑1 (MQ=255)
                                                                  agaAGCACAGAATATGGCCGGGATTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATCCGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgaga                                                                                   <  2:866003/151‑1 (MQ=255)
                                                                  agaAGCACAAAATATGGTCGTGATTTTTCGCCCCCAGCAGTCGGTAAATAAAGTCTTCTTTTTTAATCCGTGCGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgaga                                                                                   <  2:240382/151‑1 (MQ=255)
                                                                    acgCACAAAATATGGTCGTGAGTTTTCGCCCCCAGCAGTCGTTCGATAAAGTCTTCTTTTTTAATACGTGGGGAAGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                      <  2:228091/144‑1 (MQ=255)
                                                                       cacaGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCg                                                                              >  1:132428/1‑151 (MQ=255)
                                                                         caGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTc                                                                            >  1:853463/1‑151 (MQ=255)
                                                                          aGAATTGGGTCGTGAGTTTCCGCCCCCAGCGGTGGGTCGATAAAGTCTTTTTTTTTAATGCGGGGGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGc                                                                               <  1:358351/147‑1 (MQ=255)
                                                                            aaTATGGTCGTGTGTGTTCGCCCCCAGCAGTCGGTCTATAAATTCTTCTTTTTTAATACTTGCGGCATATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                           <  2:865378/149‑1 (MQ=255)
                                                                            aaTATGGTCGTGGGTGTCCGCCACCGGCAGTGGGTGGATAAAGTCTTTTTTTTTAATCGGGGGGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag                                                                                <  2:601044/144‑1 (MQ=255)
                                                                                    cGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATc                                                                 >  2:213893/1‑151 (MQ=255)
                                                                                       gtgttttCCCCCCCAGCAGTGGGTCGATAAAGTCTTCTTTTTTAATCCGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTc                                                              <  1:211129/146‑1 (MQ=255)
                                                                                          tttccccccccaGCAGCCGGTCGATAAATTCTTTTTTTTTAAGGCGGGCGGAAGATTTTCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTc                                                              <  1:944501/139‑1 (MQ=255)
                                                                                              ccccccccccattcggtcgATAAAGTCTTCTTTTTTAATACGTGCGGAAGTTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTg                                                           <  2:980254/135‑1 (MQ=255)
                                                                                              cccccccaGCAGTGGGTCAAAAAAGTCTTTTTTTTAAAAAGGGGGGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTg                                                           <  1:513959/142‑1 (MQ=255)
                                                                                              cGCCCCCAGCAGTCGGTCGATAAATTCTTCTTCTTTAATCCGTGCGGAAGATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTcc                                                             <  2:1284358/145‑1 (MQ=255)
                                                                                              cGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                       >  2:1150318/1‑151 (MQ=255)
                                                                                                cccccagaggtcgttcaATAATTTCTTCTTTTTAAATAGGGGCGGAAATTTTTCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGAt                                                     <  1:885977/142‑1 (MQ=255)
                                                                                                         gtcggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATAGTTAACGTAGCCCTGGTGAGAGAGCGTCACTACCACATCTTCCTGGGTGATCAGATCt                                              >  1:433903/1‑149 (MQ=255)
                                                                                                          tcggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                           >  1:1127844/1‑151 (MQ=255)
                                                                                                            ggtcgAAAAAGTTTTCTTCTTTAATCCGTGGGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCtt                                             <  1:1200192/147‑1 (MQ=255)
                                                                                                                gATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCt                                                                >  1:573626/1‑124 (MQ=255)
                                                                                                                 aTAAAGTCTTCTTCTTTAATCCGTGCGGCAGATTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCt                                                                <  2:573647/123‑1 (MQ=255)
                                                                                                                   aaaGTTTTTTTCTTGATTCCGCGCGGCAGTTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt                                      <  1:664219/147‑1 (MQ=255)
                                                                                                                   aaaGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg                                  >  1:658765/1‑151 (MQ=255)
                                                                                                                    aattcttcttctTTAATATGTGCGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGt                                 <  2:429686/148‑1 (MQ=255)
                                                                                                                           cttcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCt                          >  2:847536/1‑151 (MQ=255)
                                                                                                                              cttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGAGAGAGAGCGTCACGACCACATCTTCCTGGGGGATCAGATGTTCCAGGGTGATGTCTGCGCTGt                        >  1:1234194/1‑150 (MQ=255)
                                                                                                                               tttAATCCGGGCGGCAGATTTCCTTCTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGtt                       <  1:1193864/150‑1 (MQ=255)
                                                                                                                               tttAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                      >  2:1106044/1‑151 (MQ=255)
                                                                                                                                ttAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTgg                     >  2:1060166/1‑151 (MQ=255)
                                                                                                                                 tAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGc                    >  2:161789/1‑151 (MQ=255)
                                                                                                                                 tAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGGCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGc                    >  1:749720/1‑151 (MQ=255)
                                                                                                                                  attACGGGCGGAATTTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                      <  1:751817/146‑1 (MQ=255)
                                                                                                                                     aCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCgg                  <  2:459162/149‑1 (MQ=255)
                                                                                                                                     aCGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTg                >  1:96049/1‑151 (MQ=255)
                                                                                                                                         tGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAtt             >  2:988888/1‑151 (MQ=255)
                                                                                                                                            ggCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAtt             <  1:961620/148‑1 (MQ=255)
                                                                                                                                               agttttCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAg         <  2:922213/146‑1 (MQ=255)
                                                                                                                                                gATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTac      >  1:403154/1‑151 (MQ=255)
                                                                                                                                                    tACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacgacg  >  1:559584/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACG  >  NC_000913/2337575‑2337872

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