Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 790,931 C→G 9.5% W33C (TGG→TGC)  galT ← galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913790,9310CG 9.5% 211.4 / 11.2 63W33C (TGG→TGC) galTgalactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base G (4/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCA  >  NC_000913/790783‑791074
                                                                                                                                                    |                                                                                                                                                
tCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGc                                                                                                                                                  <  4:213713/149‑1 (MQ=255)
 cAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGcc                                                                                                                                                 <  1:216935/149‑1 (MQ=255)
  aaaGCCGCGAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCACCCTGCCa                                                                                                                                                <  2:666161/149‑1 (MQ=255)
  aaaGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCa                                                                                                                                                >  2:294586/1‑149 (MQ=255)
  aaaGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCa                                                                                                                                                >  3:243233/1‑149 (MQ=255)
  aaaGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCa                                                                                                                                                >  3:21557/1‑149 (MQ=255)
     gCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑ggg                                                                                                                                            <  1:652516/149‑1 (MQ=255)
                aTTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCG‑CCCCTGGCA‑gggggcgca                                                                                                                                      <  2:378411/143‑6 (MQ=255)
                aTTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCG‑CCCCTGGCA‑gggggcgca                                                                                                                                      >  1:378411/1‑138 (MQ=255)
                 ttAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGcgcg                                                                                                                                <  4:589411/149‑1 (MQ=255)
                     tGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCgcggtgc                                                                                                                            <  2:672118/149‑1 (MQ=255)
                     tGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCgcggtgc                                                                                                                            >  2:426263/1‑149 (MQ=255)
                      gAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCgcggtgcg                                                                                                                           >  1:189365/1‑149 (MQ=255)
                       aaaaCGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCgcggtgcgg                                                                                                                          <  2:116950/149‑1 (MQ=255)
                       aaaaCGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCgcggtgcgg                                                                                                                          <  1:179950/149‑1 (MQ=255)
                          aCGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGa                                                                                                                       >  1:193748/1‑149 (MQ=255)
                            gTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCG‑CCCCTGGCA‑gggg                                                                                                                                           <  2:637680/126‑1 (MQ=255)
                            gTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCG‑CCCCTGGCA‑gggg                                                                                                                                           >  1:637680/1‑126 (MQ=255)
                                 tCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCgcggtgcgg                                                                                                                          <  1:103004/139‑1 (MQ=255)
                                 tCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCgcggtgcgg                                                                                                                          >  2:103004/1‑139 (MQ=255)
                                 tCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGa                                                                                                                <  3:668443/149‑1 (MQ=255)
                                      gtgtAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCa                                                                                                           >  2:723438/1‑149 (MQ=255)
                                           aTCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGcc                                                                                                      >  2:43709/1‑149 (MQ=255)
                                            tCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGccc                                                                                                     <  4:344290/149‑1 (MQ=255)
                                            tCGGGGTTTTTACCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGccc                                                                                                     >  4:380669/1‑149 (MQ=255)
                                                  tttttATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGCATCCATTGCCCGGTGAg                                                                                               >  2:265291/1‑149 (MQ=255)
                                                        tCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGtt                                                                                         >  4:719302/1‑149 (MQ=255)
                                                              gTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAgcg                                                                                   >  2:250171/1‑149 (MQ=255)
                                                                   ccGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAt                                                                              >  4:752164/1‑149 (MQ=255)
                                                                    cGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAtg                                                                             <  4:594073/149‑1 (MQ=255)
                                                                       acaTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTgg                                                                          <  4:34627/149‑1 (MQ=255)
                                                                         aTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGgat                                                                        >  2:562546/1‑149 (MQ=255)
                                                                             ccTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATc                                                                    >  1:663678/1‑149 (MQ=255)
                                                                                  gcAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACgg                                                               >  2:745843/1‑149 (MQ=255)
                                                                                     gagGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGAt                                                            <  4:243233/149‑1 (MQ=255)
                                                                                        gAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTaa                                                         >  4:525266/1‑149 (MQ=255)
                                                                                        gAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTaa                                                         >  2:194515/1‑149 (MQ=255)
                                                                                                          cgcAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAAAGGGATTAAATTGCGTCATGGTCGTTc                                       <  2:663678/149‑1 (MQ=255)
                                                                                                           gcAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTcc                                      >  2:306772/1‑149 (MQ=255)
                                                                                                            cAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCt                                     <  1:482292/149‑1 (MQ=255)
                                                                                                             aGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCtt                                    >  3:157855/1‑149 (MQ=255)
                                                                                                               gTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTaa                                  <  2:88577/149‑1 (MQ=255)
                                                                                                                  acacCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCg                               <  2:671347/149‑1 (MQ=255)
                                                                                                                   cacCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCgg                              <  1:294586/149‑1 (MQ=255)
                                                                                                                    acCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCggg                             >  3:90964/1‑149 (MQ=255)
                                                                                                                    acCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCggg                             >  4:750455/1‑149 (MQ=255)
                                                                                                                     ccTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGa                            <  1:562546/149‑1 (MQ=255)
                                                                                                                              ctggcGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCtgt                   <  4:322825/149‑1 (MQ=255)
                                                                                                                                   gTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATg              >  4:604558/1‑149 (MQ=255)
                                                                                                                                    tttCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATgg             >  4:119787/1‑149 (MQ=255)
                                                                                                                                    tttCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATgg             <  2:249453/149‑1 (MQ=255)
                                                                                                                                     ttCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGc            <  3:719302/149‑1 (MQ=255)
                                                                                                                                      tCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCg           <  4:492783/149‑1 (MQ=255)
                                                                                                                                       ccTGCGCCCCCTGCCA‑GGGGCGCTTAGCGAGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGt          <  4:90964/149‑1 (MQ=255)
                                                                                                                                        taagcgccccTGGCAGGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGt          >  2:206258/7‑149 (MQ=255)
                                                                                                                                            cgccccTGGCAGGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACt      >  4:433845/3‑149 (MQ=255)
                                                                                                                                             cccccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgc    >  4:472784/1‑149 (MQ=255)
                                                                                                                                             cccccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgc    <  1:603962/149‑1 (MQ=255)
                                                                                                                                              ccccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcc   >  3:607220/1‑149 (MQ=255)
                                                                                                                                              ccccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcc   >  2:510997/1‑149 (MQ=255)
                                                                                                                                               cccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTGATCGGGATATCCCTGTGGATGGCGTGACtgcca  >  1:550102/1‑149 (MQ=255)
                                                                                                                                               cccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcca  <  1:426263/149‑1 (MQ=255)
                                                                                                                                               cccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcca  >  2:470352/1‑149 (MQ=255)
                                                                                                                                               cccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcca  <  2:367567/149‑1 (MQ=255)
                                                                                                                                               cccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcca  <  1:723438/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCA  >  NC_000913/790783‑791074

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