Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,908 G→T 16.8% intergenic (+15/+20) mog → / ← satP molybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,9080GT16.8% 327.5 / 54.1 126intergenic (+15/+20)mog/satPmolybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand):  ref base G (59/45);  new base T (12/9);  total (72/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

Only 100 of 128 total aligned reads displayed.

GCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGC  >  NC_000913/9809‑10007
                                                                                                   |                                                                                                   
gCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCACTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAg                                                                                                     >  1:555250/1‑100 (MQ=255)
    gAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggg                                                                                                 <  1:1025633/100‑1 (MQ=255)
    gAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggg                                                                                                 >  1:2300050/1‑100 (MQ=255)
     aaGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                                >  1:111371/1‑100 (MQ=255)
     aaGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                                <  1:1313600/100‑1 (MQ=255)
         gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCAAAAAAAGAAAAAAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                            <  1:349463/100‑1 (MQ=255)
         gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGAACGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                            <  2:2272022/100‑1 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGGGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                           <  1:1548140/100‑1 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                           <  2:1515915/100‑1 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                           <  1:167735/100‑1 (MQ=255)
           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttt                                                                                          <  2:818810/100‑3 (MQ=255)
           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttt                                                                                          <  1:2026047/100‑3 (MQ=255)
           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttt                                                                                          <  2:2281186/100‑3 (MQ=255)
            aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttt                                                                                         <  1:2228311/100‑4 (MQ=255)
             tACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                        >  1:643997/1‑100 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                       <  1:489655/100‑6 (MQ=255)
                gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                     <  1:178631/100‑1 (MQ=255)
                 ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttttat                                                                                    <  1:1832172/100‑9 (MQ=38)
                 ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                    >  1:1825046/1‑100 (MQ=255)
                 ttGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAAc                                                                                    <  1:1934448/100‑1 (MQ=255)
                  tGAAACGGCACCGCGAGTGGTTGCGGCATTCAGAAAGACGAGGGCAAGAGGCGACGTGAGCGCATAAAGAAATTCCCCCGAGCGGGGGGAgctcgaca                                                                                     >  2:1419492/1‑94 (MQ=255)
                    aaaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAt                                                                                 <  1:1843581/100‑1 (MQ=255)
                     aaCGGCACCGGACGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATCAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAtt                                                                                >  1:2397744/1‑100 (MQ=255)
                       cGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAg                                                                              >  1:2358489/1‑100 (MQ=255)
                       cGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAg                                                                              >  2:1281476/1‑100 (MQ=255)
                        ggCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGt                                                                             >  2:601223/1‑100 (MQ=255)
                         gCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                            >  2:2089279/1‑100 (MQ=255)
                          cACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTgg                                                                           <  1:1078289/100‑1 (MQ=255)
                          cACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTgg                                                                           >  2:1221669/1‑100 (MQ=255)
                              ggAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTCGCGCATAAAAAAATCCCCCCGAACGGGGGGATCTCAAAACAATTAGTGGGAtt                                                                       >  2:1564148/1‑100 (MQ=255)
                                  gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                   <  1:2090006/100‑1 (MQ=255)
                                  gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                   <  1:2452832/100‑1 (MQ=255)
                                   tGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCa                                                                  <  1:61417/100‑1 (MQ=255)
                                   tGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCa                                                                  >  2:2475303/1‑100 (MQ=255)
                                    ggTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCaa                                                                 <  1:192604/100‑1 (MQ=255)
                                      ttGCAGCATTCAGACCGAAGCGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATc                                                               <  2:1110783/100‑1 (MQ=255)
                                      ttGCAGCATTCAGACCGAAGAGTGCAAGCCGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATc                                                               <  1:1817694/100‑1 (MQ=255)
                                        gcagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCgg                                                             >  2:2238079/1‑100 (MQ=255)
                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                            >  1:1241865/1‑100 (MQ=255)
                                         cagcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                            >  1:1240237/1‑100 (MQ=255)
                                          agcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                           <  1:263999/100‑1 (MQ=255)
                                           gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAg                                                          >  2:2325727/1‑100 (MQ=255)
                                           gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAg                                                          <  2:836008/100‑1 (MQ=255)
                                           gcaTTCAGAACGAAGAGTGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGCTCTCAAAACAATTAGTGGGATTCACCACGCGGc                                                            >  2:1249063/1‑98 (MQ=255)
                                             aTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCCCAAAACAATTAGTGGGATTCACCAATCGGCAGaa                                                        >  2:1882894/1‑100 (MQ=255)
                                              ttCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAc                                                       >  2:1164057/1‑100 (MQ=255)
                                                cAGACCGACGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCGCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                                     <  2:2383537/100‑1 (MQ=255)
                                                cAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                                     >  2:923260/1‑100 (MQ=255)
                                                 aGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGt                                                    >  1:1574626/1‑100 (MQ=255)
                                                 aGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGt                                                    >  1:1574596/1‑100 (MQ=255)
                                                  gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                   >  2:294249/1‑100 (MQ=255)
                                                  gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                   >  1:445418/1‑100 (MQ=255)
                                                  gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAGCGGCAGAACGGTg                                                   >  2:2455717/1‑100 (MQ=255)
                                                   aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGc                                                  <  1:1449615/100‑1 (MQ=255)
                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGCACGGTGCg                                                 >  2:1655723/1‑100 (MQ=255)
                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                 <  1:1221677/100‑1 (MQ=255)
                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                               <  1:2302388/100‑1 (MQ=255)
                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                               <  1:1816171/100‑1 (MQ=255)
                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                               <  2:2190093/100‑1 (MQ=255)
                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                               >  2:1863672/1‑100 (MQ=255)
                                                          agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACGaaa                                           >  2:2233916/1‑100 (MQ=255)
                                                          agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCaaa                                           <  1:981325/100‑1 (MQ=255)
                                                           gTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                          <  1:1216880/100‑1 (MQ=255)
                                                             gCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTg                                        >  2:342100/1‑100 (MQ=255)
                                                              cAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGc                                       >  2:1600846/1‑100 (MQ=255)
                                                                 gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                    <  2:2457618/100‑1 (MQ=255)
                                                                 gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                    <  2:365012/100‑1 (MQ=255)
                                                                   cgcgACGTTAGCGAATAAAAAAATCACCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAACGGCAGAACGGTGCGACCAAACTGCTCg                                    >  2:665360/1‑98 (MQ=255)
                                                                    gcgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                 <  1:192465/100‑1 (MQ=255)
                                                                     cgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                >  1:1451591/1‑100 (MQ=255)
                                                                      gACGTTAGCGAAAAAAAAAATCCCCCCGAGCTGGGGGATCTCAAAACGAGTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                    <  2:577831/95‑1 (MQ=255)
                                                                       aCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                              <  1:1490505/100‑1 (MQ=255)
                                                                        cGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                             <  2:1278103/100‑1 (MQ=255)
                                                                            aGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                         >  1:1913707/1‑100 (MQ=255)
                                                                             gCGAATAAAAAAAGCCCCCCGAGCGGGGGGATATCAAAACAATTAGTAGGATTCACCAATCGGCAGAACGGTGCGAGCAAACTGCTCGTTAAGTACt                           >  2:1961823/1‑97 (MQ=255)
                                                                              cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCGCCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                       >  1:1041999/1‑100 (MQ=255)
                                                                               ttgttttgagaTCCCCCCGCTCGGGGGGATTTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                      >  1:97538/11‑100 (MQ=37)
                                                                                 gttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGCACTCACCAATCGGCAGAACGGAGCGACAAAACTGCTCGTTCAGTACTCCa                        >  2:588966/9‑96 (MQ=25)
                                                                                  tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                   >  1:2097774/1‑100 (MQ=255)
                                                                                  tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                   >  1:89062/1‑100 (MQ=255)
                                                                                  tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                   >  1:2417748/1‑100 (MQ=255)
                                                                                   tttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                  >  1:2113335/7‑100 (MQ=39)
                                                                                   aaaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                  >  2:1365162/1‑100 (MQ=255)
                                                                                   aaaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                  >  2:1334295/1‑100 (MQ=255)
                                                                                    ttgagaTCCCCCCGCTCGGGGGGATTTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                 >  1:474888/6‑100 (MQ=38)
                                                                                    ttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                 >  1:1922106/6‑100 (MQ=255)
                                                                                     tctgaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc                <  2:2514369/96‑1 (MQ=255)
                                                                                      gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc               >  1:362388/4‑100 (MQ=255)
                                                                                      gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc               >  2:2543731/4‑100 (MQ=255)
                                                                                      aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc               <  1:486398/100‑1 (MQ=255)
                                                                                        gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag             >  1:781736/2‑100 (MQ=255)
                                                                                        gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag             >  2:1066861/2‑100 (MQ=255)
                                                                                        gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag             >  2:1304396/2‑100 (MQ=255)
                                                                                           ccccccGAGCGTGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAGCGCCagata          <  2:1720922/100‑1 (MQ=255)
                                                                                           ccccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata          >  1:797160/1‑100 (MQ=255)
                                                                                           ccccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagata          <  2:2097760/100‑1 (MQ=255)
                                                                                            cccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatag         >  1:1263495/1‑100 (MQ=255)
                                                                                             ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagataga        >  1:1381084/1‑100 (MQ=255)
                                                                                                cGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCGGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTg     >  2:2537116/1‑100 (MQ=255)
                                                                                                 gAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgc    <  1:88969/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGC  >  NC_000913/9809‑10007

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