Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,984,257 C→T 7.5% V271V (GTG→GTA araG ← L‑arabinose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,984,2570CT 7.5% 250.0 / 8.4 80V271V (GTG→GTAaraGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (41/33);  new base T (2/4);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

ATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTG  >  NC_000913/1984111‑1984404
                                                                                                                                                  |                                                                                                                                                   
atcgGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTacgca                                                                                                                                                    >  2:160669/1‑148 (MQ=255)
atcgGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTAGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTacgca                                                                                                                                                    >  1:283505/1‑148 (MQ=255)
     cTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTGCGTACGcca                                                                                                                                               <  1:924381/148‑2 (MQ=255)
     cTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTGCGTACGcca                                                                                                                                               <  1:781098/148‑2 (MQ=255)
     cTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCt                                                                                                                                               <  2:941653/148‑1 (MQ=255)
     cTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTg                                                                                                                                              <  4:107707/149‑1 (MQ=255)
       gTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGt                                                                                                                                            >  3:879048/1‑149 (MQ=255)
        ttGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTGCGTACGccaataa                                                                                                                                           <  4:356358/149‑6 (MQ=255)
           gTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCtt                                                                                                                                        <  3:706837/149‑1 (MQ=255)
             cGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCttt                                                                                                                                       >  1:468600/1‑148 (MQ=255)
               atatAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTcac                                                                                                                                    >  4:941098/1‑149 (MQ=255)
               atatAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTca                                                                                                                                     >  2:563366/1‑148 (MQ=255)
                 ataAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTGCGTACGccaat                                                                                                                                             >  3:642804/1‑135 (MQ=255)
                 ataAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTGCGTACGccaat                                                                                                                                             <  4:642804/138‑4 (MQ=255)
                 ataAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTcaca                                                                                                                                   <  2:267293/148‑1 (MQ=255)
                 ataAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCACCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTcaca                                                                                                                                   >  2:247480/1‑148 (MQ=255)
                   aaaCCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCa                                                                                                                                >  3:925215/1‑149 (MQ=255)
                   aaaCCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCa                                                                                                                                >  3:279181/1‑149 (MQ=255)
                        tGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCaa                                                                                                                            <  2:335552/148‑1 (MQ=255)
                         gACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAg                                                                                                                           <  1:936815/148‑1 (MQ=255)
                         gACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAg                                                                                                                           <  2:177208/148‑1 (MQ=255)
                          aCCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGa                                                                                                                          <  1:901633/148‑1 (MQ=255)
                             ggcggTGGTTTGAGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTa                                                                                                                      >  4:478386/1‑149 (MQ=255)
                              gcggTGATTTGCGTCCCGCCAAACCTGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACACCCCAACAATTTCATCACTGCGAACCGCCAGACTTATTGGCCTACGCACGCCTGGTGCTTTCACAGCACCAAGACGTa                                                                                                                      <  2:881512/148‑1 (MQ=255)
                              gcggTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAg                                                                                                                     <  3:578511/149‑1 (MQ=255)
                               cggTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGTCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAg                                                                                                                     <  2:283505/148‑1 (MQ=255)
                                     tttGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGctc                                                                                                               >  2:625811/1‑148 (MQ=255)
                                          gTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAGCAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGcc                                                                                                         >  4:441217/1‑149 (MQ=255)
                                            cccGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCa                                                                                                        <  1:28357/148‑1 (MQ=255)
                                            cccGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCa                                                                                                        >  2:28357/1‑148 (MQ=255)
                                                   aaCATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTAc                                                                                                 >  2:962410/1‑148 (MQ=255)
                                                     cATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCCCTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTAcgcg                                                                                              <  4:67341/149‑1 (MQ=255)
                                                      aTGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTAcgcg                                                                                              <  1:896662/148‑1 (MQ=255)
                                                       tGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCgg                                                                                             <  1:919702/148‑1 (MQ=255)
                                                        gCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAACCAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGGGGtt                                                                                           <  4:714287/149‑1 (MQ=255)
                                                        gCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGtt                                                                                           <  4:739257/149‑1 (MQ=255)
                                                           tttCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGcc                                                                                        <  4:879048/149‑1 (MQ=255)
                                                            ttCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCgg                                                                                             >  2:363624/1‑143 (MQ=255)
                                                            ttCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCgg                                                                                             <  1:363624/143‑1 (MQ=255)
                                                               aTTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGcc                                                                                    >  3:836198/1‑149 (MQ=255)
                                                                ttAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGcc                                                                                    >  1:28601/1‑148 (MQ=255)
                                                                ttAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGcc                                                                                    >  2:963305/1‑148 (MQ=255)
                                                                ttAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGcc                                                                                    >  2:17842/1‑148 (MQ=255)
                                                                 tAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCg                                                                                   >  2:315694/1‑148 (MQ=255)
                                                                    ttCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAg                                                                                <  2:331317/148‑1 (MQ=255)
                                                                        cTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGatat                                                                            >  1:897126/1‑148 (MQ=255)
                                                                         tACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCg                                                                          >  4:43381/1‑149 (MQ=255)
                                                                           cgccccgcccCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGcc                                                                        <  3:757075/149‑1 (MQ=255)
                                                                             ccccgccccTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCa                                                                       >  1:908832/1‑148 (MQ=255)
                                                                               ccgccccTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATg                                                                    >  3:581161/1‑149 (MQ=255)
                                                                                    ccTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTcgc                                                                >  2:701109/1‑148 (MQ=255)
                                                                                    ccTACAAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGCGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTcgcg                                                               >  4:215417/1‑149 (MQ=255)
                                                                                      tACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGcc                                                             >  4:501852/1‑149 (MQ=255)
                                                                                       aCCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGcc                                                             >  2:14044/1‑148 (MQ=255)
                                                                                       aCCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGcc                                                             >  2:80382/1‑148 (MQ=255)
                                                                                          aGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGa                                                          >  2:139382/1‑148 (MQ=255)
                                                                                                 aCAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGc                                                   >  2:860903/1‑148 (MQ=255)
                                                                                                  cAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCt                                                 <  3:424073/149‑1 (MQ=255)
                                                                                                   aGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGTACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTg                                                >  3:506791/1‑149 (MQ=255)
                                                                                                      ccAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCAc                                             <  4:668485/149‑1 (MQ=255)
                                                                                                        aacaaTTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCa                                           >  3:163061/1‑149 (MQ=255)
                                                                                                         acaaTTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAg                                          >  4:363344/1‑149 (MQ=255)
                                                                                                            aTTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAAt                                                                     <  4:917873/119‑1 (MQ=255)
                                                                                                            aTTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAAt                                                                     >  3:917873/1‑119 (MQ=255)
                                                                                                             tttCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAgcgc                                       >  1:787961/1‑148 (MQ=255)
                                                                                                                 accacTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGcgtcgt                                  >  3:390751/1‑149 (MQ=255)
                                                                                                                  ccacTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGcgtcgt                                  >  1:294301/1‑148 (MQ=255)
                                                                                                                  ccacTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTg                                 <  3:941098/149‑1 (MQ=255)
                                                                                                                   cacTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTg                                 <  1:639245/148‑1 (MQ=255)
                                                                                                                        cGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCaa                            <  1:563366/148‑1 (MQ=255)
                                                                                                                          aaCCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGAcc                                                        <  4:893052/118‑1 (MQ=255)
                                                                                                                          aaCCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGAcc                                                        >  3:893052/1‑118 (MQ=255)
                                                                                                                          aaCCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACct                         >  4:276976/1‑149 (MQ=255)
                                                                                                                               ccAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCtgctgc                    >  3:171715/1‑149 (MQ=255)
                                                                                                                                 aGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGcc                                                  <  4:904373/117‑1 (MQ=255)
                                                                                                                                 aGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGcc                                                  >  3:904373/1‑117 (MQ=255)
                                                                                                                                        ttGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTa           <  3:441217/149‑1 (MQ=255)
                                                                                                                                         tGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTaa          <  3:363344/149‑1 (MQ=255)
                                                                                                                                          ggCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTaa          <  1:854522/148‑1 (MQ=255)
                                                                                                                                           gCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCGCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTaaa         <  1:821359/148‑1 (MQ=255)
                                                                                                                                            cGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAg        >  1:671396/1‑148 (MQ=255)
                                                                                                                                             gTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGt      >  3:509773/1‑149 (MQ=255)
                                                                                                                                               acgcacgcCTGGTGCTTTCACAGCATCAAGCCGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGttt    <  4:836198/149‑1 (MQ=255)
                                                                                                                                               acgcacgcCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGttt    <  3:43381/149‑1 (MQ=255)
                                                                                                                                                cgcacgcCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCACCATACCTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGCCCATCGTCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGtttt   <  4:59591/149‑1 (MQ=255)
                                                                                                                                                 gcacgcCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTg  <  3:276976/149‑1 (MQ=255)
                                                                                                                                                 gcacgcCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTg  >  3:913930/1‑149 (MQ=255)
                                                                                                                                                  cacgcCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGtt     >  1:294865/1‑145 (MQ=255)
                                                                                                                                                  cacgcCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGtt     <  2:294865/145‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTG  >  NC_000913/1984111‑1984404

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