Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,231,608 G→C 8.2% Q141H (CAG→CAC)  umuC → translesion error‑prone DNA polymerase V subunit; DNA polymerase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,231,6080GC 8.2% 270.8 / 9.9 81Q141H (CAG→CAC) umuCtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base C (4/3);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGACGGCGGA  >  NC_000913/1231462‑1231756
                                                                                                                                                   |                                                                                                                                                    
ccccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGAc                                                                                                                                                    >  5:136041/1‑149 (MQ=255)
 cccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑TATCGCCCAGAc                                                                                                                                                    >  6:53141/1‑148 (MQ=255)
  ccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGAcc                                                                                                                                                   >  5:76340/1‑148 (MQ=255)
     cgTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATC‑CCCACACCaaca                                                                                                                                               <  3:171118/149‑3 (MQ=255)
       tCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACg                                                                                                                                             >  3:87088/1‑149 (MQ=255)
          agaTTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTg                                                                                                                                          <  7:95693/149‑1 (MQ=255)
           gaTTTACAGTATTGATGAGGCTTTCTGCGATCTGACAGGTGTGCGTTATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTg                                                                                                                                          <  6:365762/148‑1 (MQ=255)
               tACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATC‑CCCACACCaaca                                                                                                                                               <  7:109343/139‑3 (MQ=255)
               tACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATC‑CCCACACCaaca                                                                                                                                               >  8:109343/1‑137 (MQ=255)
                aCAGTATTGATGCGGCATTCTGCGATCTGACAGGGGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTaa                                                                                                                                     <  1:273691/148‑1 (MQ=255)
                   gTATTGATGAGGCATTCTGCGAGCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCtt                                                                                                                                 <  8:212420/149‑1 (MQ=255)
                            aGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                                         <  6:76340/148‑1 (MQ=255)
                            aGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                                         <  2:250418/148‑1 (MQ=255)
                            aGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                                         <  8:182337/148‑1 (MQ=255)
                             ggCATTCTGCGATCTGACATGTGTGCGTAATTGTCGCGATCTGACTGATTTTGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATg                                                                                                                       <  8:1206/149‑1 (MQ=255)
                             ggCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATC‑CCCACACCaacag                                                                                                                                              <  4:55894/126‑4 (MQ=255)
                             ggCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATC‑CCCACACCaacag                                                                                                                                              >  3:55894/1‑123 (MQ=255)
                              gCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACATCGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGc                                                                                                                      <  2:346421/149‑1 (MQ=255)
                              gCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGc                                                                                                                      >  3:21940/1‑149 (MQ=255)
                               cATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACTGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCg                                                                                                                     <  1:253696/149‑1 (MQ=255)
                                 ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATc                                                                                                                          <  3:59393/142‑1 (MQ=255)
                                 ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATc                                                                                                                          >  4:59393/1‑142 (MQ=255)
                                 ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                                         >  4:81772/1‑143 (MQ=255)
                                 ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                                         <  3:81772/143‑1 (MQ=255)
                                   cTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCa                                                                                                                  >  2:271203/1‑148 (MQ=255)
                                    tGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaa                                                                                                                 >  6:207271/1‑148 (MQ=255)
                                        aTCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                                                                                            >  2:300577/1‑149 (MQ=255)
                                             aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                                                                                            <  8:339729/144‑1 (MQ=255)
                                             aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                                                                                            <  3:266965/144‑1 (MQ=255)
                                             aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAAt                                                                                                           <  6:323158/145‑1 (MQ=255)
                                             aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATg                                                                                                          <  5:207271/146‑1 (MQ=255)
                                             aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCa                                                                                                       <  1:300577/149‑1 (MQ=255)
                                              cAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAg                                                                                                      <  1:403380/149‑1 (MQ=255)
                                               aGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGc                                                                                                     <  8:219964/149‑1 (MQ=255)
                                                    tgCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                                         >  6:344899/1‑124 (MQ=255)
                                                    tgCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                                         <  5:344899/124‑1 (MQ=255)
                                                    tgCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAg                                                                                                >  5:271007/1‑149 (MQ=255)
                                                      cGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGAc                                                                                              >  8:318620/1‑149 (MQ=255)
                                                        tAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgg                                                                                            >  8:353210/1‑149 (MQ=255)
                                                          aTTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggt                                                                                          <  8:214804/149‑1 (MQ=255)
                                                           ttGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggtg                                                                                         <  2:189970/149‑1 (MQ=255)
                                                              tCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggtggg                                                                                       <  7:44717/148‑1 (MQ=255)
                                                              tCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggtggg                                                                                       <  4:21940/148‑1 (MQ=255)
                                                                gcgATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCaa                                                                                                                                                 >  7:352234/1‑88 (MQ=255)
                                                                gcgATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCaa                                                                                                                                                 <  8:352234/88‑1 (MQ=255)
                                                                      tgactgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCACACAAAAACGCAGGCTAAGCTACCCAATCATGCGGCAAAAAAATGGAAGCGGCATACGGCTGGGGATGTGGAt                                                                              >  3:44575/1‑149 (MQ=255)
                                                                        actgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAttt                                                                            <  4:80409/149‑1 (MQ=255)
                                                                        actgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGGTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGa                                                                               <  6:134577/146‑1 (MQ=255)
                                                                                 ggCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑TATCGCCCAGACCAAAACGCTGGCTAAGGTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATc                                                                    <  5:53141/148‑1 (MQ=255)
                                                                                 ggCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCAAAGCTTGCCAATAATGCGGCAAAAAAATGGCATCGGCAGACGGGTGGGTTGGTGGATTTATCAAAtat                                                                   >  6:280627/1‑147 (MQ=255)
                                                                                  gCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCt                                                                   <  6:136041/148‑1 (MQ=255)
                                                                                   cAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTg                                                                  >  2:403057/1‑148 (MQ=255)
                                                                                    agagAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGggtggt                                                                                  <  1:75700/131‑1 (MQ=255)
                                                                                    agagAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGggtggt                                                                                  >  2:75700/1‑131 (MQ=255)
                                                                                     gagaAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGaa                                                               >  3:10828/1‑149 (MQ=255)
                                                                                          aTTCGCGCAACGGTGCTACATGGTACGCATCTTAATGTTTGTTTGGG‑GATCGCACAGACCAAATCGCGTGCTATGATTGCCATTCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTTGTGGATTTATGATATCTGGAACGCCa                                                          <  4:44575/149‑1 (MQ=255)
                                                                                          aTTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGCCCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGc                                                            >  8:336709/1‑147 (MQ=255)
                                                                                                gcAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTa                                                     >  2:50471/1‑148 (MQ=255)
                                                                                                   aCGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGAAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAAt                                                  >  1:294980/1‑148 (MQ=255)
                                                                                                     ggTGCTACCACGTATCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTaa                                               <  5:183375/149‑1 (MQ=255)
                                                                                                      gTGCTACAACGAACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAAt                                              <  3:9531/149‑1 (MQ=255)
                                                                                                       tGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATg                                             >  6:403056/1‑149 (MQ=255)
                                                                                                                      aTCTTACAGTTGGTCTGGG‑GATATACCACACCAAACCTCTGGCTAAGCTTGAAAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGCGAGGTGTATTTATCAAATCTCTAACGCCAGCGTAAATTAATTTATGCTCACCCAGTg                              >  2:279478/1‑149 (MQ=255)
                                                                                                                       tcataatGTTGGTGTGGC‑GATCGCCCAGACCAAAACGCAGGCTAAGCTTGCCAATCATTCGGCAAAAAAATGGCAGCGGCAGACGGGTGTGGTGGTGGATTTATCAAATCTGTAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTgg                             <  2:294980/143‑1 (MQ=255)
                                                                                                                       tCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTg                              >  7:54113/1‑148 (MQ=255)
                                                                                                                       tCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTg                              >  7:71572/1‑148 (MQ=255)
                                                                                                                           aCTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGa                            <  7:353210/146‑1 (MQ=255)
                                                                                                                            cTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGa                         <  8:71572/148‑1 (MQ=255)
                                                                                                                               ttGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTc                     >  7:64604/1‑149 (MQ=255)
                                                                                                                               ttGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTc                     <  7:318620/149‑1 (MQ=255)
                                                                                                                                   tgtgGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAtt                                                                             >  1:40264/1‑89 (MQ=255)
                                                                                                                                   tgtgGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAtt                                                                             <  2:40264/89‑1 (MQ=255)
                                                                                                                                    gtgGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAACATGGCAGCGGCAGCCGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCGCCGTGGAGGCCGTCTgggg                >  4:145798/1‑149 (MQ=255)
                                                                                                                                    gtgGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTgggg                >  1:401229/1‑149 (MQ=255)
                                                                                                                                     tgGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCAGTGGATGACGTCTgggg                >  8:98410/1‑148 (MQ=255)
                                                                                                                                     tgGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTATGCTCTCCCCGTGGATGACGTCTgggg                <  3:140273/148‑1 (MQ=255)
                                                                                                                                     tgGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAAGGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTgggg                >  1:110934/1‑148 (MQ=255)
                                                                                                                                          gATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTg           >  8:390836/1‑149 (MQ=255)
                                                                                                                                               cccAGCCCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGAcgg      >  5:197621/1‑149 (MQ=255)
                                                                                                                                                 cAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGC‑AAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGAt             >  4:216890/1‑139 (MQ=255)
                                                                                                                                                 cAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGC‑AAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGAt             <  3:216890/139‑1 (MQ=255)
                                                                                                                                                   gACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGACGGCGGa  <  1:403057/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGG‑GATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGACGGCGGA  >  NC_000913/1231462‑1231756

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