Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 366,726 A→G 5.0% L262P (CTG→CCG)  lacI ← lactose‑inducible lac operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,7260AG 5.0% 245.7 / 2.5 80L262P (CTG→CCG) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (43/32);  new base G (2/2);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

GGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGTTCCGCTAT  >  NC_000913/366580‑366874
                                                                                                                                                  |                                                                                                                                                    
ggCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGc                                                                                                                                                    >  4:66074/1‑149 (MQ=255)
 gCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGTGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGc                                                                                                                                                    <  5:74017/148‑1 (MQ=255)
   ccTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTATGCCCCAGCAGGCGAAAATCCAGTTTGATGCTCGTTAACGGCGGGATTTAACATTAGCTGTCTTCGGTATCGTCGTATCCCACTACCGCGATATCCGCACCAACGCGCAGCCCg                                                                                                                                                 >  3:282485/1‑149 (MQ=255)
    cTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCgg                                                                                                                                                >  2:336218/1‑149 (MQ=255)
      gagagaGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGAc                                                                                                                                              >  8:112270/1‑149 (MQ=255)
         agagTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCg                                                                                                                                           <  7:112270/149‑1 (MQ=255)
             ttGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTaa                                                                                                                                       >  5:71082/1‑149 (MQ=255)
              tGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAAt                                                                                                                                      >  4:31091/1‑149 (MQ=255)
                 agcaAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGgc                                                                                                                                   <  5:206081/149‑1 (MQ=255)
                  gcaAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGATGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCATCCCGGACTCGGTAATGgcg                                                                                                                                  >  8:164464/1‑149 (MQ=255)
                   catgCGGTCCACGCTGGTTTGACCCAGCAGGCGAAAATCCTGTTTGATGGAGGTTAACGGCGGGATATAACATGAGCGGGCTTCGGTCTCGGCGTCTCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGgcgc                                                                                                                                 <  5:96289/146‑1 (MQ=255)
                   caAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGgcgc                                                                                                                                 <  1:336218/149‑1 (MQ=255)
                               gCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGcc                                                                                                                       >  5:110193/1‑147 (MQ=255)
                                 tGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAgc                                                                                                                   >  8:58608/1‑149 (MQ=255)
                                       gCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCCACGCGCCGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATc                                                                                                             >  4:26339/1‑149 (MQ=255)
                                        ccccAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATc                                                                                                             <  6:71082/148‑1 (MQ=255)
                                         cccAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCt                                                                                                            >  1:309582/1‑148 (MQ=255)
                                           catcagGCGAAAATCCTGTTTGATGGTGGTTAACGGCGTGATATAACATGAGCTGTCTTCTGTATCGTCGTATCCCACTACCGAGATATCCGCCCCAACGCGCAGCCCGCACTCGGTAATGGCGAGCATTGCGCCCAGCGACATCTGAt                                                                                                         >  8:154763/4‑149 (MQ=255)
                                           cagcagGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGAt                                                                                                         <  2:30925/149‑1 (MQ=255)
                                            agcagGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATc                                                                                                        >  4:45782/1‑149 (MQ=255)
                                             gcagGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCg                                                                                                       >  2:281840/1‑149 (MQ=255)
                                               agGCGAAAATCCTGTTTGATGGTGGTTAACGGCGTGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGtt                                                                                                     <  7:164464/149‑1 (MQ=255)
                                                ggCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTg                                                                                                    <  6:175193/149‑1 (MQ=255)
                                                 gCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACt                                                                                                                                             >  7:88389/1‑107 (MQ=255)
                                                 gCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACt                                                                                                                                             <  8:88389/107‑1 (MQ=255)
                                                   gAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGc                                                                                                  >  6:327896/1‑148 (MQ=255)
                                                      aaTCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAAc                                                                                               >  1:270784/1‑148 (MQ=255)
                                                       aTCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCa                                                                                             <  2:309582/149‑1 (MQ=255)
                                                       aTCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCa                                                                                             <  5:327896/149‑1 (MQ=255)
                                                           tGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCAt                                                                                         >  4:318804/1‑149 (MQ=255)
                                                            gTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATc                                                                                        >  7:158347/1‑149 (MQ=255)
                                                                 aTGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAg                                                                                    <  6:110193/148‑1 (MQ=255)
                                                                   ggtggtTAACGGCGGGATATAACATGAGCTGTCTTCTGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGAGCAGCCAGGACTCTGTAATTGCGCTCATTGCTCCGAGAGCCATCTGATCGTTGGCAACCAGCATCGCAGTg                                                                                  <  5:62300/148‑1 (MQ=255)
                                                                   ggtggtTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTg                                                                                  <  2:175398/148‑1 (MQ=255)
                                                                   gggggtTAACGGCGGCAGATACCTTGCGCTGGCTGCGGCATCGTCCTAGCCCACTACCGAGATCTCCGCCCCAACGCGCAGCCCGGACGCGCTAATGGCGCGCATTGCGCCCAGCGCCATCGGATCGTAGGCAACCAGCATCGCAGTgg                                                                                 <  4:282485/146‑1 (MQ=255)
                                                                       gtTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAAc                                                                             >  2:248963/1‑149 (MQ=255)
                                                                           aCGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATg                                                                         <  1:281840/149‑1 (MQ=255)
                                                                            cggcggGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGc                                                                        <  8:118479/149‑1 (MQ=255)
                                                                              gcggGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGccc                                                                      <  4:26647/149‑1 (MQ=255)
                                                                                gggATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTc                                                                    <  7:74803/149‑1 (MQ=255)
                                                                                  gATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCAt                                                                  >  7:372314/1‑149 (MQ=255)
                                                                                    tatgACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTc                                                                <  4:178505/149‑1 (MQ=255)
                                                                                       aaCATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAg                                                              >  4:301130/1‑148 (MQ=255)
                                                                                        ctaTGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTAGGTAATGGCTCGCATTGGGCCCAGCGCCATCTGATCTTTGGCAACCAGCATCGAAGTTTTAACGATTCCCTCAtt                                                                 <  6:266548/142‑1 (MQ=255)
                                                                                        aCATGAGATGTCTTCGGTATCGTCGTATCCCACAACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGCCGAGCATTGAGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCTCAGTGGGAACGATGCCCTCAtt                                                                 >  5:266548/1‑144 (MQ=255)
                                                                                         cATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCAt                                                           >  7:91256/1‑149 (MQ=255)
                                                                                           tGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCAttt                                                         <  1:248963/149‑1 (MQ=255)
                                                                                           tGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCAttt                                                         <  3:318804/149‑1 (MQ=255)
                                                                                             accTGTCTTCGGTATCGTCGTATCCCACTACCGCGATATCCGCACCACCGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCAGCGCAGTGGGAACGATGCCCTCATTCAGCATTTg                                                        <  2:29736/146‑1 (MQ=255)
                                                                                                  tCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGt                                                  >  6:209152/1‑149 (MQ=255)
                                                                                                    ttCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCCACCCGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGttt                                                >  4:364625/1‑149 (MQ=255)
                                                                                                      cGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTttg                                               <  2:270784/148‑1 (MQ=255)
                                                                                                       ggTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTttgtt                                             <  3:66074/149‑1 (MQ=255)
                                                                                                             gtcgtATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCCTTGCGCCCAGCGCCATCTGCTCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAAc                                       >  2:14981/1‑149 (MQ=255)
                                                                                                              tcgtATCCCACTACCGAGATATCCGCACCAACGCGCGGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATg                                                    <  7:269008/135‑1 (MQ=255)
                                                                                                              tcgtATCCCACTACCGAGATATCCGCACCAACGCGCGGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATg                                                    <  7:262894/135‑1 (MQ=255)
                                                                                                              tcgtATCCCACTACCGAGATATCCGCACCAACGCGCGGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATg                                                    >  8:262894/1‑135 (MQ=255)
                                                                                                              tcgtATCCCACTACCGAGATATCCGCACCAACGCGCGGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATg                                                    >  8:269008/1‑135 (MQ=255)
                                                                                                              tcgtATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAAc                                       >  4:339418/1‑148 (MQ=255)
                                                                                                                gtATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCg                                     <  3:119730/148‑1 (MQ=255)
                                                                                                                 tATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGa                                   <  6:212833/149‑1 (MQ=255)
                                                                                                                            cGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTc                         >  6:33460/1‑148 (MQ=255)
                                                                                                                              agaTATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAg                      >  5:22041/1‑149 (MQ=255)
                                                                                                                              agaTATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAg                      >  7:366121/1‑149 (MQ=255)
                                                                                                                                atatCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTc                    <  5:33460/149‑1 (MQ=255)
                                                                                                                                  atCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGc                  >  5:295125/1‑149 (MQ=255)
                                                                                                                                  atCCGCACCAACGCGCAGCCCGGACTCGGAAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGACCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGACGc                  >  4:289929/1‑149 (MQ=255)
                                                                                                                                   tCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGcc                 <  3:31091/149‑1 (MQ=255)
                                                                                                                                   tCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGcc                 >  5:282410/1‑149 (MQ=255)
                                                                                                                                   tCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGcc                 >  4:181467/1‑149 (MQ=255)
                                                                                                                                   tCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCAGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACAACAGTCGcc                 >  6:152265/1‑149 (MQ=255)
                                                                                                                                    ccGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCACCCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCt                >  7:171135/1‑149 (MQ=255)
                                                                                                                                    ccGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCt                >  4:343400/1‑149 (MQ=255)
                                                                                                                                     cGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCGCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCt                >  5:166708/1‑148 (MQ=255)
                                                                                                                                       cACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTcc             >  4:58880/1‑149 (MQ=255)
                                                                                                                                          cAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGt          >  1:271933/1‑149 (MQ=255)
                                                                                                                                             cgcgCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGTTCCg      >  7:124570/1‑149 (MQ=255)
                                                                                                                                                gcAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCATCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGTTCCGCt    >  6:131586/1‑149 (MQ=255)
                                                                                                                                                gcAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGtt                                                 <  3:97797/104‑1 (MQ=255)
                                                                                                                                                gcAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGtt                                                 >  4:97797/1‑104 (MQ=255)
                                                                                                                                                  aGCCCGGACTCGGTAATGGCGCGCATTGCTCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGTTCCGCTAt  <  5:143555/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGTTCCGCTAT  >  NC_000913/366580‑366874

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: