Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,231,594 G→C 7.3% V137L (GTG→CTG)  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,5940GC 7.3% 261.8 / 9.1 79V137L (GTG→CTG) umuCtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base C (3/3);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGG  >  NC_000913/1231448‑1231742
                                                                                                                                                  |                                                                                                                                                     
tGGAAGAGCTATCGCCCCGCGTCGAGAATTCCAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCTGAATTGTCGCGATCTGACTGATTTTGTCAGAGAAATTCGAGAAACTGTGCTAGATCGTACCCATCATACTTTTGTtgtg                                                                                                                                                     >  1:280276/1‑149 (MQ=255)
  gAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTggg                                                                                                                                                   <  7:279736/149‑1 (MQ=255)
  gAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTggg                                                                                                                                                   >  7:219964/1‑149 (MQ=255)
  gAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTgg                                                                                                                                                    <  5:192872/148‑1 (MQ=255)
     gagCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGAt                                                                                                                                               <  1:46389/149‑1 (MQ=255)
     gagCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGAt                                                                                                                                               >  2:46389/1‑149 (MQ=255)
         tATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGcc                                                                                                                                           >  7:214804/1‑149 (MQ=255)
          aTCGCCACGCGACGAGATTTACACTATTGATGAGCCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTCGTGTGGGGATCGcc                                                                                                                                           >  2:403380/1‑148 (MQ=255)
           tCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGcc                                                                                                                                           >  3:80409/1‑147 (MQ=255)
           gcGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTTCGATCTGCCAGGGGGGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACACCGTACCCATCTTACTGTTGGTGTGGGGATCGCCCa                                                                                                                                         <  8:80841/148‑1 (MQ=255)
              ccccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGAc                                                                                                                                      >  5:136041/1‑149 (MQ=255)
               cccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGTATCGCCCAGAc                                                                                                                                      >  6:53141/1‑148 (MQ=255)
                ccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGgtgt                                                                                                                                                      >  1:84784/1‑132 (MQ=255)
                ccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGgtgt                                                                                                                                                      <  2:84784/132‑1 (MQ=255)
                ccGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGAcc                                                                                                                                     >  5:76340/1‑148 (MQ=255)
                   cgTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATCCCCACACCaaca                                                                                                                                 <  3:171118/149‑3 (MQ=255)
                     tCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACg                                                                                                                               >  3:87088/1‑149 (MQ=255)
                        agaTTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTg                                                                                                                            <  7:95693/149‑1 (MQ=255)
                         gaTTTACAGTATTGATGAGGCTTTCTGCGATCTGACAGGTGTGCGTTATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTg                                                                                                                            <  6:365762/148‑1 (MQ=255)
                             tACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATCCCCACACCaaca                                                                                                                                 >  8:109343/1‑137 (MQ=255)
                             tACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATCCCCACACCaaca                                                                                                                                 <  7:109343/139‑3 (MQ=255)
                              aCAGTATTGATGCGGCATTCTGCGATCTGACAGGGGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTaa                                                                                                                       <  1:273691/148‑1 (MQ=255)
                                 gTATTGATGAGGCATTCTGCGAGCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCtt                                                                                                                   <  8:212420/149‑1 (MQ=255)
                                          aGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                           <  2:250418/148‑1 (MQ=255)
                                          aGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                           <  8:182337/148‑1 (MQ=255)
                                          aGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                           <  6:76340/148‑1 (MQ=255)
                                           ggCATTCTGCGATCTGACATGTGTGCGTAATTGTCGCGATCTGACTGATTTTGTCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATg                                                                                                         <  8:1206/149‑1 (MQ=255)
                                           ggCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATCCCCACACCaacag                                                                                                                                >  3:55894/1‑123 (MQ=255)
                                           ggCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTCTGGGCGATCCCCACACCaacag                                                                                                                                <  4:55894/126‑4 (MQ=255)
                                            gCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACATCGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGc                                                                                                        <  2:346421/149‑1 (MQ=255)
                                            gCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGc                                                                                                        >  3:21940/1‑149 (MQ=255)
                                             cATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACTGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCg                                                                                                       <  1:253696/149‑1 (MQ=255)
                                               ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATc                                                                                                            <  3:59393/142‑1 (MQ=255)
                                               ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATc                                                                                                            >  4:59393/1‑142 (MQ=255)
                                               ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                           >  4:81772/1‑143 (MQ=255)
                                               ttCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                           <  3:81772/143‑1 (MQ=255)
                                                 cTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCa                                                                                                    >  2:271203/1‑148 (MQ=255)
                                                  tGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaa                                                                                                   >  6:207271/1‑148 (MQ=255)
                                                      aTCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                                                                              >  2:300577/1‑149 (MQ=255)
                                                           aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                                                                              <  8:339729/144‑1 (MQ=255)
                                                           aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCaaaaaaa                                                                                              <  3:266965/144‑1 (MQ=255)
                                                           aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAAt                                                                                             <  6:323158/145‑1 (MQ=255)
                                                           aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATg                                                                                            <  5:207271/146‑1 (MQ=255)
                                                           aCAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCa                                                                                         <  1:300577/149‑1 (MQ=255)
                                                            cAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAg                                                                                        <  1:403380/149‑1 (MQ=255)
                                                             aGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGc                                                                                       <  8:219964/149‑1 (MQ=255)
                                                                  tgCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                           >  6:344899/1‑124 (MQ=255)
                                                                  tgCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCa                                                                                                           <  5:344899/124‑1 (MQ=255)
                                                                  tgCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAg                                                                                  >  5:271007/1‑149 (MQ=255)
                                                                    cGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGAc                                                                                >  8:318620/1‑149 (MQ=255)
                                                                      tAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgg                                                                              >  8:353210/1‑149 (MQ=255)
                                                                        aTTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggt                                                                            <  8:214804/149‑1 (MQ=255)
                                                                         ttGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggtg                                                                           <  2:189970/149‑1 (MQ=255)
                                                                            tCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggtggg                                                                         <  4:21940/148‑1 (MQ=255)
                                                                            tCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACgggtggg                                                                         <  7:44717/148‑1 (MQ=255)
                                                                              gcgATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCaa                                                                                                                                   <  8:352234/88‑1 (MQ=255)
                                                                              gcgATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCaa                                                                                                                                   >  7:352234/1‑88 (MQ=255)
                                                                                    tgactgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCACACAAAAACGCAGGCTAAGCTACCCAATCATGCGGCAAAAAAATGGAAGCGGCATACGGCTGGGGATGTGGAt                                                                >  3:44575/1‑149 (MQ=255)
                                                                                      actgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAttt                                                              <  4:80409/149‑1 (MQ=255)
                                                                                      actgaTTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGGTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGa                                                                 <  6:134577/146‑1 (MQ=255)
                                                                                               ggCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGTATCGCCCAGACCAAAACGCTGGCTAAGGTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATc                                                      <  5:53141/148‑1 (MQ=255)
                                                                                               ggCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCAAAGCTTGCCAATAATGCGGCAAAAAAATGGCATCGGCAGACGGGTGGGTTGGTGGATTTATCAAAtat                                                     >  6:280627/1‑147 (MQ=255)
                                                                                                gCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCt                                                     <  6:136041/148‑1 (MQ=255)
                                                                                                 cAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTg                                                    >  2:403057/1‑148 (MQ=255)
                                                                                                  agagAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGggtggt                                                                    >  2:75700/1‑131 (MQ=255)
                                                                                                  agagAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGggtggt                                                                    <  1:75700/131‑1 (MQ=255)
                                                                                                   gagaAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGaa                                                 >  3:10828/1‑149 (MQ=255)
                                                                                                        aTTCGCGCAACGGTGCTACATGGTACGCATCTTAATGTTTGTTTGGGGATCGCACAGACCAAATCGCGTGCTATGATTGCCATTCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTTGTGGATTTATGATATCTGGAACGCCa                                            <  4:44575/149‑1 (MQ=255)
                                                                                                        aTTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGCCCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGc                                              >  8:336709/1‑147 (MQ=255)
                                                                                                              gcAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTa                                       >  2:50471/1‑148 (MQ=255)
                                                                                                                 aCGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGAAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAAt                                    >  1:294980/1‑148 (MQ=255)
                                                                                                                   ggTGCTACCACGTATCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTaa                                 <  5:183375/149‑1 (MQ=255)
                                                                                                                    gTGCTACAACGAACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAAt                                <  3:9531/149‑1 (MQ=255)
                                                                                                                     tGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATg                               >  6:403056/1‑149 (MQ=255)
                                                                                                                                    aTCTTACAGTTGGTCTGGGGATATACCACACCAAACCTCTGGCTAAGCTTGAAAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGCGAGGTGTATTTATCAAATCTCTAACGCCAGCGTAAATTAATTTATGCTCACCCAGTg                >  2:279478/1‑149 (MQ=255)
                                                                                                                                     tcataatGTTGGTGTGGCGATCGCCCAGACCAAAACGCAGGCTAAGCTTGCCAATCATTCGGCAAAAAAATGGCAGCGGCAGACGGGTGTGGTGGTGGATTTATCAAATCTGTAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTgg               <  2:294980/143‑1 (MQ=255)
                                                                                                                                     tCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTg                >  7:71572/1‑148 (MQ=255)
                                                                                                                                     tCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTg                >  7:54113/1‑148 (MQ=255)
                                                                                                                                         aCTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGa              <  7:353210/146‑1 (MQ=255)
                                                                                                                                          cTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGa           <  8:71572/148‑1 (MQ=255)
                                                                                                                                             ttGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTc       >  7:64604/1‑149 (MQ=255)
                                                                                                                                             ttGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTc       <  7:318620/149‑1 (MQ=255)
                                                                                                                                                 tgtgGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAtt                                                               <  2:40264/89‑1 (MQ=255)
                                                                                                                                                 tgtgGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGAtt                                                               >  1:40264/1‑89 (MQ=255)
                                                                                                                                                  gtgGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAACATGGCAGCGGCAGCCGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCGCCGTGGAGGCCGTCTgggg  >  4:145798/1‑149 (MQ=255)
                                                                                                                                                  gtgGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTgggg  >  1:401229/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGGAAGAGCTATCGCCCCGCGTCGAGATTTACAGTATTGATGAGGCATTCTGCGATCTGACAGGTGTGCGTAATTGTCGCGATCTGACTGATTTTGGCAGAGAAATTCGCGCAACGGTGCTACAACGTACCCATCTTACTGTTGGTGTGGGGATCGCCCAGACCAAAACGCTGGCTAAGCTTGCCAATCATGCGGCAAAAAAATGGCAGCGGCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGG  >  NC_000913/1231448‑1231742

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: