Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,484 A→G 5.8% intergenic (+30/+16) 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,4840AG 5.8% 187.5 / 2.5 70intergenic (+30/+16)umuC/dsbBDNA polymerase V catalytic protein/protein thiol:quinone oxidoreductase DsbB
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base G (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.43e-01

GAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCC  >  NC_000913/1232345‑1232626
                                                                                                                                           |                                                                                                                                              
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)
                                                                                                                                           aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGCACTCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATGCCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATcc  >  1:271896/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCC  >  NC_000913/1232345‑1232626

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: