Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,786,725 C→T 24.9% pseudogene (329/2253 nt) ygaQ → uncharacterized protein; putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,786,7250CT24.9% 248.9 / 78.6 141pseudogene (329/2253 nt)ygaQuncharacterized protein; putative enzyme
Reads supporting (aligned to +/- strand):  ref base C (44/59);  new base T (13/22);  total (60/81)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 144 total aligned reads displayed.

ATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAA  >  NC_000913/2786626‑2786824
                                                                                                   |                                                                                                   
aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAAc                                                                                                     <  1:919276/100‑1 (MQ=255)
  tataAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAatat                                                                                                   >  1:139217/1‑100 (MQ=255)
     aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAgg                                                                                                <  1:560905/100‑1 (MQ=255)
     aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAgg                                                                                                >  2:1712413/1‑100 (MQ=255)
     aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAgg                                                                                                <  2:2363167/100‑1 (MQ=255)
         tttCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAATCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTa                                                                                            <  2:659756/100‑1 (MQ=255)
           tCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATc                                                                                          <  1:329552/100‑1 (MQ=255)
           tCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATc                                                                                          >  1:1498603/1‑100 (MQ=255)
            cccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCa                                                                                         >  1:1993836/1‑100 (MQ=255)
            cccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCa                                                                                         >  1:67710/1‑100 (MQ=255)
             ccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCaa                                                                                        <  2:2020165/100‑1 (MQ=255)
             ccTTTCTCCAGATACAACTTACCGAGAGTTGAAGAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATTGGTTTATCaa                                                                                        >  2:1570281/1‑100 (MQ=255)
              cTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAAc                                                                                       >  2:345901/1‑100 (MQ=255)
               tttCTACAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGAAAAAAGAAAAAAAAATATAGGTTTATCAAc                                                                                       <  2:2224096/99‑1 (MQ=255)
                ttCTCCAGATACAACTTACCGAGAGTTGAATAAATCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGCTTTATCAACga                                                                                     <  2:1255395/100‑1 (MQ=255)
                ttCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAAAATAGGTTTATCAACga                                                                                     >  2:1218730/1‑100 (MQ=255)
                ttCTCCAGATACAACTAACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACg                                                                                      <  1:1675038/99‑1 (MQ=255)
                ttCTCCAGACAAACCTTACCGAGAGATGCATAAAGCACCACAGCGTCGCCTTAAACCGCCCGACAGGCTCCGAGAAAGACAAACATAGGTTTATCAACg                                                                                      <  2:2004360/99‑1 (MQ=255)
                   tcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGAt                                                                                  <  1:1041066/100‑1 (MQ=255)
                   tcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGAt                                                                                  <  1:567843/100‑1 (MQ=255)
                     cAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATAc                                                                                >  1:557194/1‑100 (MQ=255)
                         tACAACTTACCGAGAGTTGAATAAAGCAGCACAGTGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATTGGTTTATCAACGAGATACAAt                                                                             <  2:1005851/99‑1 (MQ=255)
                          aCAACTTGCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATg                                                                            <  1:915126/99‑1 (MQ=255)
                           cAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGAt                                                                          <  1:2185802/100‑1 (MQ=255)
                             aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGAt                                                                          <  2:1508506/98‑1 (MQ=255)
                             aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTg                                                                        >  1:1567476/1‑100 (MQ=255)
                              cTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGa                                                                       >  1:118949/1‑100 (MQ=255)
                               ttACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGaa                                                                      <  1:935283/100‑1 (MQ=255)
                                tACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATAGaaa                                                                     >  2:1062438/1‑100 (MQ=255)
                                 aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAag                                                                    >  1:178679/1‑100 (MQ=255)
                                 aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAag                                                                    <  2:1660921/100‑1 (MQ=255)
                                    gagaGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAaga                                                                   >  2:2084041/1‑98 (MQ=255)
                                    gagaGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAgg                                                                 >  2:1014347/1‑100 (MQ=255)
                                      gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTg                                                               >  1:2322622/1‑100 (MQ=255)
                                       agTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTg                                                               <  1:1695468/99‑1 (MQ=255)
                                        gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACAAGATACAATGATTAAAAGAGGtaa                                                              >  1:1926783/1‑97 (MQ=255)
                                        gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGAt                                                             >  1:2119558/1‑100 (MQ=255)
                                         ttGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATa                                                            <  2:270390/100‑1 (MQ=255)
                                         ttGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATa                                                            <  2:460249/100‑1 (MQ=255)
                                          tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATaa                                                           >  2:1220269/1‑100 (MQ=255)
                                           gAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAAc                                                          >  2:1842617/1‑100 (MQ=255)
                                           gAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAAc                                                          >  1:620193/1‑100 (MQ=255)
                                              tAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGc                                                       >  2:84458/1‑100 (MQ=255)
                                              tAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGc                                                       <  2:340722/100‑1 (MQ=255)
                                               aaaaCAACACAACGTAGCATTACCGCGCCAAACAGGATAGGAGAAAGAAAAACATAATTCTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGc                                                       <  2:2381430/99‑1 (MQ=255)
                                                  gcagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGgtg                                                   >  2:50030/1‑100 (MQ=255)
                                                   cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGgtg                                                   <  1:2119515/99‑1 (MQ=255)
                                                   cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGgtgt                                                  >  1:1072808/1‑100 (MQ=255)
                                                    agcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGtt                                                 >  2:1258683/1‑100 (MQ=255)
                                                    agcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGtt                                                 <  1:2288401/100‑1 (MQ=255)
                                                    agcaCAACGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGtt                                                 <  1:1883084/100‑1 (MQ=255)
                                                     gcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGttt                                                <  1:1824477/100‑1 (MQ=255)
                                                     gcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGttt                                                <  1:313516/100‑1 (MQ=255)
                                                        caGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATc                                             <  2:816070/100‑1 (MQ=255)
                                                        caGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACAAAGGTTTATCAAAGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTa                                               >  2:2157883/1‑98 (MQ=255)
                                                         aGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                              <  1:26352/98‑1 (MQ=255)
                                                          gCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAg                                           <  2:2045577/100‑1 (MQ=255)
                                                          gCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAg                                           <  1:1772868/100‑1 (MQ=255)
                                                          gCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAg                                           <  1:1221696/100‑1 (MQ=255)
                                                          gCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAg                                           <  1:242963/100‑1 (MQ=255)
                                                          gCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAg                                           <  1:1412912/100‑1 (MQ=255)
                                                             tAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTAt                                        <  1:1346556/100‑1 (MQ=255)
                                                             tAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTAt                                        >  1:519526/1‑100 (MQ=255)
                                                              aGCATTACAGAACCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTAAACGAGGTGATAACAGCGGTGTTTATCAGTATg                                       <  2:1579181/100‑1 (MQ=255)
                                                               gCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                      >  2:913088/1‑100 (MQ=255)
                                                                cATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATAGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                      >  2:1174190/1‑99 (MQ=255)
                                                                  ttACAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCg                                   <  1:775989/100‑1 (MQ=255)
                                                                    aCAGCGCNAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGt                                  >  2:23286/1‑99 (MQ=255)
                                                                    aCAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                 >  2:281243/1‑100 (MQ=255)
                                                                     cAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                 >  1:696591/1‑99 (MQ=255)
                                                                     cAGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTAGCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGc                                <  1:463960/100‑1 (MQ=255)
                                                                      aGCGCCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCt                               <  1:288428/100‑1 (MQ=255)
                                                                       gcgcCAGACAGGATAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTg                              <  2:693334/100‑1 (MQ=255)
                                                                        cgcCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGa                             <  2:2210857/100‑1 (MQ=255)
                                                                          ccAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                            >  2:188007/1‑99 (MQ=255)
                                                                          ccAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAc                           >  1:2258920/1‑100 (MQ=255)
                                                                             taCAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAc                           <  2:1507596/96‑1 (MQ=255)
                                                                                 ggATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                    <  1:1712103/100‑1 (MQ=255)
                                                                                 ggATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                    >  1:1574212/1‑100 (MQ=255)
                                                                                 ggATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAAATTTCa                     >  1:1392582/1‑99 (MQ=255)
                                                                                 ggATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTAAAAGAAGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                    >  2:1191052/1‑100 (MQ=255)
                                                                                  gATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                   <  1:397787/100‑1 (MQ=255)
                                                                                   aTAGGAGAAAGAAAAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAccc                  <  2:951595/100‑1 (MQ=255)
                                                                                    tAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCca                 <  1:620929/100‑1 (MQ=255)
                                                                                    tAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCca                 >  1:2057733/1‑100 (MQ=255)
                                                                                      ggAGAAATAAAATCATAGGTTTATCAACGAGATAGAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCTCCtact               >  2:1338242/1‑99 (MQ=255)
                                                                                        agaaagaaAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACAtt             >  1:361087/1‑100 (MQ=255)
                                                                                        agaaagaaAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACAt              <  1:866741/99‑1 (MQ=255)
                                                                                           aagaaaAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata          >  1:1521870/1‑100 (MQ=255)
                                                                                           aagaaaAATATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata          <  1:1658216/100‑1 (MQ=255)
                                                                                           aagaaaAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata          >  2:940867/1‑100 (MQ=255)
                                                                                           aagaaaAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata          <  2:1584626/100‑1 (MQ=255)
                                                                                           aagaaaAACATAGGTTTATCAACAAGATACAATGAATAAAAAAGGTGATAACACCGGTTATTAACAGTCTGATCGTGCTGAACATTTCACCCACATtata          >  2:1184478/1‑100 (MQ=255)
                                                                                            agaaaAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTatat         >  2:53318/1‑100 (MQ=255)
                                                                                            agaaaAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTatat         >  2:1040304/1‑100 (MQ=255)
                                                                                             gaaaAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATc        <  1:1315684/100‑1 (MQ=255)
                                                                                             gaaaAACATAGGTTTATAAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACTCACATTATATc        >  2:906124/1‑100 (MQ=255)
                                                                                              aaaaaCATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCt       <  1:1045761/100‑1 (MQ=255)
                                                                                               aaaaCATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTg      <  1:1650905/100‑1 (MQ=255)
                                                                                                   cATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa  <  2:1287502/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAA  >  NC_000913/2786626‑2786824

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: