Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,197,365 G→A 23.1% pseudogene (892/2847 nt) yehH → DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,197,3650GA23.1% 166.9 / 42.0 86pseudogene (892/2847 nt)yehHDUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base A (10/10);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AGCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTGG  >  NC_000913/2197266‑2197464
                                                                                                   |                                                                                                   
aGCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTg                                                                                                     <  2:338603/100‑1 (MQ=255)
aGCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTg                                                                                                     <  2:1760309/100‑1 (MQ=255)
 gCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGc                                                                                                    >  1:1714757/1‑100 (MQ=255)
          aGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGTTCAAAAACTTATTACCGCAATTcc                                                                                           <  2:1063163/100‑1 (MQ=255)
            cAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCca                                                                                          <  1:1244721/99‑1 (MQ=255)
             aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCcaca                                                                                        <  1:526878/100‑1 (MQ=255)
             aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                        <  1:1480985/100‑1 (MQ=255)
                atTCCAGTTTTGAATGACGCATACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACatat                                                                                     >  2:497977/1‑100 (MQ=255)
                 ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATAtt                                                                                    >  2:1486380/1‑100 (MQ=255)
                 ntCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGACTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATAtt                                                                                    >  2:750555/2‑100 (MQ=255)
                   ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCAACATATTcc                                                                                  >  1:524976/1‑100 (MQ=255)
                   ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGAAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTcc                                                                                  >  2:1528997/1‑100 (MQ=255)
                    cAGATTTTAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCa                                                                                 >  2:1073095/1‑100 (MQ=255)
                    cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCa                                                                                 >  1:1230831/1‑100 (MQ=255)
                    cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                                 >  1:926055/1‑100 (MQ=255)
                      gATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGa                                                                               >  2:773398/1‑100 (MQ=255)
                      gATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGa                                                                               >  2:1187205/1‑100 (MQ=255)
                      gATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGa                                                                               >  1:764039/1‑100 (MQ=255)
                           gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggc                                                                          >  2:164756/1‑100 (MQ=255)
                           gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATgac                                                                          >  1:988567/1‑98 (MQ=255)
                           gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggcg                                                                         <  1:1357176/100‑1 (MQ=255)
                           ctaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggc                                                                          <  1:2020705/98‑1 (MQ=255)
                             aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggcag                                                                        >  2:2020683/1‑98 (MQ=255)
                              tGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcggc                                                                       <  2:736900/100‑1 (MQ=255)
                                aCGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCaa                                                                     <  1:91348/100‑1 (MQ=255)
                                 cGCCTACGTAAACATTTATCACTGACGGAAGAGTATCTCTGACAAAAATGTGCACAAAAACGTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                    <  2:1472569/100‑1 (MQ=255)
                                  gCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAAcc                                                                   <  1:191977/100‑1 (MQ=255)
                                  gCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAAcc                                                                   <  1:1529015/100‑1 (MQ=255)
                                   ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGACAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                                  >  2:397317/1‑100 (MQ=255)
                                   ccTACGTAAACATTTATCACTGACGGAAGAGAAAAAAAAAAAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                    <  2:1068674/98‑1 (MQ=255)
                                     tACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCAtt                                                                >  2:347445/1‑100 (MQ=255)
                                      aCGTAAACATTTATCACTGACGGAAGTGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTa                                                               <  2:1156885/100‑1 (MQ=255)
                                       cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTa                                                               >  2:328580/1‑99 (MQ=255)
                                       cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTaa                                                              >  1:519064/1‑100 (MQ=255)
                                        gTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAt                                                             >  1:2204116/1‑100 (MQ=255)
                                         tAAACATTTATCAATGACGTAAGAGTCTCTCTGGCAAAAATGTGCACAAAAGCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTa                                                               <  2:1944854/97‑1 (MQ=255)
                                          aaaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                           >  1:216691/1‑100 (MQ=255)
                                          aaaCATTTATCACTGACGGAAGAGTCGCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                           <  1:2173132/100‑1 (MQ=255)
                                            aCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCt                                                         <  2:2144465/100‑1 (MQ=255)
                                             caattATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAtt                                                            >  2:2302282/4‑96 (MQ=255)
                                               tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                      <  1:1187216/100‑1 (MQ=255)
                                               tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                      <  1:1878827/100‑1 (MQ=255)
                                               tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                      <  1:2058538/100‑1 (MQ=255)
                                                ttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGt                                                     <  2:1663612/100‑1 (MQ=255)
                                                ttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGt                                                     >  2:1193585/1‑100 (MQ=255)
                                                 tATCACTGACTGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATTTTCCAGAATGGCGGCAACCATTAATGGCTTTGtn                                                    <  2:1482768/100‑2 (MQ=255)
                                                  aTCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                                   <  2:1575033/100‑1 (MQ=255)
                                                   tCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttat                                                  <  1:773407/100‑1 (MQ=255)
                                                    cACTGACGTAAGAGTCGCTCTGGCAAAAATGGGCACAAAAACTTATTGCCGCAATTCCACAGATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatt                                                 >  2:1428177/1‑100 (MQ=255)
                                                    cACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGATTGGCGGCAACCATTAATTGCTTTGttatt                                                 >  1:731275/1‑100 (MQ=255)
                                                       tGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAcc                                              >  1:141150/1‑100 (MQ=255)
                                                        gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                             >  1:1800530/1‑100 (MQ=255)
                                                         aCGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCg                                            <  2:1167530/100‑1 (MQ=255)
                                                          cGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGTTTTGTTATTACCCGa                                           >  1:1725431/1‑100 (MQ=255)
                                                            gAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaa                                         <  1:477311/100‑1 (MQ=255)
                                                            gAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaa                                         >  1:1413349/1‑100 (MQ=255)
                                                              agagTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaaaa                                       >  2:1385478/1‑100 (MQ=255)
                                                                agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                     <  1:982063/100‑1 (MQ=255)
                                                                   ctctctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                                  >  1:1730564/1‑100 (MQ=255)
                                                                    tctctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaa                                 <  1:1057457/100‑1 (MQ=255)
                                                                     ctctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCGCATATTCCAGAATGGCGGCGACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaaa                                >  2:1303003/1‑100 (MQ=255)
                                                                     ctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaaa                                <  2:69697/100‑1 (MQ=255)
                                                                        tGGTAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTg                             <  1:1320290/100‑1 (MQ=255)
                                                                        tGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTg                             >  2:2209182/1‑100 (MQ=255)
                                                                        aggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTg                             >  2:1849982/2‑100 (MQ=255)
                                                                          gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                           <  2:369223/100‑1 (MQ=255)
                                                                           cAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                          <  1:1377319/100‑1 (MQ=255)
                                                                           cAAAAATGTGCACAAAAACTTATTGCCACAATTACACATATTCCAGAATGGCGGCAACCATTAATTGCTTTATTATTACCCGAAAAACCAGAAATTGca                           >  2:1005753/1‑99 (MQ=255)
                                                                           cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                          <  2:141147/100‑1 (MQ=255)
                                                                           cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                          >  1:4937/1‑100 (MQ=255)
                                                                            aaaaaTGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcaca                         <  2:886106/100‑1 (MQ=255)
                                                                            aaaaaTGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                          <  2:2236974/99‑1 (MQ=255)
                                                                                   tgCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATT                  >  1:1169026‑M1/1‑100 (MQ=255)
                                                                                    gCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTG                 >  2:657390‑M1/1‑100 (MQ=255)
                                                                                    gCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATT                  >  1:1373676‑M1/1‑99 (MQ=255)
                                                                                     cacaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGC                >  1:2098505‑M1/1‑100 (MQ=255)
                                                                                     cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGC                >  2:1799065‑M1/1‑100 (MQ=255)
                                                                                     cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGC                >  1:2253198‑M1/1‑100 (MQ=255)
                                                                                      acaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCC               <  2:216688‑M1/100‑1 (MQ=255)
                                                                                      acaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCC               <  2:2108973‑M1/100‑1 (MQ=255)
                                                                                      acaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCC               <  1:1035233‑M1/100‑1 (MQ=255)
                                                                                      acaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCC               >  2:1197461‑M1/1‑100 (MQ=255)
                                                                                       caAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCC              >  1:1484418‑M1/1‑100 (MQ=255)
                                                                                       caAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCC              >  2:204783‑M1/1‑100 (MQ=255)
                                                                                       caAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCC              <  1:696533‑M1/100‑1 (MQ=255)
                                                                                       caAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACACGAAAAACCAGAAATTGCACATGAAATTGCCC              >  2:902477‑M1/1‑100 (MQ=255)
                                                                                        aaaaaCTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCa             <  1:2402405/100‑1 (MQ=255)
                                                                                           aaCTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCg          <  2:2204090/100‑1 (MQ=255)
                                                                                                   gCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTTCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTGg  <  2:436395‑M2/100‑2 (MQ=255)
                                                                                                   gCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTGg  <  2:2083513‑M2/100‑2 (MQ=255)
                                                                                                   |                                                                                                   
AGCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTGG  >  NC_000913/2197266‑2197464

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: