Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,899,375 A→C I15S (ATT→AGT)  ygcU ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,899,3750AC98.7% 822.7 / ‑0.7 235I15S (ATT→AGT) ygcUputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (105/126);  total (106/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 253 total aligned reads displayed.

AGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACAATTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAAGAGAAAACAG  >  NC_000913/2899278‑2899471
                                                                                                 |                                                                                                
aGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACtt                                                                                                <  1:4990644/100‑1 (MQ=255)
aGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATCACGCGATCTGCACCAACACtt                                                                                                >  2:3897049/1‑100 (MQ=255)
  tataAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTc                                                                                              <  1:1630548/100‑1 (MQ=255)
  tataAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTc                                                                                              >  1:3627281/1‑100 (MQ=255)
    taAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCt                                                                                            >  2:2124204/1‑100 (MQ=255)
     aaaTGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCtt                                                                                           <  2:2184496/100‑1 (MQ=255)
      aaTGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCttt                                                                                          >  1:5454877/1‑100 (MQ=255)
      aaTGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCttt                                                                                          <  1:233483/100‑1 (MQ=255)
      aaTGCCATGAATATCAGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTCATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCttt                                                                                          >  2:3605407/1‑100 (MQ=255)
       aTGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTa                                                                                         <  1:1160255/100‑1 (MQ=255)
       aTGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGCTTTTCTTTAATACTGTTTCATAGGTATTTCCGAGACCTCCACCAAAACTTTCCtttc                                                                                         >  2:2494891/1‑99 (MQ=255)
         gcaaTGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTACTTACGCGATCTGCACCAACACTTTCCTTTAGc                                                                                       <  2:1307351/97‑1 (MQ=255)
         gCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGc                                                                                       <  1:4308311/100‑1 (MQ=255)
         gCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGc                                                                                       <  2:1476472/100‑1 (MQ=255)
          ccATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCt                                                                                      <  1:10381/100‑1 (MQ=255)
            aTGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTgg                                                                                    <  2:4253358/100‑1 (MQ=255)
            aTGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTgg                                                                                    >  1:5083311/1‑100 (MQ=255)
            aTGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTgg                                                                                    <  2:3206281/100‑1 (MQ=255)
                 tatCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGAc                                                                               >  1:4540658/1‑100 (MQ=255)
                  atCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACa                                                                              <  2:4031271/100‑1 (MQ=255)
                  atCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACa                                                                              <  2:3068257/100‑1 (MQ=255)
                  atCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCCGCACCAACACTTTCCTTTAGCTGGTCGACa                                                                              >  1:350158/1‑100 (MQ=255)
                    ccGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAAt                                                                            >  2:3768240/1‑100 (MQ=255)
                     cGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATc                                                                           <  1:1788996/100‑1 (MQ=255)
                     cGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATc                                                                           >  1:3846786/1‑100 (MQ=255)
                         aaTTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTg                                                                       <  1:3279793/100‑1 (MQ=255)
                          aTTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTgc                                                                      >  2:1778608/1‑100 (MQ=255)
                              aCGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCga                                                                  >  1:237925/1‑100 (MQ=255)
                               cGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCgag                                                                 >  1:3617935/1‑100 (MQ=255)
                               cGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCgag                                                                 <  1:1371908/100‑1 (MQ=255)
                                gAAACCGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGCCAATCGCTGCGCgaga                                                                <  2:366838/100‑1 (MQ=255)
                                gAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCgagn                                                                <  2:4516123/100‑2 (MQ=255)
                                gAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCgaga                                                                <  1:1658738/100‑1 (MQ=255)
                                gAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCgaga                                                                <  1:4012489/100‑1 (MQ=255)
                                 aaaaCGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGAt                                                               <  1:4482919/100‑1 (MQ=255)
                                  aaaCGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATa                                                              <  2:2685620/100‑1 (MQ=255)
                                      ggTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGa                                                          >  1:96160/1‑100 (MQ=255)
                                       gTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGAc                                                         >  2:1826119/1‑100 (MQ=255)
                                       gTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGAc                                                         >  1:1313857/1‑100 (MQ=255)
                                         cAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACAt                                                       >  2:440804/1‑100 (MQ=255)
                                         cAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACAt                                                       >  1:5056902/1‑100 (MQ=255)
                                         cAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACAt                                                       <  2:1463240/100‑1 (MQ=255)
                                          aaTACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACAATTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATa                                                      <  1:1996003/100‑1 (MQ=255)
                                            tACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACGAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGt                                                    <  2:518876/100‑1 (MQ=255)
                                               tGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAt                                                 <  2:3450429/100‑1 (MQ=255)
                                                gTTTTTCTTTAATACTGTTTGATCGGGAATTACGCGATCGGCACCTACACTTTCCTTTAGCTGGTCGACAATCGCCGCGCGAGATAAAGACATAGTAAtt                                                <  1:5205733/100‑1 (MQ=255)
                                                 tttttCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttc                                               >  1:904044/1‑100 (MQ=255)
                                                  ttttCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttcc                                              <  1:3011051/100‑1 (MQ=255)
                                                   tttCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttcct                                             >  2:2409848/1‑100 (MQ=255)
                                                     tCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttccttc                                           >  1:4076301/1‑100 (MQ=255)
                                                      cTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttccttcc                                          <  2:3528008/100‑1 (MQ=255)
                                                      cTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttccttcc                                          <  1:4937505/100‑1 (MQ=255)
                                                      cTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttccttcc                                          >  2:4893076/1‑100 (MQ=255)
                                                      cTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttccttcc                                          >  1:4363344/1‑100 (MQ=255)
                                                      cTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAAttccttcc                                          <  2:2959860/100‑1 (MQ=255)
                                                       tttAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATtccttcct                                         >  1:1907987/1‑100 (MQ=255)
                                                         tAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGa                                       >  2:395109/1‑100 (MQ=255)
                                                         tAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGa                                       <  2:5026916/100‑1 (MQ=255)
                                                          aaTACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGat                                      >  2:3282410/1‑100 (MQ=255)
                                                          aaTACTGTTGCACCGGTAATTACGCGATCTGCATCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCTCGAGATAAAGACATAGTAATTCCTTCCTGat                                      <  2:1803942/100‑1 (MQ=255)
                                                            tACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGataa                                    >  2:4198902/1‑100 (MQ=255)
                                                            tACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGataa                                    >  2:4097048/1‑100 (MQ=255)
                                                            tACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGataa                                    <  2:2425079/100‑1 (MQ=255)
                                                            tACTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGataa                                    <  2:4835174/100‑1 (MQ=255)
                                                             aCTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGataat                                   <  2:4486125/100‑1 (MQ=255)
                                                              cTGTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGataata                                  <  1:4873228/100‑1 (MQ=255)
                                                                gTTTCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATaaa                                <  1:529332/100‑1 (MQ=255)
                                                                 tttCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATaaaa                               >  1:3000228/1‑100 (MQ=255)
                                                                  ttCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATaaaaa                              <  1:4055566/100‑1 (MQ=255)
                                                                   tCATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAAt                             <  1:1441681/100‑1 (MQ=255)
                                                                    cATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATg                            <  1:2484275/100‑1 (MQ=255)
                                                                    cATCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATg                            <  1:45830/100‑1 (MQ=255)
                                                                    cATCGGTAATTACGCGATCTGCACCAACACTGTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATg                            <  2:4614681/100‑1 (MQ=255)
                                                                     aTCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATgg                           <  2:628875/100‑1 (MQ=255)
                                                                     aTCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATgg                           <  2:4869926/100‑1 (MQ=255)
                                                                      tCGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGa                          >  1:345553/1‑100 (MQ=255)
                                                                       cGGTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGaa                         <  2:5337378/100‑1 (MQ=255)
                                                                        ggTAATTACGCGATCTGCACCAACTCTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGaat                        >  2:1556404/1‑100 (MQ=255)
                                                                        ggTAATTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGaat                        >  1:3094628/1‑100 (MQ=255)
                                                                           aaTTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGaataat                     >  2:782173/1‑100 (MQ=255)
                                                                            aTTACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATc                    >  2:3197533/1‑100 (MQ=255)
                                                                             ttACGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCa                   >  1:2411972/1‑100 (MQ=255)
                                                                               aCGCGATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATACTAATTCCTTCCTGATAATAAAAATGGAATAATCAtt                 <  2:5152788/100‑1 (MQ=255)
                                                                                cgcgATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATta                <  2:667587/100‑1 (MQ=255)
                                                                                cgcgATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATta                >  2:2008022/1‑100 (MQ=255)
                                                                                 gcgATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATtaa               <  2:2776762/100‑1 (MQ=255)
                                                                                  cgATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATtaat              <  2:887090/100‑1 (MQ=255)
                                                                                   gATCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATtaata             <  1:2925415/100‑1 (MQ=255)
                                                                                   gATCTGCACCAACACTTACCTTTAGCTGGTCGAAAAACGCTGCGCGAGATAAAGACATAGTAATTCCTAAATGAAAAAAAAAATGGAATAATCATtaata             <  2:2307297/100‑1 (MQ=255)
                                                                                   gATCTGAACCAATACTTTCCTTTAGTTGGTCGACACTAGCGGCGCGAGATAAAGACATAGTAATTCCTTCCAGGTATTAAAAATGGCCTAGTCATttata             >  2:530510/1‑100 (MQ=255)
                                                                                    aTCTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATtaataa            <  1:2808191/100‑1 (MQ=255)
                                                                                      cTGCACCAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAaga          <  1:2151045/100‑1 (MQ=255)
                                                                                       tGCACCCACACTTTCCTTTAGCTGGTCGACAATCGCTGCGAGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAagag         <  2:3593832/100‑1 (MQ=255)
                                                                                           ccAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAAGAGaaaa     >  1:2654023/1‑100 (MQ=255)
                                                                                           ccAACACTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAAGAGaaaa     >  2:1956025/1‑100 (MQ=255)
                                                                                            caacacTTTCCTGTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAAGAGAAAAc    <  2:5458216/100‑1 (MQ=255)
                                                                                             aacacTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAAGAGAAAACa   >  2:1596216/1‑100 (MQ=255)
                                                                                             aacacTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAAGAGAAAACa   >  2:2063015/1‑100 (MQ=255)
                                                                                             aacacTTTCCATTAGCTGGCCGACAATCGCTGAACGAGATAAAGACATTGTAATTCATAAAAAAAAAAAAAAATGGAATAATCATTAATAAGAGAAAACa   <  2:1232602/100‑1 (MQ=255)
                                                                                              acacTTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAAGAGAAAACAg  <  1:3705154/100‑1 (MQ=255)
                                                                                                 |                                                                                                
AGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACAATTTCCTTTAGCTGGTCGACAATCGCTGCGCGAGATAAAGACATAGTAATTCCTTCCTGATAATAAAAATGGAATAATCATTAATAAGAGAAAACAG  >  NC_000913/2899278‑2899471

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: