Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,003,911 (T)8→7 intergenic (+31/‑80) ycbF → / → pyrD putative chaperone/dihydro‑orotate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,003,9040T.99.4% 826.0 / ‑2.8 179intergenic (+24/‑87)ycbF/pyrDputative chaperone/dihydro‑orotate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (74/104);  total (74/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

Only 100 of 219 total aligned reads displayed.

CATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAT  >  NC_000913/1003812‑1003992
                                                                                             |                                                                                        
cATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGt                                                                                          <  2:480804/93‑1 (MQ=255)
  tGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGttt                                                                                        <  1:4544787/93‑1 (MQ=255)
  tGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGttt                                                                                        <  2:2835303/93‑1 (MQ=255)
  tGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGttt                                                                                        <  2:3478602/93‑1 (MQ=255)
  tGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGttt                                                                                        <  2:3777578/93‑1 (MQ=255)
  tGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGttt                                                                                        <  2:4864582/93‑1 (MQ=255)
  tGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTTAGTGCCAAAAAAGAGGAAAAACTGCCCTCAAACAGGGTTGGGAAAAcgtttcgttt                                                                                        >  1:590642/1‑93 (MQ=255)
   gTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGtttt                                                                                       <  2:780485/93‑1 (MQ=255)
   gTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCAGAAAACGTTTGCGtttt                                                                                       <  1:1620356/93‑1 (MQ=255)
   gTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGTGCAAAGAAGGAGGAAAAACTGCCCCTAAACAGGCTCCGAAAACGTTTGTGgttt                                                                                       >  2:1682307/1‑93 (MQ=255)
    tCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGttttt                                                                                      <  2:1124260/93‑1 (MQ=255)
    tCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGttttt                                                                                      <  1:1326178/93‑1 (MQ=255)
       ccTGGCGGGCGATAAATGATTACGGCGGGGTGTGTGCAAAGAAGGAGCAAAAACTTCCCTGAAACAAGGTCGGGAAACGTTTGCGtttttttt                                                                                   >  2:634538/1‑93 (MQ=255)
        cTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGtttttttgc                                                                                  <  2:6152907/93‑3 (MQ=255)
        cTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGtttttttgc                                                                                  <  2:2194071/93‑3 (MQ=255)
         tGGCGGGCGATAAATGATTACGGCGGGTTTAGTGCAAAAAAGGAGCAAAAACTGCCCTGAAACAGGGTCGGAAAACGTTTGCGttttttttcc                                                                                 >  2:5053545/1‑91 (MQ=255)
         tGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGtttttttgcc                                                                                 <  1:2984126/93‑4 (MQ=255)
         tGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGtttttttgcc                                                                                 <  1:3049478/93‑4 (MQ=255)
         tGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGc                                                                                 <  1:3566358/93‑1 (MQ=255)
          ggcgggcgATAAATGATTACGGCGGGTTGGGGGGAAAGAAGGAGCAAAAACTGCCCTGAAACAAGTTCGGGAAAAGGTTGTGttttttttccg                                                                                >  2:87971/1‑90 (MQ=255)
          ggcgggcgATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGtttttttgccg                                                                                <  1:2937202/93‑5 (MQ=255)
             gggcgATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAg                                                                            <  1:5963491/93‑1 (MQ=255)
                  aTAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCaa                                                                       <  1:1107942/93‑1 (MQ=255)
                   tAAATGATTACGGCGGGTTGAGTGCGAAAAAGGAGCAAAATCTGCCCCGAAACAGGTTCGGGAAAAGGTTTCGTTTTTTTTCCGCAGGGCaaa                                                                      >  2:2426363/1‑92 (MQ=255)
                   tAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAAt                                                                      <  2:3005287/93‑1 (MQ=255)
                     aaTGATTACGGCGGGTTGAGTGCAAATAATGATCAAAATCTGCCCTGAAACATGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTc                                                                    <  2:1807108/93‑1 (MQ=255)
                       tGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAACCAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTccc                                                                  <  1:2547227/93‑1 (MQ=255)
                       tGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTccc                                                                  <  1:3754711/93‑1 (MQ=255)
                       tGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTccc                                                                  <  2:4045347/93‑1 (MQ=255)
                       tGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTccc                                                                  <  1:6710172/93‑1 (MQ=255)
                        gATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCt                                                                 <  2:3427503/93‑1 (MQ=255)
                         aTTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCtt                                                                <  2:954711/93‑1 (MQ=255)
                          ttACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCttt                                                               <  2:5517366/93‑1 (MQ=255)
                          ttACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCttt                                                               <  1:324171/93‑1 (MQ=255)
                            aCGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTg                                                             <  1:3677687/93‑1 (MQ=255)
                              ggcggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGt                                                           <  1:1499884/93‑1 (MQ=255)
                               gcggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTc                                                          <  1:4821454/93‑1 (MQ=255)
                               gcggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTc                                                          >  1:4115760/1‑93 (MQ=255)
                               gcggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTc                                                          <  1:992343/93‑1 (MQ=255)
                                  ggTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTACGa                                                       >  1:5319304/1‑93 (MQ=255)
                                   gTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGaa                                                      >  2:2307965/1‑93 (MQ=255)
                                   gTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGCTCGGAAAACGTTTGCGCTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGaa                                                      >  2:5894655/1‑93 (MQ=255)
                                   gTTGACTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCACGTCAATTCCCTTTTGGTCCGaa                                                      <  1:2315307/93‑1 (MQ=255)
                                    ttGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAAc                                                     <  2:6707552/93‑1 (MQ=255)
                                       aGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCg                                                  <  1:1822874/93‑1 (MQ=255)
                                        gTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGc                                                 <  1:1032729/93‑1 (MQ=255)
                                        gTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGc                                                 >  1:650415/1‑93 (MQ=255)
                                         tGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGca                                                <  2:138990/93‑1 (MQ=255)
                                           cAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCACAGGTCAATTCCCTTTTGGTCCGAACTCGcaca                                              <  1:2678864/93‑1 (MQ=255)
                                            aaaGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACat                                             >  2:5323334/1‑93 (MQ=255)
                                            aaaGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACat                                             >  1:3843419/1‑93 (MQ=255)
                                             aagaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACata                                            >  2:6535822/1‑93 (MQ=255)
                                              agaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACataa                                           <  1:3177280/93‑1 (MQ=255)
                                              agaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACataa                                           >  1:6433867/1‑93 (MQ=255)
                                               gaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTTCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACataat                                          >  2:3565369/1‑93 (MQ=255)
                                               gaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACataat                                          <  2:2916549/93‑1 (MQ=255)
                                                aagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACataata                                         <  1:1869866/93‑1 (MQ=255)
                                                aagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACataata                                         <  2:1663897/93‑1 (MQ=255)
                                                aagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACataata                                         <  1:3479027/93‑1 (MQ=255)
                                                 agGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATAc                                        >  1:3333701/1‑93 (MQ=255)
                                                 agGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATAc                                        <  1:3964565/93‑1 (MQ=255)
                                                 agGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATAc                                        <  2:3279909/93‑1 (MQ=255)
                                                    aGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGcc                                     >  2:4263038/1‑93 (MQ=255)
                                                    aGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGcc                                     <  2:4069026/93‑1 (MQ=255)
                                                       aaaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGccccc                                  >  2:1737607/1‑93 (MQ=255)
                                                       aaaaTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGccccc                                  >  2:1064524/1‑93 (MQ=255)
                                                          aTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGt                               <  1:469353/93‑1 (MQ=255)
                                                          aTCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGt                               <  2:5794235/93‑1 (MQ=255)
                                                           tCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGtt                              <  1:3719903/93‑1 (MQ=255)
                                                           tCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGtt                              <  2:6106506/93‑1 (MQ=255)
                                                            cTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                             >  2:2185225/1‑93 (MQ=255)
                                                            cTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGttt                             <  2:4350040/93‑1 (MQ=255)
                                                             tGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTg                            >  2:2565892/1‑93 (MQ=255)
                                                              gCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGc                           >  2:1550360/1‑93 (MQ=255)
                                                              gCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGc                           >  2:1287829/1‑93 (MQ=255)
                                                               cccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGca                          <  2:6267665/93‑1 (MQ=255)
                                                                ccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGcac                         <  1:5588418/93‑1 (MQ=255)
                                                                 cTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGcaca                        >  2:5241303/1‑93 (MQ=255)
                                                                 cTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGcaca                        >  2:5661261/1‑93 (MQ=255)
                                                                  tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCacac                       >  2:4950365/1‑93 (MQ=255)
                                                                    aaaCAGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCg                     <  2:5166896/93‑1 (MQ=255)
                                                                        aGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGaa                 >  1:4910756/1‑93 (MQ=255)
                                                                        aGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGaa                 >  2:851344/1‑93 (MQ=255)
                                                                        aGGTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGaa                 <  1:6639929/93‑1 (MQ=255)
                                                                           gTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATc               >  2:3373223/1‑93 (MQ=255)
                                                                           gTTCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATc               >  1:4199812/1‑93 (MQ=255)
                                                                            ttCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGGTTGCACACCGGGAATcc              >  1:2655248/1‑93 (MQ=255)
                                                                             tCGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCa             <  2:1459108/93‑1 (MQ=255)
                                                                              cGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAg            <  1:5763398/93‑1 (MQ=255)
                                                                              cGGAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAg            <  1:1754027/93‑1 (MQ=255)
                                                                                gAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGga          <  1:6293679/93‑1 (MQ=255)
                                                                                gAAAACGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGga          >  2:1739316/1‑93 (MQ=255)
                                                                                  aaaCGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga        >  1:2025569/1‑93 (MQ=255)
                                                                                  aaaCGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCACCCGGTTTGCACACCGGGAATCCAagagg        >  2:3198647/1‑92 (MQ=255)
                                                                                   aaCGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGagag       <  2:4527661/93‑1 (MQ=255)
                                                                                   aaCGTTTGCGTTTTTTTGCCGCAGGGCAATTCCCTTTTTGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGagag       >  1:4289918/1‑93 (MQ=255)
                                                                                     cGTTTGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGtt     <  1:3294914/93‑1 (MQ=255)
                                                                                       tttGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCa   >  1:2459620/1‑93 (MQ=255)
                                                                                       tttGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCa   <  1:3796539/93‑1 (MQ=255)
                                                                                       tttGCGTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCa   <  1:1238685/93‑1 (MQ=255)
                                                                                             |                                                                                        
CATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAT  >  NC_000913/1003812‑1003992

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