Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,736,100 T→A 11.5% intergenic (+141/+21) sodB → / ← ydhP superoxide dismutase, Fe/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,736,1000TA11.5% 303.2 / 27.0 96intergenic (+141/+21)sodB/ydhPsuperoxide dismutase, Fe/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (41/44);  new base A (5/6);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

TGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGAAATTA  >  NC_000913/1735953‑1736248
                                                                                                                                                   |                                                                                                                                                    
tGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATc                                                                                                                                                     <  4:1095599/149‑1 (MQ=255)
       taaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGCTATCCTTCTGc                                                                                                                                              <  1:767733/149‑1 (MQ=255)
       taaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGc                                                                                                                                              <  3:667148/149‑1 (MQ=255)
           tGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttt                                                                                                                                          >  3:1303643/1‑149 (MQ=255)
           tGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttt                                                                                                                                          >  1:245622/1‑149 (MQ=255)
             aTGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttttt                                                                                                                                        <  1:258865/149‑1 (MQ=255)
                gCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGtt                                                                                                                                     >  1:81900/1‑149 (MQ=255)
                      gcagcaTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTa                                                                                                                               <  3:77523/149‑1 (MQ=255)
                         gcaTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCt                                                                                                                            <  3:1271205/149‑1 (MQ=255)
                           aTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGt                                                                                                                          >  2:512578/1‑149 (MQ=255)
                           aTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGt                                                                                                                          >  2:870832/1‑149 (MQ=255)
                            ttGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGtt                                                                                                                         <  1:512578/149‑1 (MQ=255)
                            ttGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGtt                                                                                                                         <  1:870832/149‑1 (MQ=255)
                               ccTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGc                                                                                                                                              <  1:390364/125‑1 (MQ=255)
                               ccTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGc                                                                                                                                              >  2:390364/1‑125 (MQ=255)
                                 tGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                                                                             >  4:1367218/1‑124 (MQ=255)
                                 tGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                                                                             <  3:1367218/124‑1 (MQ=255)
                                       gcTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaaa                                                                                                              >  1:1440857/1‑149 (MQ=255)
                                        cTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAAc                                                                                                             <  2:245622/149‑1 (MQ=255)
                                         tACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGc                                                                                                                             <  4:608239/132‑1 (MQ=255)
                                         tACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGc                                                                                                                             >  3:1134411/1‑132 (MQ=255)
                                         tACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGc                                                                                                                             >  3:608239/1‑132 (MQ=255)
                                         tACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGc                                                                                                                             <  4:1134411/132‑1 (MQ=255)
                                          aCGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                                                                                                                  >  3:1307902/1‑126 (MQ=255)
                                          aCGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                                                                                                                  <  4:1307902/126‑1 (MQ=255)
                                            gCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCt                                                                                                                            >  1:1186890/1‑130 (MQ=255)
                                            gCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCt                                                                                                                            <  2:1186890/130‑1 (MQ=255)
                                            gCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGtt                                                                                                         >  4:1153614/1‑149 (MQ=255)
                                            gCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGtt                                                                                                         <  2:1440857/149‑1 (MQ=255)
                                             cTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGttt                                                                                                        <  3:1153614/149‑1 (MQ=255)
                                                 tCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAgg                                                                                                    >  3:613894/1‑149 (MQ=255)
                                                   aGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgtt                                                                                                  <  4:1303643/149‑1 (MQ=255)
                                                        tACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGtttt                                                                                             <  2:81900/149‑1 (MQ=255)
                                                              gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAAt                                                                                                                                 >  4:420105/1‑107 (MQ=255)
                                                              gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAAt                                                                                                                                 <  3:420105/107‑1 (MQ=255)
                                                              gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCt                                                                                       >  2:601245/1‑149 (MQ=255)
                                                               atcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTg                                                                                      <  1:601245/149‑1 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttttt                                                                                                                                        >  2:1416146/1‑96 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttttt                                                                                                                                        <  1:1416146/96‑1 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGt                                                                                                                          >  2:115832/1‑110 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGt                                                                                                                          <  1:115832/110‑1 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGttt                                                                                                        <  1:1187440/128‑1 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGttt                                                                                                        >  2:1187440/1‑128 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGt                                                                                                          >  4:1035742/1‑126 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGt                                                                                                          <  3:1035742/126‑1 (MQ=255)
                                                                  atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACa                                                                                   >  3:1303667/1‑149 (MQ=255)
                                                                         tAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAacac                                                                            >  4:826739/1‑149 (MQ=255)
                                                                               ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaaa                                                                                                              <  3:1080595/109‑1 (MQ=255)
                                                                               ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaaa                                                                                                              >  4:1411251/1‑109 (MQ=255)
                                                                               ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaaa                                                                                                              <  3:1411251/109‑1 (MQ=255)
                                                                               ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaaa                                                                                                              >  4:1080595/1‑109 (MQ=255)
                                                                                 aTAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAgg                                                                                                    >  3:464916/1‑117 (MQ=255)
                                                                                 aTAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAgg                                                                                                    <  4:464916/117‑1 (MQ=255)
                                                                                 aTAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCa                                                                          <  3:273651/143‑1 (MQ=255)
                                                                                 aTAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCa                                                                          >  4:273651/1‑143 (MQ=255)
                                                                                     ggCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCt                                                                                       <  1:971168/126‑1 (MQ=255)
                                                                                     ggCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCt                                                                                       >  2:971168/1‑126 (MQ=255)
                                                                                     ggCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGc                                                                <  4:613894/149‑1 (MQ=255)
                                                                                      gCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                                                                                                       <  4:24827/77‑1 (MQ=255)
                                                                                      gCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                                                                                                       >  3:24827/1‑77 (MQ=255)
                                                                                       cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGtt                                                                                                         >  4:8709/1‑106 (MQ=255)
                                                                                       cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGtt                                                                                                         <  3:8709/106‑1 (MQ=255)
                                                                                       cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCa                                                              <  4:1303667/149‑1 (MQ=255)
                                                                                        aCTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTg                                                                                                                           >  4:80305/1‑87 (MQ=255)
                                                                                        aCTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTg                                                                                                                           <  3:80305/87‑1 (MQ=255)
                                                                                        aCTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTATCAACGCAAACTGTTTCAGgttgtt                                                                                               >  1:223911/1‑115 (MQ=255)
                                                                                        aCTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTATCAACGCAAACTGTTTCAGgttgtt                                                                                               <  2:223911/115‑1 (MQ=255)
                                                                                        aCTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGtta                                                                                                                                          >  3:252768/1‑71 (MQ=37)
                                                                                        aCTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGtta                                                                                                                                          <  4:252768/72‑2 (MQ=37)
                                                                                          tCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAaca                                                                             <  1:1111947/131‑1 (MQ=255)
                                                                                          tCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAaca                                                                             >  2:1111947/1‑131 (MQ=255)
                                                                                                 aTCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaaa                                                                               <  1:1080/122‑1 (MQ=255)
                                                                                                 aTCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaaa                                                                               >  2:1080/1‑122 (MQ=255)
                                                                                                    cGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAAt                                                 <  3:826739/149‑1 (MQ=255)
                                                                                                           aTAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                                                                                                                  <  3:294190/61‑1 (MQ=255)
                                                                                                           aTAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                                                                                                                  >  4:294190/1‑61 (MQ=255)
                                                                                                           aTAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGc                                                                                                                                              <  3:167074/49‑1 (MQ=25)
                                                                                                           aTAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGc                                                                                                                                              >  4:167074/1‑49 (MQ=25)
                                                                                                                tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCAc                                     >  2:300305/1‑149 (MQ=255)
                                                                                                                            aGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAg                                                                                                     <  3:1280087/73‑1 (MQ=255)
                                                                                                                            aGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAg                                                                                                     >  4:1280087/1‑73 (MQ=255)
                                                                                                                               aaCTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGc                                                                                        >  1:449146/1‑83 (MQ=255)
                                                                                                                               aaCTATAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGc                                                                                        <  2:449146/83‑1 (MQ=255)
                                                                                                                                  tGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCt                                                                                           <  1:1354314/77‑1 (MQ=255)
                                                                                                                                  tGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCt                                                                                           >  2:1354314/1‑77 (MQ=255)
                                                                                                                                        gCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGc                                                                <  3:570569/98‑1 (MQ=255)
                                                                                                                                        gCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGc                                                                >  4:570569/1‑98 (MQ=255)
                                                                                                                                        gCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATccc                  <  4:1060371/144‑1 (MQ=255)
                                                                                                                                        gCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATccc                  >  3:1060371/1‑144 (MQ=255)
                                                                                                                                        gCAGAAGGATAACCTTCCGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACaa                                                  >  1:521451/1‑112 (MQ=255)
                                                                                                                                        gCAGAAGGATAACCTTCCGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACaa                                                  <  2:521451/112‑1 (MQ=255)
                                                                                                                                         cAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccagccca            >  3:240049/1‑149 (MQ=255)
                                                                                                                                          aGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTc                                                           <  2:1185669/101‑1 (MQ=255)
                                                                                                                                          aGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTc                                                           >  1:1185669/1‑101 (MQ=255)
                                                                                                                                             aGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCgg        >  1:1123479/1‑149 (MQ=255)
                                                                                                                                                  aTCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGAAAtt   <  2:1123479/149‑1 (MQ=255)
                                                                                                                                                   tCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGAAATTa  >  4:1381980/1‑149 (MQ=255)
                                                                                                                                                   tCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGAAATTa  >  1:142450/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCGGAAATTA  >  NC_000913/1735953‑1736248

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