Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,734 C→G 5.8% A415A (GCC→GCG lsrA → autoinducer 2 import ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,602,7340CG 5.8% 235.1 / 5.5 69A415A (GCC→GCGlsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base G (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

TCATAACCTTCGTGGATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCTGTGCTGCTTATCTCCTCCGACCT  >  NC_000913/1602587‑1602877
                                                                                                                                                   |                                                                                                                                               
tCATAACCTTCGTGGATTCTGGGTGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTACACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGcc                                                                                                                                                 >  3:161949/1‑148 (MQ=255)
  aTAACCTTCGTGGATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCaaa                                                                                                                                              >  7:12399/1‑149 (MQ=255)
  aTAACCTTCGTGGATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCaaa                                                                                                                                              >  5:114533/1‑149 (MQ=255)
     aCCTTCGTGGATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATg                                                                                                                                            <  6:8131/148‑1 (MQ=255)
       cTTCGTGGATTCTGGGTGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCtt                                                                                                                                         <  8:33689/149‑1 (MQ=255)
           gTGGATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGaa                                                                                                                                     <  4:129400/149‑1 (MQ=255)
           gTGGATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGaa                                                                                                                                     <  5:86516/149‑1 (MQ=255)
            tGGATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCGAAATGCTTGGAAg                                                                                                                                    >  3:307985/1‑149 (MQ=255)
             ggATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGc                                                                                                                                   >  7:208069/1‑149 (MQ=255)
             ggATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACAAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGc                                                                                                                                   >  2:277487/1‑149 (MQ=255)
                       cGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCGAAATGCTTGGAAGCTTCGCCGCaa                                                                                                                         <  4:307985/149‑1 (MQ=255)
                       cGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCaa                                                                                                                         <  8:208069/149‑1 (MQ=255)
                          aaaCCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATAATAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGAAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTa                                                                                                                      >  7:377652/1‑149 (MQ=255)
                            aCCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTatt                                                                                                                    <  6:400800/149‑1 (MQ=255)
                               gcgAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACCAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTcttgat                                                                                                                 <  2:126056/149‑1 (MQ=255)
                               gcgAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCGAAATGCTTGGAAGCTTCGCCGCAAGTattgat                                                                                                                 >  7:166272/1‑149 (MQ=255)
                                cgAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTattgatt                                                                                                                >  2:38875/1‑149 (MQ=255)
                                 gAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTattgattg                                                                                                               <  2:405822/149‑1 (MQ=255)
                                   aaGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTc                                                                                                             <  8:307469/149‑1 (MQ=255)
                                           gCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCaa                                                                                                                         <  3:363594/129‑1 (MQ=255)
                                           gCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCaa                                                                                                                         >  4:363594/1‑129 (MQ=255)
                                             cACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGa                                                                                                   >  6:401519/1‑149 (MQ=255)
                                                ccTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGAcgc                                                                                                <  8:193303/149‑1 (MQ=255)
                                                  tGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTgg                                                                                                                                       >  8:341590/1‑108 (MQ=255)
                                                  tGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTgg                                                                                                                                       <  7:341590/108‑1 (MQ=255)
                                                        gTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTg                                                                                                                                        >  8:129348/1‑101 (MQ=255)
                                                        gTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTg                                                                                                                                        <  7:129348/101‑1 (MQ=255)
                                                        gTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGACGCAAGTATTGATTGTCGATGAGCCGACGCGCGCCGTg                                                                                        >  1:404945/1‑149 (MQ=255)
                                                                 gggCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCg                                                                               <  5:23571/149‑1 (MQ=255)
                                                                     gcTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGccc                                                                           <  4:353778/149‑1 (MQ=255)
                                                                         aaTTTTAAATTCAACCAACCGGATCAAGCTGCACTGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGAcgcgcg                                                                                             <  4:98550/127‑1 (MQ=255)
                                                                         aaTATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTaa                                                                       >  7:300615/1‑149 (MQ=255)
                                                                         aaTATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGCTGTCAACCAGCAAAAAATACTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTAATGATTGTCCATGAGCCGAcgcgcg                                                                                             >  3:98550/1‑127 (MQ=255)
                                                                           tatTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAAt                                                                      >  8:99538/1‑148 (MQ=255)
                                                                            atTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCaa                                                                                                                         >  3:134514/1‑96 (MQ=255)
                                                                            atTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCaa                                                                                                                         <  4:134514/96‑1 (MQ=255)
                                                                             ttAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGa                                                                    <  6:114533/148‑1 (MQ=255)
                                                                               aaaTTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGatat                                                                 >  7:83345/1‑149 (MQ=255)
                                                                                aaTTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGatat                                                                 <  2:254768/148‑1 (MQ=255)
                                                                                       ccaaccGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGc                                                         <  1:277487/149‑1 (MQ=255)
                                                                                        caaccGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCGAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCt                                                        <  8:166272/149‑1 (MQ=255)
                                                                                         aaccGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTg                                                       >  4:172210/1‑149 (MQ=255)
                                                                                                    gCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATg                                                                     >  5:114670/1‑124 (MQ=255)
                                                                                                    gCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATg                                                                     <  6:114670/124‑1 (MQ=255)
                                                                                                    gCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCAt                                            >  7:196396/1‑149 (MQ=255)
                                                                                                     cTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGAGGGCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATc                                           <  3:6858/149‑1 (MQ=255)
                                                                                                         aCGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATcgccg                                       >  1:52897/1‑149 (MQ=255)
                                                                                                         aCGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATcgccg                                       >  7:247264/1‑149 (MQ=255)
                                                                                                           ggACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCTCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGca                                     >  1:345826/1‑149 (MQ=255)
                                                                                                           ggACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGca                                     >  7:323446/1‑149 (MQ=255)
                                                                                                           ggACATTATCCGGTGGCAACCAGCAAAAAATCCTCATAGCCAAAAGCTTGGAAGCTTCGACGCAAGAATTGATTGTCGATGAGCCGACGCTCGGCGTGGATGTCTCGGCCCGTAATCATATATACCAGATGTTGCCCAGCATACCCGca                                     >  3:132394/1‑149 (MQ=255)
                                                                                                                ttATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATCAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACaaaa                                >  3:110310/1‑149 (MQ=255)
                                                                                                                   tCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAAtgt                             <  8:247264/149‑1 (MQ=255)
                                                                                                                   tCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAAtgt                             <  3:172210/149‑1 (MQ=255)
                                                                                                                     cGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGCGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTgg                           >  3:361114/1‑149 (MQ=255)
                                                                                                                     cGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTgg                           <  1:38875/149‑1 (MQ=255)
                                                                                                                       gTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCt                         <  8:30079/149‑1 (MQ=255)
                                                                                                                       gTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCt                         >  7:30079/1‑149 (MQ=255)
                                                                                                                          gCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCtgtg                      <  8:323446/149‑1 (MQ=255)
                                                                                                                          gCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCtgtg                      <  5:401519/149‑1 (MQ=255)
                                                                                                                            aaCCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCTGtgct                    >  2:252384/1‑149 (MQ=255)
                                                                                                                                   aaaaaaTCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCTGTGCTGCTTAtc             >  8:2881/1‑149 (MQ=255)
                                                                                                                                   aaaaaaTCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCTGTGCTGCTTAt              >  4:251764/1‑148 (MQ=255)
                                                                                                                                    aaaaaTCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCgca                                                >  5:253235/1‑113 (MQ=255)
                                                                                                                                    aaaaaTCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCgca                                                <  6:253235/113‑1 (MQ=255)
                                                                                                                                            tCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGGATCGCCGCACAAAATGTGGCTGTGCTGCTTATCTCCTCCGAc    >  1:95675/1‑149 (MQ=255)
                                                                                                                                              aTTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCTGTGCTGCTTATCTCCTCCGACCt  <  4:161949/149‑1 (MQ=255)
                                                                                                                                              aTTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCTGTGCTGCTTATCTCCTCCGACCt  >  4:98091/1‑149 (MQ=255)
                                                                                                                                                tGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGtt                                                     <  1:72458/96‑1 (MQ=255)
                                                                                                                                                tGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGtt                                                     >  2:72458/1‑96 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCATAACCTTCGTGGATTCTGGGCGAAAACCGCGAAAGATAATGCCACCCTGGAACGTTATCGTCGGGCGCTGAATATTAAATTCAACCAACCGGAACAAGCTGCACGGACATTATCCGGTGGCAACCAGCAAAAAATCCTCATTGCCAAATGCTTGGAAGCTTCGCCGCAAGTATTGATTGTCGATGAGCCGACGCGCGGCGTGGATGTCTCGGCCCGTAATGATATCTACCAGCTGTTGCGCAGCATCGCCGCACAAAATGTGGCTGTGCTGCTTATCTCCTCCGACCT  >  NC_000913/1602587‑1602877

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