Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 9,908 G→T 15.4% 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,9080GT15.4% 563.5 / 74.0 206intergenic (+15/+20)mog/satPmolybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand):  ref base G (89/84);  new base T (17/15);  total (107/99)
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.97e-01

Only 100 of 209 total aligned reads displayed.

GCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGC  >  NC_000913/9809‑10007
                                                                                                   |                                                                                                   
   ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                  >  1:589374/1‑100 (MQ=255)
   ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                  >  1:3199989/1‑100 (MQ=255)
   ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                  >  1:1602795/1‑100 (MQ=255)
   ggAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                  >  2:2846264/1‑100 (MQ=255)
         gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCGCGGGGGGAt                                                                                            <  1:3081728/100‑1 (MQ=255)
         gCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAt                                                                                            >  2:2277850/1‑100 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                           <  1:3998418/100‑1 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                           >  1:2174720/1‑100 (MQ=255)
          ccATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                           >  2:4221922/1‑100 (MQ=255)
           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                          >  2:3483670/1‑100 (MQ=255)
           cATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                           >  1:3499613/1‑99 (MQ=255)
            aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttt                                                                                         <  2:364717/100‑4 (MQ=255)
            aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtctc                                                                                         <  2:793729/100‑1 (MQ=255)
            aTACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTAAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtct                                                                                          <  1:3918236/99‑1 (MQ=255)
             tACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCCAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCa                                                                                        <  1:2054996/100‑1 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                       <  2:1825215/100‑6 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                       <  2:2495792/100‑6 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                       <  1:1075876/100‑6 (MQ=255)
              aCGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGCGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttcttt                                                                                       <  1:325928/100‑6 (MQ=255)
               cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAttttttt                                                                                      <  2:4157512/100‑7 (MQ=39)
               cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaa                                                                                      <  2:2592951/100‑1 (MQ=255)
               cGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaa                                                                                      >  1:820525/1‑100 (MQ=255)
               cGTTGAAACGGCACCGGAAGTGGTTGCAGCATCCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaa                                                                                      <  1:1336779/100‑1 (MQ=255)
                gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATATAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                     >  1:2590481/1‑100 (MQ=255)
                gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                     <  1:598802/100‑1 (MQ=255)
                gTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaaaa                                                                                     >  2:3699248/1‑100 (MQ=255)
                  tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAATAAAATCCCCCCGCTCGGGGGGAtttttttatt                                                                                   <  2:4184729/100‑10 (MQ=38)
                  tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa                                                                                   <  2:120318/100‑1 (MQ=255)
                  tGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa                                                                                   >  1:2602209/1‑100 (MQ=255)
                  tGAAACGGCACCGGAAGTGGTTGCAGCATACAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttttat                                                                                    <  2:1424821/99‑9 (MQ=38)
                     aaCGGCACCGGAAGTGTTTGCAGCAATCAGACTGACTAGTGCAAGACGCGACTTTAGCGAATAAAAAAATCCCCCCGCTCGGGGGGAtttttt                                                                                       <  1:1858188/93‑6 (MQ=25)
                     aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAtt                                                                                <  1:4240632/100‑1 (MQ=255)
                       cGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAg                                                                              >  1:342895/1‑100 (MQ=255)
                        ggCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGCGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAg                                                                              >  1:1111152/1‑99 (MQ=255)
                         gCACCGGCAGTGATTGCAGAAGTCATACAGATGATTGCAAGAAGAAAAAATAAAAAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                            <  1:1645167/100‑1 (MQ=255)
                         gCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                            <  2:3517276/100‑1 (MQ=255)
                           aCCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTggg                                                                          <  2:2941165/100‑1 (MQ=255)
                            ccGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTTGCGAATAAAAAAATACCCCCTAGCGGGGGGATCTCAAAACAATTAGTGGGa                                                                         >  2:571326/1‑100 (MQ=255)
                             cGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                        >  1:3845216/1‑100 (MQ=255)
                             cGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                        <  2:1125991/100‑1 (MQ=255)
                                 aGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAc                                                                    >  1:2299075/1‑100 (MQ=255)
                                  gTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCAcc                                                                   <  2:2891942/100‑1 (MQ=255)
                                     gTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                                >  2:1154311/1‑100 (MQ=255)
                                       tGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                              <  2:2226552/100‑1 (MQ=255)
                                       tGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCg                                                              <  2:1528366/100‑1 (MQ=255)
                                          agcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCa                                                           >  1:4064048/1‑100 (MQ=255)
                                           gccttCACACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCCCGCGGGGGGATCCCACAACCATTAGTGGCACTCACCAATCGGCAg                                                          >  1:2649055/4‑100 (MQ=255)
                                           gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAg                                                          >  2:3829439/1‑100 (MQ=255)
                                           gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAg                                                          <  2:4054161/100‑1 (MQ=255)
                                            caTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGa                                                         >  2:1059334/1‑100 (MQ=255)
                                              gacAGACCGAAGATGGCAAGACACAACATTAGCGAATAAAAAAATCCCTCCGAACGGGGGTATCTAACAAAAATTAACAGAATTCAACAATCGGCAGaa                                                        >  1:1561013/3‑99 (MQ=255)
                                                  gACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTg                                                   <  2:1645057/100‑1 (MQ=255)
                                                   aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGc                                                  >  2:4240545/1‑100 (MQ=255)
                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCg                                                 >  1:224242/1‑100 (MQ=255)
                                                    ccGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGACAGAACGGTGCg                                                 >  1:1340123/1‑100 (MQ=255)
                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                               >  2:3093110/1‑100 (MQ=255)
                                                      gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAc                                               <  2:187863/100‑1 (MQ=255)
                                                        agagTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCa                                             <  2:1602802/100‑1 (MQ=255)
                                                         gagTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAAAGGTGc                                                  >  1:346913/1‑94 (MQ=255)
                                                           gTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAAc                                          <  2:2174727/100‑1 (MQ=255)
                                                             gCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTg                                        >  2:3200942/1‑100 (MQ=255)
                                                             gCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATATCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTg                                        >  1:1353078/1‑100 (MQ=255)
                                                                 gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCACAACAATTACTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                    <  2:874453/100‑1 (MQ=255)
                                                                 gACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCg                                    <  2:1458259/100‑1 (MQ=255)
                                                                   cgcgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                                  >  2:3916388/1‑100 (MQ=255)
                                                                    gcgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                 <  1:2137531/100‑1 (MQ=255)
                                                                     cgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                <  2:1580386/100‑1 (MQ=255)
                                                                       aCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                              >  2:3739805/1‑100 (MQ=255)
                                                                       aCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                              >  2:3237360/1‑100 (MQ=255)
                                                                         gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCACTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                            <  1:1412160/100‑1 (MQ=255)
                                                                         gTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                              <  1:3647564/98‑1 (MQ=255)
                                                                            aGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                         <  1:3395624/100‑1 (MQ=255)
                                                                                tgttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                     >  2:1352822/10‑100 (MQ=38)
                                                                                aaaaaaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                      <  1:860922/96‑1 (MQ=255)
                                                                                aaTAAAACAAGCCCACCGAGCGTTGGGATCTCAAGACAATTATGTGGCTTCACCAATCGGCAGAGCGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                       <  1:3684567/98‑1 (MQ=255)
                                                                                aaTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                     <  1:1138138/100‑1 (MQ=255)
                                                                                 aTAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                    >  1:4400943/1‑100 (MQ=255)
                                                                                  ttttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                   >  2:2638959/8‑100 (MQ=38)
                                                                                  tAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                   >  2:1734028/1‑100 (MQ=255)
                                                                                   aaaaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                  <  2:3329656/100‑1 (MQ=255)
                                                                                    ttgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                 >  1:4201436/6‑100 (MQ=255)
                                                                                    aaaaaaTCCCCCCGAGCGGGGGGATATCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                  >  2:1547178/1‑99 (MQ=255)
                                                                                     tgagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc                >  2:2519403/5‑100 (MQ=255)
                                                                                     ngagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc                >  2:3271477/5‑100 (MQ=255)
                                                                                     aaaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGCACTTCACCCATCGc                >  1:796612/1‑100 (MQ=255)
                                                                                      gagaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc               >  1:3158601/4‑100 (MQ=255)
                                                                                      aaaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc               >  1:354485/1‑100 (MQ=255)
                                                                                      aaaaTCCCCCCGAGCGGGGGGATCTCAAAAAAAATAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTAACCCATCGcc               >  1:2899932/1‑100 (MQ=255)
                                                                                       aaaTCCCCCCGAGCGGGGGGATCTCCAAACACTTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa              <  2:606764/100‑1 (MQ=255)
                                                                                       aaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa              <  2:389/100‑1 (MQ=255)
                                                                                       aaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa              <  2:1026968/100‑1 (MQ=255)
                                                                                       aaaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa              <  2:1258878/100‑1 (MQ=255)
                                                                                        gaTCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGCTTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag             >  1:3506994/2‑100 (MQ=255)
                                                                                          tctccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagat           <  2:2409621/97‑1 (MQ=255)
                                                                                          tctccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagat           <  2:201405/97‑1 (MQ=255)
                                                                                          tCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACCGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCag             <  1:1012952/98‑1 (MQ=255)
                                                                                             ncccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagataga        >  1:201405/2‑100 (MQ=255)
                                                                                             ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagataga        >  1:2409610/1‑100 (MQ=255)
                                                                                               ccGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCACACTGCTCGTTCAGTACTTCACCCACCGCCAGATAGAtt      <  1:2721540/100‑1 (MQ=255)
                                                                                                 gAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgc    <  1:3737866/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: