Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,754,874 G→A 5.5% intergenic (‑446/+99) yiaW ← / ← aldB DUF3302 family inner membrane protein/aldehyde dehydrogenase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,754,8740GA 5.5% 790.7 / 22.7 234intergenic (‑446/+99)yiaW/aldBDUF3302 family inner membrane protein/aldehyde dehydrogenase B
Reads supporting (aligned to +/- strand):  ref base G (112/109);  new base A (7/6);  total (119/115)
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.92e-01

Only 100 of 242 total aligned reads displayed.

GCCCGCTCACCATCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTG  >  NC_000913/3754726‑3755021
                                                                                                                                                    |                                                                                                                                                   
 cccGCTCACCATCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCgtgt                                                                                                                                                    <  3:358114/149‑1 (MQ=255)
     cgaaCCATCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGc                                                                                                                                                <  4:100618/146‑1 (MQ=255)
       cACCATCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCt                                                                                                                                              >  4:411139/1‑149 (MQ=255)
        aCCATCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCtt                                                                                                                                             >  1:571635/1‑149 (MQ=255)
         ccATCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTa                                                                                                                                            <  2:571635/149‑1 (MQ=255)
           aTCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGccc                                                                                                                                              <  4:1439683/145‑2 (MQ=255)
           aTCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGccc                                                                                                                                              >  3:1439683/1‑144 (MQ=255)
               aaaCGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGc                                                                                                                                      <  1:511127/149‑1 (MQ=255)
                  cgcgAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCgg                                                                                                                                       >  4:63007/1‑145 (MQ=255)
                       gaTTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTa                                                                                                                                                   <  3:636958/128‑1 (MQ=255)
                        aTTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGcc                                                                                                                                               <  1:762198/131‑1 (MQ=255)
                              cTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTcc                                                                                                                       <  3:623345/149‑1 (MQ=255)
                              cTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTcc                                                                                                                       <  1:1446769/149‑1 (MQ=255)
                               tCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTccc                                                                                                                      >  3:668897/1‑149 (MQ=255)
                               tCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATAAGGCGTACACGccgcatc                                                                                                                                          >  4:760235/1‑124 (MQ=255)
                               tCCTGACAAAAAACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTATCAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCCTATCCGGCCTACGGTCAGGTTccc                                                                                                                      >  3:963800/1‑149 (MQ=255)
                                  tGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCgtgt                                                                                                                                                    >  4:533713/1‑116 (MQ=255)
                                         acacGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGa                                                                                                            <  3:411139/149‑1 (MQ=255)
                                           acGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGTCCTACGGTCAGGTTCCCGTAGGCATGATa                                                                                                          >  4:906102/1‑149 (MQ=255)
                                               cccAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCg                                                                                                                                        <  3:76963/115‑1 (MQ=255)
                                                ccAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGTCCTACGGTCAGGTTCCCGTAGGCATGATAAGAcg                                                                                                     <  3:906102/149‑1 (MQ=255)
                                                ccAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTccc                                                                                                                      <  2:903734/132‑1 (MQ=255)
                                                ccAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGAGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTccc                                                                                                                      >  1:903734/1‑132 (MQ=255)
                                                 cAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGc                                                                                                                                      >  4:766283/1‑115 (MQ=255)
                                                  aGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGAcgcg                                                                                                   <  4:963800/149‑1 (MQ=255)
                                                     aaCGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATcc                                                                                                                                         >  3:1389207/1‑108 (MQ=255)
                                                       cGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGc                                                                                              <  4:663294/149‑1 (MQ=255)
                                                            ttCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTc                                                                                           <  1:218282/147‑1 (MQ=255)
                                                            ttCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTc                                                                                           >  2:218282/1‑147 (MQ=255)
                                                                gTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCa                                                                                     >  2:1114618/1‑149 (MQ=255)
                                                                    gCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCa                                                                                 >  2:1241891/1‑149 (MQ=255)
                                                                          tCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCa                                                                                                               <  1:152414/113‑1 (MQ=255)
                                                                          tCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCa                                                                                                               >  2:152414/1‑113 (MQ=255)
                                                                            tGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAGCAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTc                                                                                           <  3:224999/131‑1 (MQ=255)
                                                                            tGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGcc                                                                                                                                               >  2:1155389/1‑79 (MQ=255)
                                                                            tGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATAc                                                                         >  3:433465/1‑149 (MQ=255)
                                                                              ggTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGc                                                                                              <  2:322512/126‑1 (MQ=255)
                                                                              ggTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGc                                                                                              >  1:322512/1‑126 (MQ=255)
                                                                               gTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGAcgc                                                                                                    <  4:1377618/119‑1 (MQ=255)
                                                                                   ttGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTc                                                                                           >  4:341487/1‑124 (MQ=255)
                                                                                       gcACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTc                                                                                           <  2:787067/120‑1 (MQ=255)
                                                                                           gTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGa                                                                                                       >  1:276073/1‑104 (MQ=255)
                                                                                           gTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCACTGCGAc                                                               <  2:1219237/144‑1 (MQ=255)
                                                                                            tAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGctt                                                         >  1:1216661/1‑149 (MQ=255)
                                                                                             aGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGAt                                                                                                           >  4:1075432/1‑98 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTACCAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAgg                                                                                                                           <  2:397738/80‑1 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAgg                                                                                                                           >  1:397738/1‑80 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTccc                                                                                                                      <  2:1140945/85‑1 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCAt                                                                                                              >  1:95344/1‑93 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCAt                                                                                                              <  2:95344/93‑1 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATaa                                                                                                         <  3:320085/98‑1 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATaa                                                                                                         >  4:320085/1‑98 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGAc                                                                                                      >  1:1145001/1‑101 (MQ=255)
                                                                                               gTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGcttctt                                                      <  2:820600/149‑1 (MQ=255)
                                                                                                tttAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATaa                                                                                                         >  3:1098900/1‑97 (MQ=255)
                                                                                                   aaaaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAg                                                                                               >  4:1182964/1‑104 (MQ=255)
                                                                                                   aaaaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAg                                                                                               <  3:1182964/104‑1 (MQ=255)
                                                                                                    aaaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTc                                                                                           <  4:614940/107‑1 (MQ=255)
                                                                                                    aaaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGAcc                                                              <  2:585279/136‑1 (MQ=255)
                                                                                                     aaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCa                                                >  3:1142265/1‑149 (MQ=255)
                                                                                                      aCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCTTACGCGTc                                                                                           >  1:864139/1‑105 (MQ=255)
                                                                                                        gTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGaa                                             >  2:645515/1‑149 (MQ=255)
                                                                                                         tAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAg                                                                                                                  <  3:1060370/79‑1 (MQ=255)
                                                                                                         tAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCaa                                                                     >  3:187072/1‑124 (MQ=255)
                                                                                                         tAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCaa                                                                     <  4:187072/124‑1 (MQ=255)
                                                                                                         tAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAAc                                            >  1:62776/1‑149 (MQ=255)
                                                                                                            aTTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGAc                                                                                                      >  3:764190/1‑88 (MQ=255)
                                                                                                             ttGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAg                                                                                                                  <  4:439604/75‑1 (MQ=255)
                                                                                                             ttGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGt                                                                                            >  2:1018799/1‑97 (MQ=255)
                                                                                                                tAACAACGTGCGCTTTGTTTATGTCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGTaa                                                                                <  2:1418384/106‑1 (MQ=255)
                                                                                                                    aacGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGccc                                       >  4:256934/1‑143 (MQ=255)
                                                                                                                      cGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGtt                               <  1:484442/149‑1 (MQ=255)
                                                                                                                        tGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTa                             <  4:1142265/149‑1 (MQ=255)
                                                                                                                         gcgcTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTAt                            <  1:645515/149‑1 (MQ=255)
                                                                                                                          cgcTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATc                           <  2:496562/149‑1 (MQ=255)
                                                                                                                            cTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGtt                                                                                                                         >  1:1003351/1‑53 (MQ=255)
                                                                                                                            cTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGAc                                                                                                      >  3:1041986/1‑72 (MQ=255)
                                                                                                                            cTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGAc                                                                                                      <  4:1041986/72‑1 (MQ=255)
                                                                                                                                gTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGc                                                           >  3:607004/1‑111 (MQ=255)
                                                                                                                                 tttATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATg                                                                                                             <  3:341471/60‑1 (MQ=255)
                                                                                                                                 tttATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATAccc                                                                       >  1:803448/1‑98 (MQ=255)
                                                                                                                                  ttATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGc                                                                                              <  1:1024635/74‑1 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTta                                                     >  2:447676/1‑111 (MQ=37)
                                                                                                                                      gCCGGATAAGGCGTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTta                                                     <  1:447676/111‑1 (MQ=37)
                                                                                                                                      gCCGGATAAGGCGTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGcccc                                      >  1:1192195/1‑126 (MQ=38)
                                                                                                                                       ccGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGtt                               <  2:717595/132‑1 (MQ=255)
                                                                                                                                        cGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATc                           <  3:826383/135‑1 (MQ=255)
                                                                                                                                         ggATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAG            >  3:5169‑M1/1‑149 (MQ=255)
                                                                                                                                          gataaggcgtacacGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGC           >  2:878253‑M1/6‑149 (MQ=255)
                                                                                                                                           ataaggcgTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGcttc                                                        <  2:1291559/99‑1 (MQ=37)
                                                                                                                                           ataaggcgTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCa                                                <  4:1465217/107‑1 (MQ=38)
                                                                                                                                           ataaggcgTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGa                        <  2:353363/131‑1 (MQ=39)
                                                                                                                                           ataaggcgTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGa                        >  1:353363/5‑135 (MQ=39)
                                                                                                                                            taaggcgTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGt                                <  1:1240572/123‑1 (MQ=38)
                                                                                                                                            tGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGGCGCGTAAGCGTCGCATCCGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACcagcag         >  1:1373138/1‑149 (MQ=255)
                                                                                                                                             gcggcgTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAAt                                                                    >  2:199496/1‑89 (MQ=255)
                                                                                                                                             gcggcgTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAAt                                                                    <  1:199496/89‑1 (MQ=255)
                                                                                                                                               ggcgTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCa                                                                                     >  2:230912/1‑70 (MQ=255)
                                                                                                                                                gcgTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTc                                                                                           >  1:1119249/1‑63 (MQ=255)
                                                                                                                                                gcgTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAAt                                                                    <  4:1063203/86‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCCGCTCACCATCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTG  >  NC_000913/3754726‑3755021

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