Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 860,412 C→T G473E (GGG→GAG)  ybiW ← putative formate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,4120CT99.1% 367.8 / ‑5.3 107G473E (GGG→GAG) ybiWputative formate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (42/64);  total (43/64)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 126 total aligned reads displayed.

GCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGTT  >  NC_000913/860322‑860501
                                                                                          |                                                                                         
gcgTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAAGAATTTGTCGCTGGCGGGATCATGCTCg                                                                                         >  2:516692/1‑93 (MQ=255)
 cgTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGGACACTTTTCCGCTGGGTGCATCATGCTCGc                                                                                        >  2:373842/1‑93 (MQ=255)
 cgTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACCTTGCCGCTGGTGGCATCATGCTCGc                                                                                        >  2:1794779/1‑93 (MQ=255)
  gTCCATCACTTCATCGAAGTTGTTGAAGTTCCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGcc                                                                                       <  2:3073423/93‑1 (MQ=255)
  gTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGcc                                                                                       <  2:1032669/93‑1 (MQ=255)
   tCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCt                                                                                      <  1:2388804/93‑1 (MQ=255)
     cATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACCCTTTGCCGCTGGTGGCATCATGCTCGCCTTc                                                                                    >  2:2145388/1‑93 (MQ=255)
     cATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTc                                                                                    >  1:3369109/1‑93 (MQ=255)
       tCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTccc                                                                                  >  1:3106787/1‑92 (MQ=255)
         aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACCTTTACGCTGGTGGCATCATGCGCGCCTTCCCgc                                                                                >  2:2835717/1‑93 (MQ=255)
           ttCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGGTCGCCTTCCAgcgc                                                                              >  2:3814559/1‑93 (MQ=255)
           ttCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAgcgc                                                                              <  2:625538/93‑1 (MQ=255)
              aTCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAAAACTTTGTCGCTGGGGGCATCATGCTCGCCTTCCCGcgccgc                                                                           >  2:3102005/1‑93 (MQ=255)
               tCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGGATCATGCTCGCCTTTCAGCgccgcc                                                                          >  2:3816090/1‑93 (MQ=255)
               tCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCCGCgccgcc                                                                          >  2:267970/1‑93 (MQ=255)
                 gAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGGACACTTTGACGCTGGTGGGATCCTGCTCGCCCTCCAGCgccgcccg                                                                        >  2:4497046/1‑91 (MQ=255)
                 gAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCCTCATGCTCGCCTTCCCGCgccgcccg                                                                        >  2:2901105/1‑91 (MQ=255)
                 gAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCgccgcccg                                                                        >  2:4027833/1‑91 (MQ=255)
                  aaGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGccgcccgc                                                                       >  1:4257841/1‑93 (MQ=255)
                  aaGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGc                                                                       >  1:948863/1‑93 (MQ=255)
                   aGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCa                                                                      <  1:1025950/93‑1 (MQ=255)
                   aGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCCccacca                                                                      >  1:1062801/1‑93 (MQ=255)
                   aGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCAAGCTCGCCCTCCAGCGCCGCCAGCa                                                                      >  1:566112/1‑93 (MQ=255)
                     ttgttgAAGTTACATGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATc                                                                    <  2:3329201/93‑1 (MQ=255)
                       gttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGGATCAc                                                                  >  2:1486918/1‑93 (MQ=255)
                         tgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAcgc                                                                <  1:3798713/93‑1 (MQ=255)
                          taaGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAcgcg                                                               <  1:1854791/92‑1 (MQ=255)
                           aaGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCgg                                                              <  2:580474/93‑1 (MQ=255)
                           aaGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCgg                                                              <  2:2962884/93‑1 (MQ=255)
                            aGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCCGCATCACGCggg                                                             >  2:4103987/1‑93 (MQ=255)
                            aGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCCGCATCACGCggg                                                             >  2:4161762/1‑93 (MQ=255)
                             gTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCgggg                                                            >  2:2640233/1‑92 (MQ=255)
                                aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGTGaa                                                         <  2:3734354/93‑1 (MQ=255)
                                  cTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGt                                                       >  1:216478/1‑93 (MQ=255)
                                   tGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGtt                                                      <  2:4529570/93‑1 (MQ=255)
                                     ccGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGa                                                    <  2:2917454/93‑1 (MQ=255)
                                       gACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACCGGGGGGAAGTTGGTa                                                  >  2:3470728/1‑93 (MQ=255)
                                        aCAACGCTTTTTCTTGTGGCAGGAACACTTTTCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGGTGATaa                                                 >  2:2658706/1‑93 (MQ=255)
                                        aCAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATaa                                                 >  2:1840214/1‑93 (MQ=255)
                                        aCAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATaa                                                 <  1:1589794/93‑1 (MQ=255)
                                         cAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTGGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATaaa                                                <  1:3304682/93‑1 (MQ=255)
                                         cAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATaaa                                                <  2:408883/93‑1 (MQ=255)
                                         cAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATaaa                                                <  1:415725/93‑1 (MQ=255)
                                          aaCGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAg                                               <  1:2214410/93‑1 (MQ=255)
                                           aCGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGc                                              >  1:3775333/1‑93 (MQ=255)
                                           aCGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGc                                              <  1:3511191/93‑1 (MQ=255)
                                           aCGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGc                                              >  2:2266969/1‑93 (MQ=255)
                                            cGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCt                                             <  2:916776/93‑1 (MQ=255)
                                             gCTTTTTCTTGTGGCAGGCACACTTTGCCGCTGGTGGCCTCCTGCTCGCCTTCCAGCGCCCCCAGCATCACGCGGGCGGCGTTGATAAAGCTc                                            >  2:2184954/1‑93 (MQ=255)
                                              cTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCa                                           >  2:1306556/1‑93 (MQ=255)
                                              cTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCa                                           <  1:1319764/93‑1 (MQ=255)
                                              cTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCa                                           <  2:2983108/93‑1 (MQ=255)
                                              cTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCa                                           <  1:2817130/93‑1 (MQ=255)
                                              cTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGGGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCa                                           >  2:1406406/1‑93 (MQ=255)
                                               tttttCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCAt                                          <  1:23706/93‑1 (MQ=255)
                                               tttttCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCAt                                          <  1:2646264/93‑1 (MQ=255)
                                               tttttCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCAt                                          <  1:320207/93‑1 (MQ=255)
                                                ttttCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATg                                         <  1:455246/93‑1 (MQ=255)
                                                   tCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCg                                      <  1:675893/93‑1 (MQ=255)
                                                   tCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCg                                      <  1:2374796/93‑1 (MQ=255)
                                                     ttGTGGCAGGAACACTTTTCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGt                                    >  2:1389668/1‑93 (MQ=255)
                                                     ttGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAAGCGGGCGGAGTTGATAAAGCTCATGCCGGt                                    >  2:3434074/1‑93 (MQ=255)
                                                      tgtgGCAAGAACGCTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTa                                   <  1:715549/93‑1 (MQ=255)
                                                          tcAGGAACACTTTGCCGCTGGAGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGc                               <  1:2513551/92‑1 (MQ=255)
                                                          gcagcaACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGc                               <  2:1722695/93‑1 (MQ=255)
                                                          gCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGc                               <  1:1363912/93‑1 (MQ=255)
                                                          gCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGc                               <  1:2718387/93‑1 (MQ=255)
                                                            aGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGa                             >  1:2699549/1‑93 (MQ=255)
                                                            aGGAACACTTGGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGa                             <  2:3150606/93‑1 (MQ=255)
                                                             ggAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGAt                            >  1:1057206/1‑93 (MQ=255)
                                                             ggAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGAt                            >  2:2805853/1‑93 (MQ=255)
                                                             ggAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGAt                            >  1:400946/1‑93 (MQ=255)
                                                             ggAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGAt                            >  1:4046602/1‑93 (MQ=255)
                                                             ggAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGAt                            >  1:964346/1‑93 (MQ=255)
                                                             ggAACACTTACCAGCTGGTGGCATCATGCCCGCCTTCCAGCACCGCCATCCACACGCGTGAGAGCTGGATAAACCTCATGCCGGGATAGCGAt                            >  2:3607726/1‑93 (MQ=255)
                                                              gAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATa                           <  2:933283/93‑1 (MQ=255)
                                                              gAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATa                           <  1:1020696/93‑1 (MQ=255)
                                                               aaCACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAg                          >  1:776982/1‑93 (MQ=255)
                                                               aaCACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAg                          <  1:3877053/93‑1 (MQ=255)
                                                               aaCACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAg                          <  1:3772126/93‑1 (MQ=255)
                                                               aaCACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAg                          >  1:3545342/1‑93 (MQ=255)
                                                                acacTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGc                         >  2:2783614/1‑93 (MQ=255)
                                                                 cacTTTGCCGCTGGTGGCATCATGCTCGCCTTGCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGcc                        <  1:4488369/93‑1 (MQ=255)
                                                                 cacTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGcc                        <  1:1447599/93‑1 (MQ=255)
                                                                    tttGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCa                     >  1:4517924/1‑93 (MQ=255)
                                                                     ttGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCAt                    <  2:529625/93‑1 (MQ=255)
                                                                      tGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCAtt                   <  2:3987807/93‑1 (MQ=255)
                                                                       tccGCTGGTCGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCAttt                  <  2:569325/92‑1 (MQ=255)
                                                                        ccGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATTCCGGTACAGCGATAGCCCCATTTg                 >  1:4684616/1‑93 (MQ=255)
                                                                         cGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGc                <  1:3452043/93‑1 (MQ=255)
                                                                            tggtggCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCAc             >  2:4171486/1‑93 (MQ=255)
                                                                               tggCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGa          <  2:2851544/93‑1 (MQ=255)
                                                                                ggCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAc         >  2:1595661/1‑93 (MQ=255)
                                                                                ggCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAc         <  2:2597704/93‑1 (MQ=255)
                                                                                 gCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcg        >  1:182198/1‑93 (MQ=255)
                                                                                    tcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcggcg     >  1:4153839/1‑93 (MQ=255)
                                                                                    tcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcggcg     <  1:98962/93‑1 (MQ=255)
                                                                                      atGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGt   <  2:1125153/93‑1 (MQ=255)
                                                                                       tGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGtt  >  2:1421571/1‑93 (MQ=255)
                                                                                       tGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGtt  >  2:215719/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGTT  >  NC_000913/860322‑860501

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑