Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,966 T→C 5.7% H559R (CAC→CGC) ‡ cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,9660TC 5.7% 580.9 / 17.2 176H559R (CAC→CGC) ‡cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (81/85);  new base C (2/8);  total (83/93)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

Only 100 of 177 total aligned reads displayed.

CCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGA  >  NC_000913/448818‑449113
                                                                                                                                                    |                                                                                                                                                   
 cAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTg                                                                                                                                                    <  3:821338/149‑1 (MQ=255)
  aGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGa                                                                                                                                                   >  2:661133/1‑149 (MQ=255)
         ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                                                                                                            <  1:1060039/149‑3 (MQ=255)
         ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                                                                                                            <  2:1626702/149‑3 (MQ=255)
         ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                                                                                                            <  2:216357/149‑3 (MQ=255)
         ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                                                                                                            <  3:1509144/149‑3 (MQ=255)
          gCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCgg                                                                                                                                           <  4:1003525/149‑1 (MQ=255)
              aaCCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAAc                                                                                                                                                 <  2:1084455/139‑1 (MQ=255)
              aaCCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAAc                                                                                                                                                 >  1:1084455/1‑139 (MQ=255)
                ccGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTAc                                                                                                                                     >  4:1548124/1‑149 (MQ=255)
                  gAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTga                                                                                                                                             >  4:369356/1‑138 (MQ=255)
                  gAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTga                                                                                                                                             <  3:369356/139‑2 (MQ=255)
                      aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCaaa                                                                                                                               <  1:1542126/149‑1 (MQ=255)
                      aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCaaa                                                                                                                               <  2:1512407/149‑1 (MQ=255)
                         gtggAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGtt                                                                                                                            <  3:173842/149‑1 (MQ=255)
                         gtggAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGtt                                                                                                                            >  3:973255/1‑149 (MQ=255)
                           ggAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTtat                                                                                                                          <  3:908632/149‑1 (MQ=255)
                            gagaAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTtata                                                                                                                         >  1:71196/1‑149 (MQ=255)
                               aaaGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGaa                                                                                                                      <  3:347673/149‑1 (MQ=255)
                                aaGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAAc                                                                                                                     <  3:436037/149‑1 (MQ=255)
                                 aGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACg                                                                                                                    <  4:352447/149‑1 (MQ=255)
                                 aGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACg                                                                                                                    <  4:1328035/149‑1 (MQ=255)
                                   cTGCGATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGa                                                                                                                  <  4:424092/149‑1 (MQ=255)
                                   cTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGa                                                                                                                  <  1:279151/149‑1 (MQ=255)
                                        aTGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCgg                                                                                                             <  3:471216/149‑1 (MQ=255)
                                            cGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGaa                                                                                                         <  4:772799/149‑1 (MQ=255)
                                                aCCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAgg                                                                                                     <  2:916503/149‑1 (MQ=255)
                                                  cTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGtt                                                                                                   >  4:588362/1‑149 (MQ=255)
                                                    gCACCGCTGTTTTTCGGCATATGAATTTCTTCATCGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGc                                                                                                 <  3:1012747/149‑1 (MQ=255)
                                                    gCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGc                                                                                                 >  3:464710/1‑149 (MQ=255)
                                                      aCCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGccc                                                                                               >  3:640155/1‑149 (MQ=255)
                                                        cGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCAc                                                                                             >  4:1629410/1‑149 (MQ=255)
                                                        cGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCAc                                                                                             >  2:1642464/1‑149 (MQ=255)
                                                         gCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACt                                                                                            >  2:742575/1‑149 (MQ=255)
                                                            gTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTTCGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCCCTGCa                                                                                         >  1:1315223/1‑149 (MQ=255)
                                                            gTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCa                                                                                         <  1:258693/149‑1 (MQ=255)
                                                                 tCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTa                                                                                    <  2:1210479/149‑1 (MQ=255)
                                                                  cGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTAc                                                                                   <  1:1219408/149‑1 (MQ=255)
                                                                   ggCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACg                                                                                  >  1:1385491/1‑149 (MQ=255)
                                                                   ggCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCATGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACg                                                                                  <  4:422023/149‑1 (MQ=255)
                                                                     cATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGc                                                                                <  3:1126189/149‑1 (MQ=255)
                                                                      atatGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGcc                                                                               <  4:326816/149‑1 (MQ=255)
                                                                        atGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTTGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCg                                                                                                              >  3:1330913/1‑116 (MQ=255)
                                                                         tGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCgggg                                                                                                           <  4:749491/118‑1 (MQ=255)
                                                                         tGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCgggg                                                                                                           >  3:749491/1‑118 (MQ=255)
                                                                          gAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCAccc                                                                           >  2:1041574/1‑149 (MQ=255)
                                                                               tctttATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACgg                                                                      <  4:1461770/149‑1 (MQ=255)
                                                                               tcttcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACgg                                                                      <  1:181331/149‑1 (MQ=255)
                                                                               tcttcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACgg                                                                      <  4:397995/149‑1 (MQ=255)
                                                                               tcttcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACgg                                                                      <  3:1561368/149‑1 (MQ=255)
                                                                                cttcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACggg                                                                     <  3:487237/149‑1 (MQ=255)
                                                                                cttcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACggg                                                                     >  2:1528333/1‑149 (MQ=255)
                                                                                  tcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCa                                                                                         <  3:385239/127‑1 (MQ=255)
                                                                                  tcATAGTGGGCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTTCTGTACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTc                                                                   >  2:452273/1‑149 (MQ=255)
                                                                                    aTAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGc                                                                 <  2:337596/149‑1 (MQ=255)
                                                                                      aGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCggg                                                                                                            <  3:726013/104‑1 (MQ=255)
                                                                                      aGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCggg                                                                                                            >  4:726013/1‑104 (MQ=255)
                                                                                        tGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGt                                                             <  2:1215617/149‑1 (MQ=255)
                                                                                         ggTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTc                                                            <  3:1629410/149‑1 (MQ=255)
                                                                                           tCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACgg                                                                      >  1:574443/1‑137 (MQ=255)
                                                                                           tCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACgg                                                                      <  2:574443/137‑1 (MQ=255)
                                                                                                cTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCAc                                                     >  4:492113/1‑149 (MQ=255)
                                                                                                 tttttATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACg                                                    >  2:1118040/1‑149 (MQ=255)
                                                                                                    ttATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGtt                                                 <  1:742575/149‑1 (MQ=255)
                                                                                                      ataCGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCt                                               <  1:482020/149‑1 (MQ=255)
                                                                                                       taCGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTg                                              <  2:1385491/149‑1 (MQ=255)
                                                                                                          gCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtc                                           >  2:355974/1‑149 (MQ=255)
                                                                                                            ttCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgc                                         >  2:859703/1‑149 (MQ=255)
                                                                                                               gCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTc                                                                   >  3:1065190/1‑120 (MQ=255)
                                                                                                               gCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTc                                                                   <  4:1065190/120‑1 (MQ=255)
                                                                                                                 cTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGATTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgcggtcg                                    >  4:1349170/1‑149 (MQ=255)
                                                                                                                  tttCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgcggtcgc                                   >  1:879021/1‑149 (MQ=255)
                                                                                                                  tttCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgcggtcgc                                   >  2:774823/1‑149 (MQ=255)
                                                                                                                  tttCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgcggtcgc                                   <  1:1008081/149‑1 (MQ=255)
                                                                                                                  tttCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgcggtcgc                                   >  4:122457/1‑149 (MQ=255)
                                                                                                                     cactTTCATTACCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGaa                                <  4:687291/147‑1 (MQ=255)
                                                                                                                       ctTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATa                              <  2:667468/149‑1 (MQ=255)
                                                                                                                        tttCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAg                             >  3:53557/1‑149 (MQ=255)
                                                                                                                        tttCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAg                             >  3:25774/1‑149 (MQ=255)
                                                                                                                        tttCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAg                             >  3:1670100/1‑149 (MQ=255)
                                                                                                                          tCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGaa                           >  4:546551/1‑149 (MQ=255)
                                                                                                                           cATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGaaa                          <  3:780997/149‑1 (MQ=255)
                                                                                                                           cATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGaaa                          <  1:661133/149‑1 (MQ=255)
                                                                                                                            aTTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTc                                                <  4:1371404/126‑1 (MQ=255)
                                                                                                                            aTTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAAc                         >  1:1451366/1‑149 (MQ=255)
                                                                                                                            aTTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAAc                         <  4:197718/149‑1 (MQ=255)
                                                                                                                             tttCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACg                        >  1:358441/1‑149 (MQ=255)
                                                                                                                             tttCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACg                        >  1:1595875/1‑149 (MQ=255)
                                                                                                                             tttCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACg                        >  3:784186/1‑149 (MQ=255)
                                                                                                                               tCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTa                      >  3:407811/1‑149 (MQ=255)
                                                                                                                                  cAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACAt                   >  4:884080/1‑149 (MQ=255)
                                                                                                                                    gAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCt                 >  2:1315241/1‑149 (MQ=255)
                                                                                                                                           tCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc          >  3:610057/1‑149 (MQ=255)
                                                                                                                                           tCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc          >  3:1619424/1‑149 (MQ=255)
                                                                                                                                           tCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc          >  1:990594/1‑149 (MQ=255)
                                                                                                                                             aCGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGa        >  4:81750/1‑149 (MQ=255)
                                                                                                                                                 tCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa    <  1:1642464/149‑1 (MQ=255)
                                                                                                                                                   gTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTAc                                                                                   <  3:1396435/68‑1 (MQ=255)
                                                                                                                                                   gTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCaga  <  4:1670100/149‑1 (MQ=255)
                                                                                                                                                   gTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCaga  >  3:167470/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGA  >  NC_000913/448818‑449113

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