Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,337,7230AC42.4% 107.5 / 23.8 66G566G (GGT→GGG) gyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (34/4);  new base C (0/28);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.05e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTAC  >  NC_000913/2337574‑2337868
                                                                                                                                                     |                                                                                                                                                  
ccGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAc                                                                                                                                                   <  2:501232/151‑1 (MQ=255)
       cgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTccc                                                                                                                                           >  1:534661/1‑151 (MQ=255)
          cACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGcc                                                                                                                                        >  2:117534/1‑151 (MQ=255)
          cACGAGAGGCTTCCGGCAACTGATAAACTTTCCTCGAATAGCCGCGACCACGGCTGGAGAAGACCAGAATATGGTCGTGAGTGTTCGCCAGCAGCAGTCGGTCGATAAAGCCTTCTTCTTTAATACGTGCGGCAGAATTAC‑CTTTCCCGc                                                                                                                                         >  1:476068/1‑150 (MQ=255)
            cGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCac                                                                                                                                      >  2:691886/1‑151 (MQ=255)
             gAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCacg                                                                                                                                     >  2:339760/1‑151 (MQ=255)
                  gCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCt                                                                                                                                >  2:1132782/1‑151 (MQ=255)
                  gCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCt                                                                                                                                >  2:9197/1‑151 (MQ=255)
                   cTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTg                                                                                                                               >  1:286730/1‑151 (MQ=255)
                    ttCCAGAAAGTGATAAACTTTCATCGAATAGACGCGACCACGCCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACAAGTAGCCGGTCGATAAAGTCTTTTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCTCCACGACGCt                                                                                                                                >  1:924119/1‑149 (MQ=255)
                           aaCTAATAACCTTTCATCAAATAGACGCGCCCACGGCTGGAGAAGCACAGAATAGGGTCGTGAGTTTTCCCCCCCAGCAGTCGGTCAATAAAGTCTTTTTCTTTAATACGTGCGGCAGTTTTCC‑CTTTCCCGCCACGACGCTGCGCTTCg                                                                                                                        <  2:784761/150‑1 (MQ=255)
                            aCTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTa                                                                                                                      <  2:431386/151‑1 (MQ=255)
                                aaaaaCTTTCATCAAATAGACGCGACCACGGCTGGAGAAGCAAAGAATATGGTCGTGAGTGTCCGCCACCAGCAGTCGGTCGATAAATTCTTCTTTTTTAATCCGGGCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                   <  2:285244/148‑1 (MQ=255)
                                  aaaCTTTCATCAAATAGACGCGCCCCGGGTGGGAGAAGCACAGAATAGGGTCGGGAGTGTTCGCCACCAGCAGTCGGTCGAAAAAGTCTTCTTCTTTAATCCGTGCGGAAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGa                                                                                                                <  1:548743/151‑1 (MQ=255)
                                      tttCATCAAAAAGCCGCGCCACCGGCTGGAAAAGCAAGGAATATGGTCGTGAGTGTCCGCCACAAGCAGGCGGTCGATAAAGTCTTCTTCTTAAATACGTGCGGCAGTTTTCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAg                                                                                                             <  1:682348/150‑1 (MQ=255)
                                         cATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGc                                                                                                         >  2:21158/1‑151 (MQ=255)
                                          aTCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                        >  1:754854/1‑151 (MQ=255)
                                            cGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTg                                                                                                       >  2:798776/1‑150 (MQ=255)
                                            cGAATAGACGCGACCACGGCGGGAGAAGCACAGAATATGGTCGTGAGTGTCCGCCACCAGCGGTCGGTCGATAAATTCTTCTTTTTTAATAGGTGCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                        <  1:745667/149‑1 (MQ=255)
                                              aaTAGACGCGACCACGGCTGGAGAAGCACAGAATAGGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCC‑TTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTg                                                                                                       <  2:341168/148‑1 (MQ=255)
                                                  gACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACtt                                                                                                 >  1:464155/1‑150 (MQ=255)
                                                       gACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                            >  2:767514/1‑150 (MQ=255)
                                                         ccccggcggGAAAAGCAAAAAATATGGTCGGGAGTGTTCGCCCCCAGCAGTCGGTGAAAAAAGTCTTCTTCTTTAATCCGTGCGGCAGATTTCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                            <  2:920117/145‑1 (MQ=255)
                                                         ccACGGCTGGAGAAGCACAAAATATGGTCGTGAGTTTTCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGc                                                                                         <  2:26792/151‑1 (MQ=255)
                                                             ggCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTg                                                                                     >  2:977404/1‑151 (MQ=255)
                                                                tGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                  >  2:1064568/1‑151 (MQ=255)
                                                                 ggAGAAGCACAGAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGGGCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                  <  2:724875/150‑1 (MQ=255)
                                                                 ggAAAAGCACAGAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCAAAAAAGTCTTTTTTTTTAATCCGTGCGGCAGATTTCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                   <  2:478025/149‑1 (MQ=255)
                                                                                  ggTCGTGAGTGTTCGCCCCCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGAc                                                                   <  2:562085/148‑1 (MQ=255)
                                                                                  ggTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcac                                                                >  2:39714/1‑151 (MQ=255)
                                                                                        gAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTc                                                          >  2:706405/1‑151 (MQ=255)
                                                                                          gggtTCGCCCCCAGCATTCGTTCGTAAAAGTCTTTTTCTTTAATCCGGGCGGCAGTTTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGTGAGAGCGTCACGACCACATc                                                             <  1:949390/144‑1 (MQ=255)
                                                                                           ttttCGCCCCCACCGGTCGGTCGAAAAAGTCTTTTTCTTAAATCCGGGGGGCAGTTTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt                                                        <  1:234865/148‑1 (MQ=255)
                                                                                           tgtcccccaccaGCAGTCGGTCGATAAAGTCTTTTTCTTTAATACGTGAGGCATATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCtt                                                           <  1:1050375/141‑1 (MQ=255)
                                                                                             ttCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                     >  1:253763/1‑151 (MQ=255)
                                                                                              tCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGt                                                    >  1:106722/1‑151 (MQ=255)
                                                                                              tCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGt                                                    >  1:788123/1‑151 (MQ=255)
                                                                                               cGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGTGCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                     <  2:260795/149‑1 (MQ=255)
                                                                                               cGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCCCGCCCCCCTCTTCCTGGGTg                                                   >  1:886317/1‑151 (MQ=255)
                                                                                               cGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                   >  1:687541/1‑151 (MQ=255)
                                                                                               cGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTAAATGCGTGCGGCAGATTTCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                   <  1:521056/151‑1 (MQ=255)
                                                                                                  caccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                >  2:830018/1‑151 (MQ=255)
                                                                                                    ccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCa                                               <  2:613660/150‑1 (MQ=255)
                                                                                                        cggtcggtcgATAAAGTTTTCTTTTTTAATAGGTGCGGAATATTTCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTgg                                                      <  1:117386/137‑1 (MQ=255)
                                                                                                          gtcgggcgATAAAGTCTTTTTTTTTAATGCGGGGGGCAGATTTCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATc                                           <  2:678962/148‑1 (MQ=255)
                                                                                                           tcggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                       >  2:773682/1‑151 (MQ=255)
                                                                                                              gtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                    >  2:1146241/1‑151 (MQ=255)
                                                                                                                  aTAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCAACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTg                                 >  1:260935/1‑150 (MQ=255)
                                                                                                                    aaaGTCTTCTTTTTTAATCCGTGCGGCAGATTTCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGa                                <  2:108904/149‑1 (MQ=255)
                                                                                                                    aaaGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg                              >  1:169218/1‑151 (MQ=255)
                                                                                                                     aaGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGt                             >  2:1155052/1‑151 (MQ=255)
                                                                                                                       gttttcttctaTATGTCGTGCGG‑AAATTTACTCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGAt                               <  1:691558/144‑1 (MQ=255)
                                                                                                                         cttcttcttTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgc                         >  1:990999/1‑151 (MQ=255)
                                                                                                                           tctttttTAATCCGTGCGGCAGTTTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCt                           <  1:310626/147‑1 (MQ=255)
                                                                                                                                tttAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                  >  2:996600/1‑151 (MQ=255)
                                                                                                                                 ttAATAGGGGCGGCAGATTTCC‑TTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGt                    <  1:781548/148‑1 (MQ=255)
                                                                                                                                  tAATACGTGCGGCAGATTTCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                  <  2:788465/149‑1 (MQ=255)
                                                                                                                                    aTACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCgg              >  2:1114203/1‑151 (MQ=255)
                                                                                                                                     tACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCgg              <  2:754016/150‑1 (MQ=255)
                                                                                                                                        gTGCGGCAGATTTCC‑TTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGt             <  2:714652/148‑1 (MQ=255)
                                                                                                                                        gTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAt          >  2:910253/1‑151 (MQ=255)
                                                                                                                                          gCGGTAAATTTAC‑TTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAttt        <  1:414992/151‑1 (MQ=255)
                                                                                                                                          gCGGCAGATTTCC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAt          <  1:184829/149‑1 (MQ=255)
                                                                                                                                             gCAGATTTCC‑CCTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAg     <  2:279400/151‑1 (MQ=255)
                                                                                                                                                gATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTac  >  2:188169/1‑151 (MQ=255)
                                                                                                                                                  tttCC‑CTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTa   <  1:1064151/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAC‑CTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTAC  >  NC_000913/2337574‑2337868

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