Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,699,320 T→C A369A (GCT→GCC bcsG → DUF3260 family cellulose production inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,699,3200TC98.4% 196.8 / ‑3.9 62A369A (GCT→GCCbcsGDUF3260 family cellulose production inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (34/26);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CCTGATGATGGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCTGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAG  >  NC_000913/3699173‑3699458
                                                                                                                                                   |                                                                                                                                          
ccTGATGATGGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCcgtg                                                                                                                                         <  2:443859/151‑1 (MQ=255)
 cTGATGATGGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCcgtg                                                                                                                                         >  2:759296/1‑150 (MQ=255)
     tgatgGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTa                                                                                                                                     <  1:38892/150‑1 (MQ=255)
     tgatgGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAAAAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGGTCGCCGGTTTATGACGATAACGCCGTGCTTa                                                                                                                                     >  1:451393/1‑150 (MQ=255)
      gatgGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCTTCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCTGTGCTTAAc                                                                                                                                   <  2:480758/151‑1 (MQ=255)
      gatgGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGGTTATGACGATACCGCCGTGCTTaa                                                                                                                                    >  1:342764/1‑150 (MQ=255)
       atgGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAAcc                                                                                                                                  <  2:102950/151‑1 (MQ=255)
        tgGGACATAACCGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGACGGCATGCAGGGCGAATTGATGGATCAAACAAATCTCCCGGTTGTTTTGCTGGGCTTTGATGATTCGCCGGTTTATGACGATACCGCCGTGCTTAAcc                                                                                                                                  <  1:1417/150‑1 (MQ=255)
         gggACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCg                                                                                                                                 <  1:139748/150‑1 (MQ=255)
             cATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCgctgg                                                                                                                             >  1:601421/1‑150 (MQ=255)
             cATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCgctgg                                                                                                                             <  1:583728/150‑1 (MQ=255)
                aaCGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGctggctgg                                                                                                                         <  2:137302/151‑1 (MQ=255)
                              cggTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGaaaaa                                                                                                           <  2:475213/151‑1 (MQ=255)
                               ggTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAg                                                                                                          <  1:714384/151‑1 (MQ=255)
                                     ttGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATaaa                                                                                                     >  2:824531/1‑150 (MQ=255)
                                     ttGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATaaa                                                                                                     >  2:383395/1‑150 (MQ=255)
                                      tGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATaaa                                                                                                     <  1:554895/149‑1 (MQ=255)
                                      tGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATaaa                                                                                                     >  2:554933/1‑149 (MQ=255)
                                        aaaGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACa                                                                                                 <  1:473238/151‑1 (MQ=255)
                                        aaaGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACa                                                                                                 <  1:241110/151‑1 (MQ=255)
                                         aagaagTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAg                                                                                                <  2:345809/151‑1 (MQ=255)
                                          agaagTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGc                                                                                               >  1:228297/1‑151 (MQ=255)
                                            aagTTCGCGAAAATGGCGGCATGCAGGGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGCGCTTTGATGGGTCGCCGGTTTATGCCGATACCGCCGTGCTTAACCGCTGGCTGGACGATACCGAAAAAGATAAAAAc                                                                                                  >  1:685508/1‑146 (MQ=255)
                                               ttCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTa                                                                                           <  2:479397/150‑1 (MQ=255)
                                                    gAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGcc                                                                                      >  2:374242/1‑150 (MQ=255)
                                                    gAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCa                                                                                     <  1:546555/151‑1 (MQ=255)
                                                      aaaTGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACg                                                                                   >  2:298009/1‑151 (MQ=255)
                                                          ggcggcATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCt                                                                               <  2:110646/151‑1 (MQ=255)
                                                           gcggcATGCAGAGCGACTTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATTAAAACAGCCGCAGCGCCACGTTCt                                                                               >  1:643904/1‑150 (MQ=255)
                                                             ggcATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTaca                                                                            <  2:394201/151‑1 (MQ=255)
                                                                  gCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACg                                                                       >  1:255274/1‑151 (MQ=255)
                                                                   cAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTCCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAacac                                                                        <  2:482158/149‑1 (MQ=255)
                                                                     gagCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGCCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGc                                                                      >  2:628993/1‑149 (MQ=255)
                                                                      agCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGTTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTc                                                                   >  2:811479/1‑151 (MQ=255)
                                                                       gCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTcc                                                                  >  1:44480/1‑151 (MQ=255)
                                                                        cGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTcc                                                                  >  1:26579/1‑150 (MQ=255)
                                                                        agtactGATCGATCAAACAAATCTGACGTCTATTTGGTGGGCCTTTGATGGCTCGCCGGTTGATGCCGAAACAGCAGGACTTAACCGCTGGCCGGGCGGTCCCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTc                                                                   <  1:534062/144‑1 (MQ=255)
                                                                         gAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCAc                                                                >  2:548175/1‑151 (MQ=255)
                                                                           aTTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTg                                                              >  1:483695/1‑151 (MQ=255)
                                                                            ttGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTGGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTg                                                              <  1:388288/150‑1 (MQ=255)
                                                                              gATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCa                                                            <  2:557137/150‑1 (MQ=255)
                                                                                 ggATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGa                                                         <  1:120953/150‑1 (MQ=255)
                                                                                     caaacaaaTCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGc                                                     >  1:844407/1‑150 (MQ=255)
                                                                                         caaaTCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAAcc                                                 >  2:640675/1‑150 (MQ=255)
                                                                                              cTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATc                                           >  2:433134/1‑151 (MQ=255)
                                                                                                  cGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCgcg                                       >  2:693754/1‑149 (MQ=255)
                                                                                                    gTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCggggg                                     >  1:173004/1‑151 (MQ=255)
                                                                                                       aTTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTc                                   >  1:760904/1‑150 (MQ=255)
                                                                                                       aTTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCa                                  <  2:647642/151‑1 (MQ=255)
                                                                                                        ttttGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTc                                   >  2:67608/1‑149 (MQ=255)
                                                                                                           ttcttggttttGATGGCTCGTTGTTTTATAATGATAACGTCATTCTTAAACGCTGGCTGGACGTTACCGAATAATATAAAACCAGCAGTAATCCCCCGTTCAACAACATGGTTTCACTTCATGATGGCAACCATGATCCTGGTGTTAGCa                               >  1:749721/9‑150 (MQ=255)
                                                                                                            gCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCaa                              >  2:293696/1‑150 (MQ=255)
                                                                                                                ggCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAAc                           >  1:147125/1‑149 (MQ=255)
                                                                                                                   tttGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCgg                      >  2:607720/1‑151 (MQ=255)
                                                                                                                      gATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGAt                    <  2:601460/150‑1 (MQ=255)
                                                                                                                        tGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTa                  <  2:164020/150‑1 (MQ=255)
                                                                                                                         ggTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTAc                 >  2:72804/1‑150 (MQ=255)
                                                                                                                              gCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAgcg           <  2:769527/151‑1 (MQ=255)
                                                                                                                               ccGGTTTATGACGATACCCCCGTGCTTAACCGCTGGTTGGCCGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAgcgc          >  1:2361/1‑151 (MQ=255)
                                                                                                                                cGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGcgcg         >  2:788644/1‑151 (MQ=255)
                                                                                                                                  gTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCggg       >  1:765394/1‑151 (MQ=255)
                                                                                                                                     tATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGgcg     <  2:173027/150‑1 (MQ=255)
                                                                                                                                        gACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGgcgc    >  2:37105/1‑148 (MQ=255)
                                                                                                                                        gACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAg  <  1:206738/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CCTGATGATGGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCTGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAG  >  NC_000913/3699173‑3699458

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