Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,728,444 C→A 100% T33N (ACC→AAC)  rnt → RNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,728,4440CA99.3% 528.8 / ‑4.1 156T33N (ACC→AAC) rntRNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (83/71);  total (84/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

GATGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCAT  >  NC_000913/1728346‑1728539
                                                                                                  |                                                                                               
gATGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAcc                                                                                                <  2:2535431/100‑1 (MQ=255)
   gTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGAt                                                                                             <  2:1017662/100‑1 (MQ=255)
   gTCCGATAACACTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGAt                                                                                             <  1:2805057/100‑1 (MQ=255)
    tCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATg                                                                                            >  1:1150242/1‑100 (MQ=255)
     ccGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCACAAACGATgc                                                                                           >  1:3779537/1‑100 (MQ=255)
       gATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATgcgc                                                                                         >  2:2565724/1‑100 (MQ=255)
           aCGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCgctgct                                                                                     <  1:3078896/100‑1 (MQ=255)
           aCGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCgctgct                                                                                     >  2:4108066/1‑100 (MQ=255)
            cGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCtt                                                                                    >  2:2213403/1‑100 (MQ=255)
             gCTCACCTTACCGGTCTGTACGACCGTTTACGTGGCTTCTATCCTGTTGTGATCGACGTTGAAACAGCCGGACTTAACGCCAAAAACGATGCTCAGCtt                                                                                    >  1:120595/1‑99 (MQ=255)
             gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTg                                                                                   >  2:583020/1‑100 (MQ=255)
              cTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTga                                                                                  >  1:1547479/1‑100 (MQ=255)
              cTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTATATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTga                                                                                  >  1:2273332/1‑100 (MQ=255)
              cTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGATTTTATCCTGTTGTGATCGATGTTGACACAGCCGGCTATAACGCCAAAAACCATGCGCTGCtt                                                                                    >  1:2106642/1‑98 (MQ=255)
                cAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTgaga                                                                                <  2:2497592/100‑1 (MQ=255)
                cAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTgaga                                                                                <  2:1458069/100‑1 (MQ=255)
                  aCTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGAtt                                                                              >  1:2314749/1‑100 (MQ=255)
                   gaaaacGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTg                                                                             >  1:2289964/6‑100 (MQ=255)
                   cTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTg                                                                             <  2:4118585/100‑1 (MQ=255)
                     tACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgcc                                                                           <  1:3668354/100‑1 (MQ=255)
                     tACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgcc                                                                           <  2:2993905/100‑1 (MQ=255)
                     tACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgcc                                                                           >  2:161173/1‑100 (MQ=255)
                      aCCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccg                                                                          <  2:1349072/100‑1 (MQ=255)
                       ccGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccgc                                                                         >  1:593270/1‑100 (MQ=255)
                       ccGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccgc                                                                         <  2:517472/100‑1 (MQ=255)
                        cGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccgcc                                                                        <  1:2147689/100‑1 (MQ=255)
                         ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                       >  1:3604890/1‑100 (MQ=255)
                         ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                       <  1:1904108/100‑1 (MQ=255)
                         ggTCTGTGCGACCGTTTTCGTGATTTTTTTTTTGTTGTTATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                       >  1:2482715/1‑100 (MQ=255)
                           tCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATc                                                                     >  2:2532342/1‑100 (MQ=255)
                            cTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCa                                                                    >  1:1725229/1‑100 (MQ=255)
                            cTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCa                                                                    >  1:899288/1‑100 (MQ=255)
                              gtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAcc                                                                  >  2:1466535/1‑100 (MQ=255)
                              gtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCCCCCTCAcc                                                                  >  1:106122/1‑100 (MQ=255)
                               tgCGCCCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAccc                                                                 <  2:1734825/100‑1 (MQ=255)
                               tgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACCGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGCTTGCCGCCATCAccc                                                                 >  1:32508/1‑100 (MQ=255)
                               atctaCCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAccc                                                                 <  1:2652875/96‑1 (MQ=255)
                                gCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCt                                                                >  2:1285162/1‑100 (MQ=255)
                                   aCCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAcccagat                                                             >  2:2652884/1‑96 (MQ=255)
                                   aCCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaa                                                             <  1:490120/100‑1 (MQ=255)
                                   aCCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaa                                                             <  2:2770368/100‑1 (MQ=255)
                                    ccGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaaa                                                            <  2:2663129/100‑1 (MQ=255)
                                         ttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGa                                                       >  2:1920367/1‑100 (MQ=255)
                                          tcatGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTTAAAATGGAt                                                      <  1:4079073/97‑1 (MQ=255)
                                          tCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGAt                                                      >  2:3556353/1‑100 (MQ=255)
                                          tCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGAt                                                      <  2:899291/100‑1 (MQ=255)
                                           cGTGGTTTTTATTTTTTTTTTATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGACATCACCCTGAAAATGGATg                                                     >  1:396278/1‑100 (MQ=255)
                                           cGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                     >  2:3669113/1‑100 (MQ=255)
                                            gTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGa                                                    >  1:1466674/1‑100 (MQ=255)
                                              ggTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAAc                                                  >  2:3581128/1‑100 (MQ=255)
                                              ggTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAAc                                                  <  1:3450350/100‑1 (MQ=255)
                                               gTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACa                                                 <  2:88587/100‑1 (MQ=255)
                                                tttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACaa                                                >  2:3643354/1‑100 (MQ=255)
                                                   ttATCCTGTTTTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggc                                             <  1:1037409/100‑1 (MQ=255)
                                                       ccTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGTCATCACCCTGAAAATGGATGAACAAggctggc                                         <  2:4058256/100‑1 (MQ=255)
                                                            tgtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATg                                    >  2:3273496/1‑100 (MQ=255)
                                                             gtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                   >  2:2330755/1‑100 (MQ=255)
                                                             gtgAACGATGTTGAAACAGCCGGATTTAACGCGAAAAACAATACGCTGCTTGGGATTTCCGCCACCACTCCGAAAATGGATGAACAAAGCTGTCTGa                                      >  1:595164/1‑97 (MQ=255)
                                                              tgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGcc                                  <  1:4106283/100‑1 (MQ=255)
                                                               naTCGATGGTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCATGAAAATGGATGAACAAGGATGGCTGATGCCg                                 >  1:3153893/2‑100 (MQ=255)
                                                               gATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCg                                 <  2:1452973/100‑1 (MQ=255)
                                                                aTCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                <  1:669917/100‑1 (MQ=255)
                                                                aTCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                <  1:1389461/100‑1 (MQ=255)
                                                                 tagATGTTGAAACAAACGGATTAAACGCCAAAAACGATGCGCGGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                   >  1:2484730/3‑96 (MQ=255)
                                                                  cGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGac                              <  2:121727/100‑1 (MQ=255)
                                                                    aTGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGacac                            <  2:3533913/100‑1 (MQ=255)
                                                                    aTGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGac                              <  2:1656373/98‑1 (MQ=255)
                                                                     tGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTAGAGATTGCCGCCATCACCCTGAAAACGGATGAACAAGCCTGGCTGATACCGGac                              >  1:4008411/1‑97 (MQ=255)
                                                                       ttGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcaccac                         <  2:1090342/100‑1 (MQ=255)
                                                                        tGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                        >  2:3875873/1‑100 (MQ=255)
                                                                        tGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                        >  1:42730/1‑100 (MQ=255)
                                                                            aCAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattac                    <  2:330110/100‑1 (MQ=255)
                                                                             cAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattcca                   >  1:13174/1‑100 (MQ=255)
                                                                             cAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattaca                   <  2:513336/100‑1 (MQ=255)
                                                                              aGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAAGGGATGAACAAGGCTGGCTGATGCCGGACACCacattacat                  >  1:4036239/1‑100 (MQ=255)
                                                                               gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                 >  2:516720/1‑100 (MQ=255)
                                                                               gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                 >  2:135738/1‑100 (MQ=255)
                                                                               gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                 <  2:3014374/100‑1 (MQ=255)
                                                                               gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                 >  2:3439144/1‑100 (MQ=255)
                                                                                ccGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACAttt                >  2:2262875/1‑100 (MQ=255)
                                                                                 cGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTc               >  1:2410227/1‑100 (MQ=255)
                                                                                 cGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTc               <  2:3916447/100‑1 (MQ=255)
                                                                                 cGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTc               >  1:3361104/1‑100 (MQ=255)
                                                                                  ggATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTcc              >  2:2781502/1‑100 (MQ=255)
                                                                                   gATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCa             <  1:721452/100‑1 (MQ=255)
                                                                                   gATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCa             <  1:1466531/100‑1 (MQ=255)
                                                                                    aTTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCAc            <  1:752781/100‑1 (MQ=255)
                                                                                     tttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACg           <  2:739806/100‑1 (MQ=255)
                                                                                       tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc         <  2:3135280/100‑1 (MQ=255)
                                                                                       tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc         <  2:1417147/100‑1 (MQ=255)
                                                                                       tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc         <  2:1237398/100‑1 (MQ=255)
                                                                                       tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc         <  1:2949728/100‑1 (MQ=255)
                                                                                        aaCGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCg        <  2:1351406/100‑1 (MQ=255)
                                                                                          cGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGGCGGACACCACATTACATTTCCACGTCGaa      >  2:1290797/1‑100 (MQ=255)
                                                                                          cGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGaa      <  1:2615999/100‑1 (MQ=255)
                                                                                          cGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGaa      >  1:2497796/1‑100 (MQ=255)
                                                                                            ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc    <  2:459736/100‑1 (MQ=255)
                                                                                            ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc    <  2:3421611/100‑1 (MQ=255)
                                                                                            ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc    >  2:2470855/1‑100 (MQ=255)
                                                                                              aaaaACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCAt  >  2:3596869/1‑100 (MQ=255)
                                                                                                  |                                                                                               
GATGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCAT  >  NC_000913/1728346‑1728539

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