Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,692,976 T→C 7.6% T430A (ACG→GCG)  uidB ← glucuronide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,692,9760TC 7.6% 260.4 / 10.2 92T430A (ACG→GCG) uidBglucuronide transporter
Reads supporting (aligned to +/- strand):  ref base T (36/49);  new base C (5/2);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

GGCCATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAATTG  >  NC_000913/1692828‑1693125
                                                                                                                                                    |                                                                                                                                                     
ggCCATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTg                                                                                                                                                      <  1:1551897/150‑1 (MQ=255)
   cATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGc                                                                                                                                                   >  2:322996/1‑150 (MQ=255)
    aTGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCg                                                                                                                                                  <  1:754223/150‑1 (MQ=255)
     tGGCGACTATTTTTCTCATAGTTGCAGTTTTGATGTTCTGATTATTTTTTTTGAATATTAATTAGGGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGTTTATCAATTTCAACCACGATTTCTTTTAATTTTTTATCCGTGAGCgg                                                                                                                                                 <  2:1083012/150‑1 (MQ=255)
        ccACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCTCTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATa                                                                                                                                              <  2:1097717/150‑1 (MQ=255)
         cACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATaa                                                                                                                                             <  2:1296461/150‑1 (MQ=255)
           cTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATATAAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATatat                                                                                                                                           >  1:821973/1‑147 (MQ=255)
            tATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAAc                                                                                                                                          <  2:955421/150‑1 (MQ=255)
            tATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAAc                                                                                                                                          <  1:1675761/150‑1 (MQ=255)
            tATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaa                                                                                                                                          <  2:1113132/150‑2 (MQ=255)
             aTTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAAcc                                                                                                                                         >  2:161978/1‑150 (MQ=255)
              tttttCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCa                                                                                                                                        <  1:428128/150‑1 (MQ=255)
                tttCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGa                                                                                                                                      <  1:695436/150‑1 (MQ=255)
                tttCTCATAGTTCCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTCCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGa                                                                                                                                      <  2:67807/150‑1 (MQ=255)
                  tctcATAGTTTCAACTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGata                                                                                                                                    <  1:257862/150‑1 (MQ=255)
                   ctcATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGataa                                                                                                                                   >  2:551178/1‑150 (MQ=255)
                   ctcATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTAAGATTATCAATTTAAACCACGATTTATTTTAATTTTTTATCCGTTAGCGTATAAAACCAGataa                                                                                                                                   >  1:261168/1‑150 (MQ=255)
                         tttGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGa                                                                                                                             <  2:656144/149‑1 (MQ=255)
                         gTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGa                                                                                                                             >  2:644855/1‑150 (MQ=255)
                              aCCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGcc                                                                                                                        <  1:441908/150‑1 (MQ=255)
                                cTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAg                                                                                                                      <  2:1322819/150‑1 (MQ=255)
                                cTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAg                                                                                                                      >  2:175649/1‑150 (MQ=255)
                                  ttGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTa                                                                                                                    <  1:788992/150‑1 (MQ=255)
                                   tGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAg                                                                                                                   >  1:846608/1‑150 (MQ=255)
                                     aTGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCa                                                                                                                 <  2:1456583/150‑1 (MQ=255)
                                      tGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCAC‑TTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATa                                                                                                               <  1:40162/150‑1 (MQ=255)
                                       gTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATa                                                                                                               <  2:1368921/150‑1 (MQ=255)
                                       gTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATa                                                                                                               <  2:1070748/150‑1 (MQ=255)
                                         tCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACAAGATAATAACGAATGCCAGTAGCATaaa                                                                                                             >  1:1152318/1‑150 (MQ=255)
                                              ttattTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGc                                                                                                        <  1:223920/150‑1 (MQ=255)
                                              ttattTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGc                                                                                                        <  1:1705074/150‑1 (MQ=255)
                                                 tttttATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAg                                                                                                     >  2:734592/1‑150 (MQ=255)
                                                     tATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACTAATGCCAGTAGCATAAATCCGCAAGGTAc                                                                                                 >  1:379676/1‑150 (MQ=255)
                                                        tGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTaa                                                                                              >  1:331067/1‑150 (MQ=255)
                                                        tGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTaa                                                                                              >  1:893375/1‑150 (MQ=255)
                                                           atatTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGAGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGTTACTAAGGc                                                                                           >  1:1374729/1‑150 (MQ=255)
                                                            tatTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCa                                                                                          <  1:1058463/150‑1 (MQ=255)
                                                              ttaattaGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAAt                                                                                        <  1:322996/150‑1 (MQ=255)
                                                               taattaGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAAtt                                                                                       <  1:644855/150‑1 (MQ=255)
                                                               taattaGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAAtt                                                                                       <  2:1502660/150‑1 (MQ=255)
                                                                   taGTGATATCGCTGATTAATTGCTGCTTCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGAtg                                                                                   <  2:239773/150‑1 (MQ=255)
                                                                     gTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGAtgtg                                                                                 >  2:1725068/1‑150 (MQ=255)
                                                                        atatCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCgg                                                                              <  2:369580/150‑1 (MQ=255)
                                                                        atatCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCgg                                                                              <  2:1194954/150‑1 (MQ=255)
                                                                         tatCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGa                                                                             <  2:1316924/150‑1 (MQ=255)
                                                                             gCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATATCGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGcc                                                                         >  1:268739/1‑150 (MQ=255)
                                                                               tGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTAATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCa                                                                       <  2:59368/150‑1 (MQ=255)
                                                                                 aTTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCata                                                                     >  1:1039365/1‑150 (MQ=255)
                                                                                 aTTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCata                                                                     >  1:1666839/1‑150 (MQ=255)
                                                                                   tAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCataat                                                                   <  1:630869/150‑1 (MQ=255)
                                                                                    aaTTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCataata                                                                  <  2:1170549/150‑1 (MQ=255)
                                                                                          tgctgcACTTTTTTACGTTTATCAATTTCAACCACGATTTCTTTTAATTTTTTATCCGTGAGCGGATAAAACCATATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTcc                                                            <  2:434770/150‑1 (MQ=255)
                                                                                            ctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCgg                                                          >  2:1670577/1‑150 (MQ=255)
                                                                                             tgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGc                                                         >  1:1254910/1‑150 (MQ=255)
                                                                                              gcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCg                                                        <  1:1074607/150‑1 (MQ=255)
                                                                                              gcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCg                                                        >  2:1719381/1‑150 (MQ=255)
                                                                                                  ttttttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTg                                                    <  2:800972/150‑1 (MQ=255)
                                                                                                   tttttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGc                                                   >  1:23827/1‑150 (MQ=255)
                                                                                                       tACGATTATCAATTTCAACCACGATTTATTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGATAGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCt                                               <  2:1608386/150‑1 (MQ=255)
                                                                                                              aTCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGc                                        <  2:23827/150‑1 (MQ=255)
                                                                                                              aTCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGc                                        <  1:1556151/150‑1 (MQ=255)
                                                                                                              aTCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGc                                        >  2:666205/1‑150 (MQ=255)
                                                                                                              aTCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGCCGTTTGCACCTGATTGGc                                        >  1:164213/1‑150 (MQ=255)
                                                                                                              aTCAATTTCAACCACGATTTCTTTGAATTTTTAATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGc                                        <  2:646735/150‑1 (MQ=255)
                                                                                                               tCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCg                                       <  1:1306226/150‑1 (MQ=255)
                                                                                                               tCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCg                                       <  2:331067/150‑1 (MQ=255)
                                                                                                                aaaTTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGa                                      <  2:1635448/149‑1 (MQ=255)
                                                                                                                 aaTTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGat                                     >  2:1287228/1‑150 (MQ=255)
                                                                                                                      cAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATc                                >  2:1640899/1‑150 (MQ=255)
                                                                                                                      cAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATc                                >  2:91286/1‑150 (MQ=255)
                                                                                                                      cAACCACGATTTCTTTGAATTTTTTATCCGCGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGttt                                                     <  1:259990/129‑1 (MQ=255)
                                                                                                                      cAACCACGATTTCTTTGAATTTTTTATCCGCGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGttt                                                     >  2:259990/1‑129 (MQ=255)
                                                                                                                            cGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTTCGGCTGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCAGCTTa                          >  2:222084/1‑150 (MQ=255)
                                                                                                                            cGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTa                          <  1:680148/150‑1 (MQ=255)
                                                                                                                               tttCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAAcc                       <  1:323805/150‑1 (MQ=255)
                                                                                                                                   tttGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCaa                   <  1:666205/150‑1 (MQ=255)
                                                                                                                                    ttGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTGGAATAATTCCGCAAGGTACTACGGCAATTGATGTGCGGAGGCCCATAATAACTTCCGGCGTTTGCCCCTGATTGGCGATTTATCCGCTTAACCCCaaa                  >  1:598568/1‑150 (MQ=255)
                                                                                                                                    ttGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCaaa                  <  1:91286/150‑1 (MQ=255)
                                                                                                                                     atctggttttATCCGCTCACGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCaaaa                 >  2:1364092/7‑150 (MQ=255)
                                                                                                                                       ctggttttATCCGCTCACGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATa               >  2:1100913/5‑150 (MQ=255)
                                                                                                                                       ctggttttATCCGCTCACGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATa               >  2:44140/5‑150 (MQ=255)
                                                                                                                                        aTTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATaa              >  1:370037/1‑150 (MQ=255)
                                                                                                                                        aTTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATaa              >  1:211892/1‑150 (MQ=255)
                                                                                                                                          gttttATCCGCTCACGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGACCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAAACCCAATATAAAg            >  1:209124/2‑150 (MQ=255)
                                                                                                                                           ttttATCCGTGTGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCATGGTACTAATGCAATTGATGTGCGGATGCCCAAAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAgg           <  2:273262/150‑1 (MQ=255)
                                                                                                                                           ttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAgg           >  2:427827/1‑150 (MQ=255)
                                                                                                                                            tttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGc          <  2:1019496/150‑1 (MQ=255)
                                                                                                                                               aTCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCCCCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAgg       >  2:70520/1‑150 (MQ=255)
                                                                                                                                               aTCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAACATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAgg       <  2:145135/150‑1 (MQ=255)
                                                                                                                                                  cGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAAt    >  1:750746/1‑150 (MQ=255)
                                                                                                                                                   gTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAAtt   >  2:1521695/1‑150 (MQ=255)
                                                                                                                                                   gTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAAtt   >  2:114178/1‑150 (MQ=255)
                                                                                                                                                    tGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAATTg  <  2:846608/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAGGAATTG  >  NC_000913/1692828‑1693125

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑