Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 5.9% intergenic (+27/+19) umuC → / ← dsbB DNA polymerase V catalytic protein/protein thiol:quinone oxidoreductase DsbB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,4810CT 5.9% 183.1 / 2.3 68intergenic (+27/+19)umuC/dsbBDNA polymerase V catalytic protein/protein thiol:quinone oxidoreductase DsbB
Reads supporting (aligned to +/- strand):  ref base C (36/28);  new base T (2/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

GCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAA  >  NC_000913/1232339‑1232622
                                                                                                                                              |                                                                                                                                             
gCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCCTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                                                               >  2:992249/1‑143 (MQ=255)
  agagGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGc                                                                                                                                             >  1:810399/1‑143 (MQ=255)
      gAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGACGTTGATCGATACAGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGa                                                                                                                                         >  1:241959/1‑143 (MQ=255)
         cacTCTATTTTGCCGGGCAGGGGTTCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                                                                      <  1:21755/143‑1 (MQ=255)
         cacTCTATTTGGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                                                                      <  2:911400/143‑1 (MQ=255)
           ctctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                                                                    <  1:1041370/143‑1 (MQ=255)
           accgaTTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgct                                                                                                                                    <  1:530727/139‑2 (MQ=255)
            tctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTAGACAACGCGACGCTATGATCTACTGAGGGGCAAATAAATATAGCGGCAGGAAACAAGCGCTCACGCAGGAGGGCCg                                                                                                                                   >  2:372761/1‑143 (MQ=255)
              tATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGaa                                                                                                                                 >  1:354736/1‑143 (MQ=255)
                 tttGCCGGGCGGGGGATCCAGCAGCAAGGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGGAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttc                                                                                                                              <  2:921870/143‑8 (MQ=255)
                  ttGCCGGGAAGTGGATGCAGCGAGAATGGCAGATGAAGCGAGCCATGCTATCACCACGTTATACAACGCGAAGGTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGa                                                                                                                             <  2:465855/143‑1 (MQ=255)
                  gtGCCGGGCAGGGGATTCAGATAGTATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGa                                                                                                                             <  2:4672/142‑1 (MQ=255)
                     ccGGGCAGGGGATTCAGGAACAATGGCAGACGTAGCGAGCCATGGTTTCACCGCGGCAGACAACGCGAAGGTCTCATTTGCGGAGGGTGAAAGAAATAGAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                                                          <  2:472190/143‑1 (MQ=255)
                     ccGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                                                          >  1:797439/1‑143 (MQ=255)
                     ccGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                                                          <  1:416171/143‑1 (MQ=255)
                     ccGGGCAGGGGAACCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGAGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                                                          <  1:204527/143‑1 (MQ=255)
                      cGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCCGGAAAAAAGCGCCCCCGCAGGAGCGCCGAATGGACTAg                                                                                                                         >  2:182655/1‑143 (MQ=255)
                           aGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGGAGGAAAAACGCGCTGCCGCAGGACCGCCGAATGGATTAGCgcac                                                                                                                    >  1:1073934/1‑140 (MQ=255)
                           aGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAcc                                                                                                                    >  1:152061/1‑143 (MQ=255)
                            ggggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                                                                   >  1:676202/1‑143 (MQ=255)
                                   cAGCAACCATGGGAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTGACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                                                            <  1:414588/143‑1 (MQ=255)
                                   cAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                                                            >  1:1003033/1‑143 (MQ=255)
                                   cAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTAGCGCATCATTGCCGAATGGAGTAGCGACCGAACAGAt                                                                                                            >  1:830333/1‑143 (MQ=255)
                                    aGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGGGGGTCAAATAAATATAGCGGGAGGAAAAAAGCGCACCCGCAGGAGCGCCGAATCAACAAGGGACCGAAcggggc                                                                                                           >  2:569902/1‑138 (MQ=255)
                                    aGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                                                           >  2:349684/1‑143 (MQ=255)
                                          aaTGGCAGATGAAGCGAGCCATGCTGTCACCACGTTATACAACGCGAAGTTCGGATTAACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttt                                                                                                     <  1:39186/143‑1 (MQ=255)
                                            tGGCAGATGAAGCGAGCCATGCTTCCACCAAGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttttt                                                                                                   <  1:639072/143‑1 (MQ=255)
                                                aGATGAAGCGAGCCATGCTTTCACCACGTGATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                                                               <  1:764034/143‑1 (MQ=255)
                                                 gATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCt                                                                                              >  2:151291/1‑143 (MQ=255)
                                                 gATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCt                                                                                              <  2:478797/143‑1 (MQ=255)
                                                   tGAAGCGAGCCATGCTTTCACCACGATATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCttt                                                                                            >  1:847001/1‑143 (MQ=255)
                                                       gCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGCTCACGTTTTTTCGCTTTAAAc                                                                                        >  2:952360/1‑143 (MQ=255)
                                                             cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGGGCAGATCCCGTTTTGTCGCTTTAAACGGATgg                                                                                  >  1:236100/1‑143 (MQ=255)
                                                             cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                                                                  <  1:749670/143‑1 (MQ=255)
                                                              aTGCTTGCACCACGTTATACAACGCGAAGTTCTGCTTCACTGAGGGCCAAATAAATATAGCGGCAGGACACAAGCGCTCCCGCAGGAGCGCCGAAGGGATGAGCGACCGACCAGATCACGTTTTTTCGCTTTAAACGGccagc                                                                                 >  2:234176/1‑139 (MQ=255)
                                                                  tttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAt                                                                             <  1:832723/143‑1 (MQ=255)
                                                                    tctcgacGTGATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                                                           <  2:941986/138‑1 (MQ=255)
                                                                        cacGTTATACAACGCGAAGTTCTGAGTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcaccac                                                                       <  2:817174/143‑1 (MQ=255)
                                                                                aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTg                                                               <  1:125476/143‑1 (MQ=255)
                                                                                  aaCGCGAAGTTCTGCTTTACTGAGGGTCAAATAAATCTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                                                             <  1:74617/143‑1 (MQ=255)
                                                                                  aaCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCCGGGAAAACAGCACCAGCACTGCg                                                             >  2:45865/1‑143 (MQ=255)
                                                                                    cgcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACCGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCCCCAGCACGGCGAc                                                           >  1:646331/1‑143 (MQ=255)
                                                                                           ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAGTTGATCACGTTTTCTCCCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAGCAgg                                                    >  1:198686/1‑143 (MQ=255)
                                                                                            tCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAACACCACCAGCACTGCGACCATCAggg                                                   >  2:606213/1‑142 (MQ=255)
                                                                                            tCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAGCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACATCCTGCACTGCCACAATCAggg                                                   >  2:507787/1‑142 (MQ=255)
                                                                                                aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCCCTGCGACAATCAGGTAAGc                                               >  1:681936/1‑143 (MQ=255)
                                                                                                       gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAACCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaaa                                        >  1:15116/1‑143 (MQ=255)
                                                                                                           gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAGTGGATTAGCGACCGAACAGGTCACGTTTTTTCGCTTTAAACGGCTGGGAGATCACCGCCAGCACTGCGACGATCAGGTAAGAGATAAAAATAc                                    >  1:848597/1‑143 (MQ=255)
                                                                                                           gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGTGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAACACTGCGACAATCAGGTAAGCGATAAAAATAc                                    >  2:531847/1‑143 (MQ=255)
                                                                                                           gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAc                                    >  2:43735/1‑143 (MQ=255)
                                                                                                            tCAAAGAAATACAGTGGCATGAAAAGAGCGCTCCCGCAAGAGCGCCGAATGGATTAGTGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                   <  2:119490/143‑1 (MQ=255)
                                                                                                              aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAGCACCACAAGCACTGGGACAATCAGGTAAGAGATAAAAATACCGa                                 >  1:496971/1‑143 (MQ=255)
                                                                                                              aaataaacATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGa                                 <  1:221928/143‑1 (MQ=255)
                                                                                                               aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGag                                <  1:992249/143‑1 (MQ=255)
                                                                                                                 taaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagca                              >  1:210198/1‑143 (MQ=255)
                                                                                                                 taaataGAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagca                              <  1:922830/143‑1 (MQ=255)
                                                                                                                      ataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACGAGCACTGCGACAATCAGGTAAGCGATGGAAATACCGAGCAGCCAc                         >  1:280339/1‑143 (MQ=255)
                                                                                                                          cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCg                     <  1:349684/143‑1 (MQ=255)
                                                                                                                            gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATCGATTAGCCACCGAACAGAACAAGAATATTAAATATAAACGGCTGGGAAATCACCACCAGCACTGCAACACTCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGc                   >  1:813924/1‑143 (MQ=255)
                                                                                                                               ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCAtt                <  1:747335/143‑1 (MQ=255)
                                                                                                                                gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTATAAACGACTGGGAAATCACCACCAGCACTGCGACAATCAGGTCAGCGATAACAAGACCGAACAGCCACTCCGGCAatt               >  2:778247/1‑143 (MQ=255)
                                                                                                                                   aaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCa            <  1:942010/143‑1 (MQ=255)
                                                                                                                                   aaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGACTGGGAAATCACCACCAGCACTGCGACCATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCa            >  1:916260/1‑143 (MQ=255)
                                                                                                                                    tcaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg           >  1:766318/3‑143 (MQ=255)
                                                                                                                                      ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTACACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGGACCCGCTAAACATACCGAGCAGCCACTGCGGCATTTCCAGAc         >  1:1062270/2‑143 (MQ=255)
                                                                                                                                      ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCCCCAGCCCTGCGACAATCAGGGAAGCGATAAAAAAACCAACCAGCCACTGCGGCATGTCCacac         >  2:479430/2‑143 (MQ=255)
                                                                                                                                      ctcgcTCCCGGCGGGGCGCCGAATGGAGTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAc         <  2:857337/141‑1 (MQ=255)
                                                                                                                                      aGCGCGCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAc         <  1:609257/143‑1 (MQ=255)
                                                                                                                                           tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGCCACTCGGGTAAGCGATAAAAAGACCGAGCAGCCACTGCGGCATTTACGGACCTaaa    >  1:143099/1‑143 (MQ=255)
                                                                                                                                           gcccGCAGGAGCGCCTAATGGATTAGCGCCCGACCAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaa    <  2:113858/142‑1 (MQ=255)
                                                                                                                                             ccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCCCCAGCACTGCGCCAATAAGGTAAACGATAAAAAAAACAAACAGCCACTGCGGAATTTCCAGAACTaaaaa  >  1:312897/1‑143 (MQ=255)
                                                                                                                                             ccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCGGCCCTacaaa  >  1:148889/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAA  >  NC_000913/1232339‑1232622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: