Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,515,5650AG59.7% 38.6 / 99.6 77Y117H (TAC→CAC) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (14/17);  new base G (21/25);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCA  >  NC_000913/3515419‑3515711
                                                                                                                                                  |                                                                                                                                                  
gTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAACTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCg                                                                                                                                                 >  2:547764/1‑150 (MQ=255)
  aaGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGtt                                                                                                                                              >  1:482984/1‑151 (MQ=255)
   aGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGtt                                                                                                                                              >  1:509648/1‑150 (MQ=255)
     aCTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCaa                                                                                                                                           <  2:35216/151‑1 (MQ=255)
        cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAat                                                                                                                                         >  2:607534/1‑150 (MQ=255)
        cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAat                                                                                                                                         <  2:193058/150‑1 (MQ=255)
        cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAata                                                                                                                                        <  1:98353/151‑1 (MQ=255)
           aaTAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaa                                                                                                                                     <  1:190018/151‑1 (MQ=255)
              aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaaaa                                                                                                                                   <  2:362142/150‑1 (MQ=255)
              aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaa                                                                                                                                   <  1:436382/150‑1 (MQ=255)
               gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAc                                                                                                                                 <  1:10190/151‑1 (MQ=255)
               gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAAc                                                                                                                                 >  1:242194/1‑151 (MQ=255)
               gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAAc                                                                                                                                 >  1:289086/1‑151 (MQ=255)
               gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAAc                                                                                                                                 >  2:570550/1‑151 (MQ=255)
                   aaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGt                                                                                                                              >  1:360909/1‑150 (MQ=255)
                    aGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTa                                                                                                                             >  2:303432/1‑150 (MQ=255)
                          tCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGccc                                                                                                                       >  1:79360/1‑150 (MQ=255)
                          tCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGccc                                                                                                                       <  1:133661/150‑1 (MQ=255)
                             gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGac                                                                                                                    <  2:382343/150‑1 (MQ=255)
                             gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCg                                                                                                                      <  1:403039/148‑1 (MQ=255)
                             gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacg                                                                                                                   <  2:229264/151‑1 (MQ=255)
                              cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacgg                                                                                                                  >  2:306737/1‑151 (MQ=255)
                                  ttttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaac                                                                                                               <  2:423224/150‑1 (MQ=255)
                                  ttttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaac                                                                                                               <  1:512023/150‑1 (MQ=255)
                                    ttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaacgg                                                                                                             <  2:274423/150‑1 (MQ=255)
                                       aGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTAcccccc                                                                                                                       >  1:347963/1‑137 (MQ=255)
                                       aGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGccc                                                                                                                       <  2:347956/137‑1 (MQ=255)
                                        gCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                                                                                        >  2:358557/1‑151 (MQ=255)
                                               gATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCACCCGACGGAACGGATCCTGGCtttt                                                                                                 >  1:286706/1‑151 (MQ=255)
                                                   cAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgtt                                                                                              >  1:89209/1‑150 (MQ=255)
                                                    aaCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTTAAATATAAAAACAGAACCCCCCGACGGAAAGGACCCTGGCtt                                                                                                   >  1:435713/1‑144 (MQ=255)
                                                    aaCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGa                                                                                            <  1:183257/151‑1 (MQ=255)
                                                     aCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGAACGGATCCTGGCTTttgt                                                                                               >  1:162491/1‑147 (MQ=255)
                                                          gCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAActc                                                                                        <  2:79641/149‑1 (MQ=255)
                                                              ccGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcgc                                                                                   <  1:250599/150‑1 (MQ=255)
                                                                gggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcgcg                                                                                  <  2:400185/149‑1 (MQ=255)
                                                                 ggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGaa                                                                               <  2:351283/151‑1 (MQ=255)
                                                                  gAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAAc                                                                              <  2:226998/151‑1 (MQ=255)
                                                                   aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACt                                                                             <  1:172981/151‑1 (MQ=255)
                                                                    aaTAGGGTTTTTTGTAGCGGCCGAACGGCACTTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg                                                                            <  2:195623/151‑1 (MQ=255)
                                                                       aGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                          <  1:3341/150‑1 (MQ=255)
                                                                       aGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                          <  1:260552/150‑1 (MQ=255)
                                                                       aGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                          >  2:618075/1‑150 (MQ=255)
                                                                              tttGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTaaaac                                                                  >  1:19726/1‑149 (MQ=255)
                                                                                   aGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGgc                                                             <  2:454371/151‑1 (MQ=255)
                                                                                         cGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAtt                                                       <  1:249976/151‑1 (MQ=255)
                                                                                         cGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGAACGGATCCTGGGTTTTGTTTAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAtt                                                       >  2:105653/1‑151 (MQ=255)
                                                                                         cGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAtt                                                       >  2:156802/1‑151 (MQ=255)
                                                                                                cACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAAGAGTACCCACCGACGGAACGGACCCTGGCGTTTTGTGGACTCTTCATGGAACTGACAGTAAACCTCGGCGCAATTTCttt                                                  >  2:469967/1‑149 (MQ=255)
                                                                                                 aCGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCgg                                               >  2:528984/1‑151 (MQ=255)
                                                                                                 aCGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCgg                                               <  2:458014/151‑1 (MQ=255)
                                                                                                   gTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAAATGATAGTAAACCTCGGCGCAATTTGTTTCggg                                              >  2:205673/1‑150 (MQ=255)
                                                                                                      aaaCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaaca                                          <  2:437762/151‑1 (MQ=255)
                                                                                                      aaaCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaac                                           <  2:94349/150‑1 (MQ=255)
                                                                                                           cACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaacaaac                                       <  2:20615/149‑1 (MQ=255)
                                                                                                           cACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAg                                     <  2:593907/151‑1 (MQ=255)
                                                                                                             ccGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCt                                   <  1:146116/151‑1 (MQ=255)
                                                                                                              cgcgCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTc                                  <  2:312105/151‑1 (MQ=255)
                                                                                                                 gcAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGTAATTTGTTTCGGGAACAAACAGCTcgcg                               <  1:177027/151‑1 (MQ=255)
                                                                                                                  cAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGt                              <  1:81595/151‑1 (MQ=255)
                                                                                                                   aGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTg                             <  2:276355/151‑1 (MQ=255)
                                                                                                                     aTTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAATCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCg                           >  2:10331/1‑151 (MQ=255)
                                                                                                                            cGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCGTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCCGCGCAATTTTTTTCGGGAACAAACAGCTCGCGTGCGGCCTgt                     >  1:503084/1‑150 (MQ=255)
                                                                                                                                acaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgta                <  1:569931/151‑1 (MQ=255)
                                                                                                                                acaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAATCTTCGCGGAACTGATAGTAAAACTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacg                  >  2:14360/1‑149 (MQ=255)
                                                                                                                                 caGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgtac               >  1:547211/1‑151 (MQ=255)
                                                                                                                                  aGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACt              <  1:38983/151‑1 (MQ=255)
                                                                                                                                  aGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACt              <  1:112298/151‑1 (MQ=255)
                                                                                                                                   ggCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtc             >  1:6424/1‑151 (MQ=255)
                                                                                                                                   ggCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtc             >  1:254198/1‑151 (MQ=255)
                                                                                                                                    gCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtc             >  2:337384/1‑150 (MQ=255)
                                                                                                                                      cGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc          >  2:259367/1‑151 (MQ=255)
                                                                                                                                      cGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc          >  2:596759/1‑151 (MQ=255)
                                                                                                                                      cGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc          >  1:271077/1‑151 (MQ=255)
                                                                                                                                       gTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAAACTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc          >  2:348317/1‑150 (MQ=255)
                                                                                                                                            aaCCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTAc     <  2:579658/150‑1 (MQ=255)
                                                                                                                                              ccGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCa  <  2:571053/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCA  >  NC_000913/3515419‑3515711

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