Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,286,275 A→G intergenic (+68/+35) b1228 → / ← tpr hypothetical protein/protamine‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,2750AG87.9% 548.0 / 55.8 207intergenic (+68/+35)b1228/tprhypothetical protein/protamine‑like protein
Reads supporting (aligned to +/- strand):  ref base A (13/12);  new base G (100/82);  total (113/94)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 211 total aligned reads displayed.

CTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCT  >  NC_000913/1286185‑1286364
                                                                                          |                                                                                         
cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACt                                                                                         >  1:6625360/1‑93 (MQ=255)
cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAAAGGTCCGGGACt                                                                                         >  2:6348004/1‑93 (MQ=255)
 tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCtt                                                                                        <  2:1665561/93‑1 (MQ=255)
  gACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                                                                                       <  2:2855242/93‑1 (MQ=255)
   gcGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCtttt                                                                                      <  1:5769345/92‑1 (MQ=255)
    cGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                                                                     >  1:5991266/1‑93 (MQ=255)
     gCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCATAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGa                                                                                    >  2:3803844/1‑93 (MQ=255)
     gCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGCGTTCAGGGGCTTTTGa                                                                                    >  1:5597252/1‑93 (MQ=255)
      cTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGaa                                                                                   <  2:2286891/93‑1 (MQ=255)
        cGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAg                                                                                 >  2:6051775/1‑93 (MQ=255)
        cGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAg                                                                                 >  2:5114532/1‑93 (MQ=255)
            ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAt                                                                             <  1:975440/93‑1 (MQ=255)
            ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAt                                                                             <  2:6357529/93‑1 (MQ=255)
             cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGTTTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtg                                                                            >  2:1327823/1‑93 (MQ=255)
             cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtg                                                                            <  1:638426/93‑1 (MQ=255)
             cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtg                                                                            <  1:4555041/93‑1 (MQ=255)
             cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtg                                                                            <  1:4130397/93‑1 (MQ=255)
             cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtg                                                                            <  1:2224826/93‑1 (MQ=255)
             cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtg                                                                            >  1:1025379/1‑93 (MQ=255)
              gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtgg                                                                           >  2:1524875/1‑93 (MQ=255)
              gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtgg                                                                           >  2:3957366/1‑93 (MQ=255)
              gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtgg                                                                           >  2:4644517/1‑93 (MQ=255)
               aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtggt                                                                          <  2:1466818/93‑1 (MQ=255)
               aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtggt                                                                          <  2:6139818/93‑1 (MQ=255)
               aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggt                                                                          >  2:743227/1‑93 (MQ=255)
               aaCCTTAGTCGAAGCTTATCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGCGAGTTCAGGGGCTTTTGGAAGTGAtggg                                                                          >  1:1075473/1‑92 (MQ=14)
                 ccTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtggtgg                                                                        >  2:4856390/1‑93 (MQ=255)
                  cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtggtggt                                                                       <  2:5450260/93‑1 (MQ=255)
                   ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtggtggtg                                                                      >  2:233942/1‑93 (MQ=255)
                    tAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGAtggtggtgg                                                                     >  1:1854822/1‑93 (MQ=255)
                     aGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTggg                                                                    >  2:2998681/1‑93 (MQ=255)
                      gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTgggg                                                                   >  2:4914738/1‑93 (MQ=255)
                      gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTgggg                                                                   <  1:5056267/93‑1 (MQ=255)
                          aaGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAg                                                               <  2:5367887/93‑1 (MQ=255)
                          aaGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAg                                                               <  1:5810575/93‑1 (MQ=255)
                          aaGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAg                                                               >  1:191199/1‑93 (MQ=255)
                            gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                             <  2:5780693/93‑1 (MQ=255)
                            gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                             >  1:6584573/1‑93 (MQ=255)
                            gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                             <  1:2323673/93‑1 (MQ=255)
                             cTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                            >  1:6330432/1‑93 (MQ=255)
                              ttCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                                           >  1:6218981/1‑93 (MQ=255)
                              ttCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                                           >  1:2901498/1‑93 (MQ=255)
                                ctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAc                                                         >  2:1736035/1‑93 (MQ=255)
                                 tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                                        >  1:3881298/1‑93 (MQ=255)
                                  cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                       >  1:5452022/1‑93 (MQ=255)
                                  cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                       <  2:6314152/93‑1 (MQ=255)
                                  cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                       <  2:638534/93‑1 (MQ=255)
                                   aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                                                      >  2:2485332/1‑93 (MQ=255)
                                   aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                                                      <  1:819302/93‑1 (MQ=255)
                                   aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                                                      >  1:1669689/1‑93 (MQ=255)
                                    tCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAGTACTCAg                                                     >  1:3271959/1‑93 (MQ=255)
                                       ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                  <  2:4298402/93‑1 (MQ=255)
                                       ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                  >  2:3236703/1‑93 (MQ=255)
                                          cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgc                                               <  2:1049179/93‑1 (MQ=255)
                                          cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgc                                               >  2:1753355/1‑93 (MQ=255)
                                          cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgc                                               <  2:3042121/93‑1 (MQ=255)
                                          cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgc                                               >  2:5862593/1‑93 (MQ=255)
                                           ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcg                                              >  2:5534235/1‑93 (MQ=255)
                                             gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg                                            >  1:1631769/1‑93 (MQ=255)
                                             gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg                                            <  2:3217951/93‑1 (MQ=255)
                                             gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg                                            >  1:2866380/1‑93 (MQ=255)
                                              cATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGATAGTGATGGAGGTGGGGGAATGATTACTCATCGCTgcgcgc                                           >  2:6651856/1‑93 (MQ=14)
                                               aTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt                                          >  1:249481/1‑93 (MQ=255)
                                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGGAGGATTACTCAGCGCTGCGCGCTTCg                                       >  1:4859544/1‑93 (MQ=255)
                                                   gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc                                      <  1:162010/93‑1 (MQ=37)
                                                     aGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc                                    <  2:2015476/93‑1 (MQ=255)
                                                      taaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTCTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt                                   <  1:2743797/92‑1 (MQ=255)
                                                      gAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt                                   >  2:565100/1‑93 (MQ=255)
                                                      gAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt                                   >  1:4662257/1‑93 (MQ=255)
                                                      gAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTAGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt                                   <  1:5902288/93‑1 (MQ=255)
                                                       aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCtt                                  >  1:4002503/1‑93 (MQ=255)
                                                       aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCtt                                  <  1:1243786/93‑1 (MQ=37)
                                                        atatTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                 >  1:3025358/1‑93 (MQ=255)
                                                        atatTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                 <  2:1761602/93‑1 (MQ=37)
                                                         tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                <  1:3542009/93‑1 (MQ=37)
                                                         tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                >  1:1435863/1‑93 (MQ=255)
                                                         tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                >  2:5333492/1‑93 (MQ=255)
                                                         tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                <  2:3112516/93‑1 (MQ=37)
                                                          atTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg                               >  1:2364475/1‑93 (MQ=255)
                                                           tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCggg                              >  1:5875345/1‑93 (MQ=255)
                                                            ttGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGAGGTGGGGGAAGGATTACTCAACGCTGCGCGATTCGCACTTCGGGt                             >  2:2650645/1‑93 (MQ=11)
                                                             tGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTc                            <  1:3938567/93‑1 (MQ=37)
                                                               aTTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGt                          >  2:3659540/1‑93 (MQ=255)
                                                               aTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGt                          <  2:6478717/93‑1 (MQ=255)
                                                               aTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGt                          <  2:213093/93‑1 (MQ=255)
                                                                   cGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                      <  1:4329486/93‑1 (MQ=37)
                                                                    ggATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt                                   <  1:5213092/79‑1 (MQ=255)
                                                                    ggATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                     >  2:5568670/1‑93 (MQ=25)
                                                                     gATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                    >  2:5687697/1‑93 (MQ=255)
                                                                         cGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGc                <  2:185615/93‑1 (MQ=37)
                                                                           cTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCaa              >  2:4503691/1‑93 (MQ=255)
                                                                           cTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa              >  1:2103413/1‑93 (MQ=21)
                                                                           cTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa              <  2:2064919/93‑1 (MQ=21)
                                                                             tGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACg            >  1:5342881/1‑93 (MQ=18)
                                                                              gagaGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc           <  2:5991266/93‑1 (MQ=40)
                                                                               agagTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGt                             <  2:626279/74‑1 (MQ=21)
                                                                               agagTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGct          >  1:376997/1‑93 (MQ=21)
                                                                                 agTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCtt                                  <  2:1199734/67‑1 (MQ=18)
                                                                                  gTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc       <  2:6596518/93‑1 (MQ=14)
                                                                                  gTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTAGTTGCCTGCGGCAACGctctc       >  2:1947046/1‑93 (MQ=9)
                                                                                          |                                                                                         
CTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCT  >  NC_000913/1286185‑1286364

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