Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA40.0% 139.7 / 40.7 90E384D (GAA→GATyhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (12/42);  new base A (36/0);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 4.10e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACG  >  NC_000913/3508518‑3508812
                                                                                                                                                 |                                                                                                                                                         
tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATtaca                                                                                                                                                     >  2:780530/1‑150 (MQ=255)
  tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATcacaca                                                                                                                                                   >  1:95155/1‑150 (MQ=255)
  tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGGGGAAAAAATcac                                                                                                                                                       >  2:854359/1‑147 (MQ=255)
   ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCtctcac                                                                                                                                                  <  2:809922/150‑1 (MQ=255)
   ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATtacaca                                                                                                                                                   >  1:355248/1‑149 (MQ=255)
       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTa                                                                                                                                              <  1:1079462/150‑1 (MQ=255)
       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTa                                                                                                                                              <  2:1256283/150‑1 (MQ=255)
       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATAACACACCGt                                                                                                                                               >  2:298727/1‑149 (MQ=255)
         aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGcc                                                                                                                                           <  1:149180/151‑1 (MQ=255)
          gCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGccc                                                                                                                                          <  1:136379/151‑1 (MQ=255)
           ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTaacaacaggagtccc                                                                                                                                          >  2:121251/1‑140 (MQ=255)
           ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGCGGGAAAAAACACAAACAGTCGgccc                                                                                                                                         >  2:1031769/1‑151 (MQ=255)
             ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCca                                                                                                                                       <  1:448068/151‑1 (MQ=255)
             ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGCGGCAAAAATCATACACgggataacgc                                                                                                                                        >  1:628775/1‑142 (MQ=255)
               agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCca                                                                                                                                       <  1:820456/149‑1 (MQ=255)
                 agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGGAGccccccc                                                                                                                                      >  1:467506/1‑146 (MQ=255)
                  gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCCGAAAAAATTACACCCCGGAGCCCCcac                                                                                                                                      >  1:1108984/1‑147 (MQ=255)
                     cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGGAACCCCCCCAACa                                                                                                                                  >  1:82960/1‑148 (MQ=255)
                     aaaTATGAAGATGCTCACGGGCGAGGGAGTAAGCCGTCGGATCAAAACTCATAGGCTGTCTCCTGCAAACCGGAGGGAGGAAAGACAAGAGCGACCATTTTGTGGCGGTGGGGCGCCGAAAAATTtacacaaaataccccc                                                                                                                                         >  2:958150/2‑131 (MQ=255)
                          tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGTAGgccccccccccgaaagacg                                                                                                                          >  1:1205403/1‑137 (MQ=255)
                           gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAACACCGGAGCCCCCACATCACAAAACggg                                                                                                                        >  1:328438/1‑149 (MQ=255)
                           gCAGATGCTCATGGGCGAGGCTGTATACCGTCGCATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCACTTTGCGGCGGTGGCGCGCGGAAAAATTCACACAGGGTAGCCCCCACATCAGAAAGCgg                                                                                                                         <  2:444751/150‑2 (MQ=255)
                            cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCCCGGGAAAATTTACCCACCGGAGGCCCCCCAACAGGAAGCgtgg                                                                                                                       >  1:41968/1‑151 (MQ=255)
                               aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATTAAAAACCGGAGCCCCCACACCAAAAAACgggg                                                                                                                       >  2:1056320/1‑148 (MQ=255)
                                tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTAAAAAACGTAGCCCCCACCTCAGAAAACggggtggg                                                                                                                   >  2:198311/1‑149 (MQ=255)
                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATTAAACACAGTAGCCCCCCCATCAGAAAAGGTGGTGCGCACCCCGAGc                                                                                                         >  1:632141/1‑150 (MQ=255)
                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGACCCCCCAACACAAAACGGGGGGGGcac                                                                                                                >  1:868421/1‑143 (MQ=255)
                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGTAGCCCCCCCCTCAGAAAGCGGGGGGGGCACAcc                                                                                                             >  1:555901/1‑146 (MQ=255)
                                           gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAAAACGGGAGGCCCCCCCACACAAAAAGTGGGGGGGACACCCAGc                                                                                                         >  2:531267/1‑150 (MQ=255)
                                            aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGAAAAAATAACAACCCGGGCGCCCCACAACACAAAGCGTGGTGTGCGCAcc                                                                                                             >  2:250340/1‑145 (MQ=255)
                                                    aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGGGGGGGAAAAATCACACCCGGGAGGCCCCCCATCACAAAAACCGGGGCGGACACCGAGCCGCCGAGCg                                                                                                >  2:381131/1‑150 (MQ=255)
                                                     cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                              <  1:217856/151‑1 (MQ=255)
                                                     cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                              <  1:1068594/151‑1 (MQ=255)
                                                          cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGTAGCCCCCCCAACAGAAAACGTGGTGGGCACACCGAGCCGCGTAGCGATCATcc                                                                                        >  1:130672/1‑151 (MQ=255)
                                                               cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                                   <  2:144316/151‑1 (MQ=255)
                                                                aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCCTccc                                                                                       >  1:308019/1‑146 (MQ=255)
                                                                 aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGTAGCCCCCCCACCAGAAAACGGGGTGGGCACACCCAGCCGCGTAGCGAACATcccct                                                                                     >  1:840392/1‑145 (MQ=255)
                                                                 aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                                 <  2:869596/151‑1 (MQ=255)
                                                                  aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATCACACAAAGGAACCCCCCCCTCCGAAAACGAGGTGGGGACACCGAGCTGCGTAGCGAGCAACCCTTGc                                                                                   >  2:1276045/1‑148 (MQ=255)
                                                                      cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                                                                             <  2:996842/150‑1 (MQ=255)
                                                                      cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                                                                             <  1:1227065/150‑1 (MQ=255)
                                                                      cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGGCCGGAAAAATCACACACGGGAGGCCCCCCCCCAGAAAAAGGGGGGCGCACACCGAAGGGCGGAGGGGGCATCCCCTGCTGGTaa                                                                             >  1:47200/1‑150 (MQ=255)
                                                                        tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGGGAAAATCACACACCGGAGCCCCCCCATCAGAAAACGGGGGGGGCACACCGAGCTGCGGAGCGCTCATCCCCTGCTGATAAAcc                                                                          >  1:921133/1‑151 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAGCGGGGGGCGCCCACCCAGCTGCGTAGCGGACCTCCCCTGCTGATAAAcc                                                                          >  1:238447/1‑150 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCCGCGGAGCGCTCATCCCCTTCTGATaaaa                                                                           >  1:140052/1‑148 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGCAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                         <  2:410452/151‑1 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                         <  1:239086/151‑1 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGCGGGGAAAAAATCACACACGGGAGCCCCCCCATACGAGAAGAGGGGGGGGACACCGCGCTGCGGAGGGGTCATCCCGTGCTGaaaa                                                                             >  2:506447/1‑147 (MQ=255)
                                                                          gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAAACGTAGCCCCCACAACAGAAAAGGGGGGGGGCCCACCCAGCTGCGGAGCGATCATCCCCTGCTGATaaa                                                                            >  1:838256/1‑147 (MQ=255)
                                                                           ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAAAACGGTAGCCCCCCCCTCAGAAAACGGGGGGGGGACACCCAGCTGCGTAGCGATCATCCCCTTCTGTTAAACCa                                                                         >  1:1220896/1‑149 (MQ=255)
                                                                              gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                                     <  1:854390/150‑1 (MQ=255)
                                                                              gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGGGGGAAAATTTTGACAAGGGAACCCCCACCTCACAAAAAGCGGTGCGGAGCCCGAGCAGCGGGGCCAACATCCCTTGCTTGCAAACCAGc                                                                       >  1:109397/1‑148 (MQ=255)
                                                                                ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATAACACACAGGAACCCCCACCACAGAAAAAGGGGGGGGCAACCCCAGCCGCGGAGCGATAATCCCCTGCTTATAAAACAGCACttg                                                                  >  2:579677/1‑149 (MQ=255)
                                                                                        aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACtt                                                          <  2:578439/151‑1 (MQ=255)
                                                                                         aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                         <  1:1082601/151‑1 (MQ=255)
                                                                                            ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                      <  1:687657/151‑1 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACCGTAACCCCCCCATCAGAAAGCGGGGGGGGCACACCCAGCTGCGGAGCGATCATCCCCTTCTTATAAAACAACACCGGCACCACCTCCCggg                                                    >  1:1165285/1‑151 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                                    <  2:625876/151‑1 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                                    <  2:26859/151‑1 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGGGGGAAAAATAACACACCGGAACCCCCACCTCAGAAAACGGGGGGCGCCCACCGAGCTGCGGAGAGGTCATCCCTTTCTGGTAAAAAACCACCGGCACCACTTCGCggg                                                    >  2:1034310/1‑151 (MQ=255)
                                                                                               gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                                    <  2:446997/150‑1 (MQ=255)
                                                                                               gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                                   <  1:601951/151‑1 (MQ=255)
                                                                                                    gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                              <  1:530125/151‑1 (MQ=255)
                                                                                                     aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAAACACACACGGTAGCCCCCCCATCCGAAAACGGGGGGGGCACACCCAGCCGGGGGAGGAGCATCCCCTTCTTGTAAAACAACACTGGCACCCCTTCCCGGGGAtgggg                                             >  2:247897/1‑149 (MQ=255)
                                                                                                      aagaaaagaGCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATAAAACACCGGAGCCCCCACATCAGAAAGCGGGGGGCGCACCCCGAGCTGCGTAGCGATCATCCCTTTCTTATAAAACAGCACCGGGACCACCTCCCGGGGGtggtgg                                            >  2:686552/1‑151 (MQ=255)
                                                                                                       agaaaagaGCGACCATTTTGTGGCGGAGGCGCGCGGAAAAATTAAAAAACGGGTACCCCCCCAACAGCAAGAGCGGGGGGCACACCGCGCCCCGGAGGGGTCATTCCCTTCTTGTAAAACCACACCGGGACCACCTCGGGGGtgttg                                               >  2:57126/1‑143 (MQ=255)
                                                                                                        gaaaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATAACAAACGGCAGCCCCCCCATCAGAAAACGGGGGGGGCACACCCCGCCGCGGAGCGAGCATCCCTTGCTGGTAAAACAACACCGGCACCCCATCCCGGGGAtggtgg                                            >  2:405569/1‑149 (MQ=255)
                                                                                                            agagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGa                                      <  1:780558/151‑1 (MQ=255)
                                                                                                              agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                                    <  1:1132264/151‑1 (MQ=255)
                                                                                                               gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                                   <  1:1028995/151‑1 (MQ=255)
                                                                                                                cGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGAcc                                    >  1:29383/1‑149 (MQ=255)
                                                                                                                 gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                                 <  2:767892/151‑1 (MQ=255)
                                                                                                                  aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                                <  1:198319/151‑1 (MQ=255)
                                                                                                                     aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                             <  1:1082972/151‑1 (MQ=255)
                                                                                                                       tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                            <  1:1159495/150‑1 (MQ=255)
                                                                                                                       tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                            <  2:921095/150‑1 (MQ=255)
                                                                                                                         tgtgGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGa                          <  2:411425/150‑1 (MQ=255)
                                                                                                                             gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt                     <  1:1019408/151‑1 (MQ=255)
                                                                                                                             gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt                     <  1:1135350/151‑1 (MQ=255)
                                                                                                                              cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGCAGCGTCCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                       >  1:411088/1‑148 (MQ=255)
                                                                                                                                  ggcgCGCGGAAAACTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTCGCTGTGACCAATGGTCGGATCGTTGCCGt                >  2:429092/1‑151 (MQ=255)
                                                                                                                                  ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg                 <  1:29894/150‑1 (MQ=255)
                                                                                                                                  ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt                <  1:1225080/151‑1 (MQ=255)
                                                                                                                                   gcgcgcGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg               <  2:856242/151‑1 (MQ=255)
                                                                                                                                    cgcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg               >  2:354857/1‑150 (MQ=255)
                                                                                                                                            aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc        >  2:208705/1‑150 (MQ=255)
                                                                                                                                             aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCAGCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt      <  2:944075/151‑1 (MQ=255)
                                                                                                                                             aaaaTTCACACACAGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt      <  2:628751/151‑1 (MQ=255)
                                                                                                                                               aaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCa    <  1:79328/151‑1 (MQ=255)
                                                                                                                                                aTTCACACACGGTAACCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGCATGGCGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAacac   >  2:879325/1‑151 (MQ=255)
                                                                                                                                                 ttCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACg  >  1:926832/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                         
TGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACG  >  NC_000913/3508518‑3508812

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