Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,826,1880TC51.2% 0.2 / 107.2 82T246A (ACA→GCA) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (13/26);  new base C (17/24);  total (32/50)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

TTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTG  >  NC_000913/1826041‑1826333
                                                                                                                                                   |                                                                                                                                                 
ttACGGTAAAACGGTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGgtgt                                                                                                                                                   >  1:630146/1‑148 (MQ=255)
 tACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCtc                                                                                                                                                 >  2:1258292/1‑149 (MQ=255)
  aCGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCa                                                                                                                                              <  2:73097/151‑1 (MQ=255)
   cGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCAc                                                                                                                                             <  2:979354/151‑1 (MQ=255)
   cGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCAc                                                                                                                                             <  1:346501/151‑1 (MQ=255)
    ggTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGGGTtc                                                                                                                                                 >  2:992056/1‑146 (MQ=255)
       aaaaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATa                                                                                                                                         <  1:856881/151‑1 (MQ=255)
       aaaaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGGGTTCTTACGAt                                                                                                                                          >  1:1255501/1‑150 (MQ=255)
        aaaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTGGCGAAACCACCCGATAACGGAGCGGCTGTGCTCTCAc                                                                                                                                             <  2:1058096/146‑1 (MQ=255)
         aaCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATGCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATAcc                                                                                                                                       <  1:54190/151‑1 (MQ=255)
          aCGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCCCGATAcc                                                                                                                                       >  2:813125/1‑150 (MQ=255)
            gTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTTGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGac                                                                                                                                    <  2:153757/151‑1 (MQ=255)
               cccctcCGCCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATAGGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGTGTGGCGGTGCTATCACGATACCGACAct                                                                                                                                 <  1:483719/148‑1 (MQ=255)
               ccTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGCGCCCTCACGATAACGacac                                                                                                                                  >  2:1241136/1‑150 (MQ=255)
                 tctcCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGGAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGGGCTCTCACGGTACCGACActc                                                                                                                                >  1:632685/1‑150 (MQ=255)
                       tCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTTTCCAGATAACGACACCCTCaac                                                                                                                           >  1:1090853/1‑147 (MQ=255)
                        ccTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGa                                                                                                                        <  1:458348/151‑1 (MQ=255)
                        ccTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGa                                                                                                                        <  1:1113165/151‑1 (MQ=255)
                         cggcgcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCCAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGAt                                                                                                                       <  1:630434/149‑1 (MQ=255)
                           gcgcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGGTCTCACGATACCGACACCCCCACCAGa                                                                                                                        >  2:424182/1‑148 (MQ=255)
                                 gCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAgcagcgc                                                                                                               <  2:842406/151‑1 (MQ=255)
                                  cAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAgcagcgca                                                                                                              <  2:1105292/151‑1 (MQ=255)
                                  cAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAgcagcgca                                                                                                              <  2:713765/151‑1 (MQ=255)
                                   aaTGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCCGATAgcagcgc                                                                                                               >  2:573839/1‑149 (MQ=255)
                                    aTGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCGCTCACGGTACCCGCCCTCCCCCCCGACAgcggcgccg                                                                                                             >  1:390670/1‑148 (MQ=255)
                                      gTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCgcagca                                                                                                           <  2:63739/150‑1 (MQ=255)
                                        tCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAAtt                                                                                                        >  2:111022/1‑151 (MQ=255)
                                        tCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAAt                                                                                                         <  2:1258489/150‑1 (MQ=255)
                                         ccTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAAtt                                                                                                        <  2:214351/150‑1 (MQ=255)
                                                cAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGc                                                                                                <  2:560419/151‑1 (MQ=255)
                                                 aaaCGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGATTAATTTGCGAAACCACCCGATAACGGGGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGc                                                                                                <  2:991612/150‑1 (MQ=255)
                                                 aaaCGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCa                                                                                               <  1:349801/151‑1 (MQ=255)
                                                     ggCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTa                                                                                           <  1:972256/151‑1 (MQ=255)
                                                      gCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTAc                                                                                          <  2:459240/151‑1 (MQ=255)
                                                       cATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCCGATAGCAGCGCAGCAATTTCGCTGGCAGtt                                                                                            >  2:296961/1‑148 (MQ=255)
                                                          aaaaTTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCa                                                                                      <  2:124397/151‑1 (MQ=255)
                                                            aaTTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCaaa                                                                                    <  1:975507/151‑1 (MQ=255)
                                                            aaTTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATAACGACACTCTCACCAGATAACAGCGCAGCAATTTCGCTGGCAGTTACTGcca                                                                                     >  2:259245/1‑148 (MQ=255)
                                                            aaTTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCAAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACCCTCTCACCCGAAAACAACGCTGCAATTTCGCTTGCCGTTACAGcca                                                                                     >  2:1182971/1‑148 (MQ=255)
                                                               tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGATACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCaaa                                                                                    >  2:907031/1‑148 (MQ=255)
                                                               tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACt                                                                                  <  2:690958/150‑1 (MQ=255)
                                                                cAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTg                                                                                 <  1:1228654/150‑1 (MQ=255)
                                                                 aGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTggt                                                                               <  2:148142/151‑2 (MQ=255)
                                                                   cGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGa                                                                             <  2:1039422/151‑1 (MQ=255)
                                                                   cGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAAATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCg                                                                              >  2:1250911/1‑150 (MQ=255)
                                                                    gtgtCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGaa                                                                            <  1:61959/151‑1 (MQ=255)
                                                                    gtgtCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGGGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAAATGGCGa                                                                             >  2:721407/1‑150 (MQ=255)
                                                                        cACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGAt                                                                         <  1:1078499/150‑1 (MQ=255)
                                                                         aCTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCg                                                                       <  1:296945/151‑1 (MQ=255)
                                                                               cATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGc                                                                 <  2:984152/151‑1 (MQ=255)
                                                                               cATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGc                                                                 <  2:317609/151‑1 (MQ=255)
                                                                                 tatGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTTccc                                                               >  2:1110202/1‑151 (MQ=255)
                                                                                 tatGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGccc                                                               >  2:878070/1‑151 (MQ=255)
                                                                                   tGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGAAACCGACACCCTCACCCGATAAAAGCGCAGCAATTGCGCTGGCCGTTACTTCAAACTGGCGAAGATCGTTTTTccc                                                               >  2:1058387/1‑149 (MQ=255)
                                                                                     cATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAg                                                           <  2:48090/151‑1 (MQ=255)
                                                                                     cATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACAGAGCGGCGGGGACCTCACGGTACGCACACTCTCCCCCAAAAGTAGCGCCACAAATGCGGTGGAAGCTACTGAAAAATGGCGAAAAACGgttt                                                                   >  1:953905/1‑143 (MQ=255)
                                                                                       tttCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCCACACTCTCACCAGATAGCAGCGCAGCAAATTCGCTGGCAGGTACTGCAAACTGGCGAAGATCGTTTTGCCCCAAggg                                                         >  2:990725/1‑151 (MQ=255)
                                                                                         tCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGGGCCCCCCCGATACCGACACTCTCCCCAGAAAGCAGCGCAGCCATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAagg                                                         >  1:244879/1‑149 (MQ=255)
                                                                                         tCGAAGGACGGCGAGTAGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCTGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCATACTGGAGAAGATCGTTTTGCCCAAAggg                                                         <  1:923621/149‑1 (MQ=255)
                                                                                           gAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAAc                                                     <  2:1187038/151‑1 (MQ=255)
                                                                                             aGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGc                                                   <  1:1083309/151‑1 (MQ=255)
                                                                                              ggACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCt                                                  <  1:1096821/151‑1 (MQ=255)
                                                                                               gACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCtt                                                 <  1:621077/151‑1 (MQ=255)
                                                                                                aCGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCCACACTCTCCCCAGATAGCAGCGCAGCAATTTCGCCGGCAGTTACTGCAAACTGGCGAAGATCGTTTTTCCCCAAGGGCAAAGCttt                                                >  2:357324/1‑151 (MQ=255)
                                                                                                   gcgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGcc                                             <  2:1126647/151‑1 (MQ=255)
                                                                                                    cgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCa                                            <  1:640870/151‑1 (MQ=255)
                                                                                                     gagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCaa                                           <  2:8926/151‑1 (MQ=255)
                                                                                                         ggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCCCCAGATAGCAGCGCAGCAATTGCGCCGGCAGTTACTTCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGtt                                        >  1:429088/1‑150 (MQ=255)
                                                                                                            gagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAAATGGCGAAGAACGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTcc                                      >  1:577712/1‑149 (MQ=255)
                                                                                                             agtAATTTCCGAAACCACCCGATAACGGAGCCGCGGTGTACTCACGATGCCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTcc                                      <  2:630178/148‑1 (MQ=255)
                                                                                                              gtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCCCCCCACATAGCAGAGCAGCAAATGCGCTGGCAGGTTCCGCAAACTGGCGCAGATCGTTTTGGCCCAAGGGAAACGCCTTGCCAAGTTcc                                      >  1:466418/1‑147 (MQ=255)
                                                                                                              gtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCg                                  <  2:813586/151‑1 (MQ=255)
                                                                                                                   ttGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACa                              <  1:877594/150‑1 (MQ=255)
                                                                                                                         aaCCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAgc                       <  2:655249/151‑1 (MQ=255)
                                                                                                                           ccaccCGATAACGGAGCGGCGGTGCTCTCACGATACCGACACCCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCCAACTGGCGAAGATCGTTTTGCCCAAAAGGCAACGCTTTGCCAAGTTCCCGCGTACAGGCCAgcgc                     >  2:681230/1‑151 (MQ=255)
                                                                                                                            caccCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAgcgcg                    <  1:1155904/151‑1 (MQ=255)
                                                                                                                            caccCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAAATGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAgggc                          >  2:1157258/1‑143 (MQ=255)
                                                                                                                               ccGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCa                 <  1:992019/151‑1 (MQ=255)
                                                                                                                               ccGATAACGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGGTCCAGCGTACAGGCCAGCGCGcc                  >  1:537229/1‑150 (MQ=255)
                                                                                                                                    aaCGGAGCGGCGGCGTTCTCACGAAACCGGCACTCTCACCCAATAGCAGCGCCGCAATTGCGCTGGCAGGTACTGCAAACTGGGGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCCGCGCGCCaaaa              >  1:490203/1‑149 (MQ=255)
                                                                                                                                      cGGAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGta          <  1:107008/151‑2 (MQ=255)
                                                                                                                                        gAGCGGCGGTGCTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTg        <  1:178000/151‑1 (MQ=255)
                                                                                                                                               ggTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTg  >  2:473029/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTG  >  NC_000913/1826041‑1826333

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