Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,513,653 C→T 5.8% S204S (AGC→AGT nupC → nucleoside (except guanosine) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,513,6530CT 5.8% 388.0 / 6.8 120S204S (AGC→AGTnupCnucleoside (except guanosine) transporter
Reads supporting (aligned to +/- strand):  ref base C (57/56);  new base T (5/2);  total (62/58)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

Only 100 of 121 total aligned reads displayed.

AAGATATCCTCGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGCCATTATCGTTGC  >  NC_000913/2513505‑2513802
                                                                                                                                                    |                                                                                                                                                     
aaGATATCCTCGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATTCATGACCATGCTGGATCCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCa                                                                                                                                                      <  2:1598049/150‑1 (MQ=255)
        cTCGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATc                                                                                                                                              <  1:1328710/150‑1 (MQ=255)
        cTCGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATc                                                                                                                                              >  2:916532/1‑150 (MQ=255)
         tCGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTGCTGAACATGTTCAGTAccagcgcag                                                                                                                                             <  2:1626165/150‑8 (MQ=255)
         tCGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTGCTGAACATGTTCAGTAccagcgcag                                                                                                                                             <  2:1590392/150‑8 (MQ=255)
          cGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGt                                                                                                                                            >  2:675158/1‑150 (MQ=255)
          cGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGt                                                                                                                                            <  1:1421213/150‑1 (MQ=255)
          cGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGt                                                                                                                                            <  1:1414636/150‑1 (MQ=255)
          cGGCAAAATCTCCCGTAATCGTATGTAAACCATGGCTGCGACGGCAATTTCCACCGTGTCTATTTCCATCTTTGTTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGt                                                                                                                                            <  2:259569/150‑1 (MQ=255)
           ggCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGTAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTg                                                                                                                                           <  1:1551368/150‑1 (MQ=255)
           ggCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTg                                                                                                                                           <  1:354822/150‑1 (MQ=255)
            gCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTTGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGc                                                                                                                                          <  2:927270/150‑1 (MQ=255)
            gCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCCGCACCTTTATCGTGc                                                                                                                                          >  2:1392213/1‑150 (MQ=255)
            gCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGGCCATCGTTGGTGCATACATGACCATGCTGGAACAGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGc                                                                                                                                          >  1:98223/1‑150 (MQ=255)
              aaaaTCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTg                                                                                                                                        <  1:408888/150‑1 (MQ=255)
               aaaTCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGt                                                                                                                                       >  2:1301602/1‑150 (MQ=255)
                 aTCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCg                                                                                                                                     >  2:1231950/1‑150 (MQ=255)
                     cccGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGa                                                                                                                                 >  2:1300731/1‑150 (MQ=255)
                           aTCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATc                                                                                                                           <  2:466700/150‑1 (MQ=255)
                           aTCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATc                                                                                                                           >  2:82572/1‑150 (MQ=255)
                            tCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATcc                                                                                                                          <  1:538687/150‑1 (MQ=255)
                             cGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGATCCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCt                                                                                                                         <  1:327255/150‑1 (MQ=255)
                                 tGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTAcc                                                                                                                     >  1:1078419/1‑150 (MQ=255)
                                  gTACACAATGGCTGCCACGGCAAGGGCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCg                                                                                                                    <  2:275124/150‑1 (MQ=255)
                                     cacCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCgtgt                                                                                                                 >  1:1495237/1‑150 (MQ=255)
                                        cATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTTGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGa                                                                                                              <  1:1467550/150‑1 (MQ=255)
                                            gCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGcc                                                                                                          <  2:1104489/150‑1 (MQ=255)
                                               gCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGt                                                                                                       >  2:18409/1‑150 (MQ=255)
                                               gCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGt                                                                                                       >  2:1715748/1‑150 (MQ=255)
                                                 cACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTga                                                                                                     >  1:1024983/1‑150 (MQ=255)
                                                   cGGCAATGTCCAGCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTgaag                                                                                                   <  1:362398/150‑1 (MQ=255)
                                                    ggCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTgaaga                                                                                                  >  1:1117027/1‑150 (MQ=255)
                                                     gCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGtt                                                                                                                >  2:797616/1‑135 (MQ=255)
                                                     gCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGtt                                                                                                                <  1:797616/135‑1 (MQ=255)
                                                            ccACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTTGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATcc                                                                                          <  1:1043874/150‑1 (MQ=255)
                                                             cACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCa                                                                                         <  2:659377/150‑1 (MQ=255)
                                                               ccGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGa                                                                                       <  1:82572/150‑1 (MQ=255)
                                                                cGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGAt                                                                                      <  2:1364502/150‑1 (MQ=255)
                                                                   gtCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTc                                                                                   <  2:1600835/150‑1 (MQ=255)
                                                                   gtCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTc                                                                                   >  1:1526238/1‑150 (MQ=255)
                                                                    tCGATGTCCATCGTTGGTGCATACATGACCAGGCTGGAACCGAAATACGTCGTTTCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTcc                                                                                  <  1:1319147/150‑1 (MQ=255)
                                                                     cGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCa                                                                                 >  1:886155/1‑150 (MQ=255)
                                                                      gATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGCTCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCaa                                                                                >  1:1139588/1‑150 (MQ=255)
                                                                        tGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATAGGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAAcc                                                                              <  2:294434/150‑1 (MQ=255)
                                                                        tGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAg                                                                                        <  1:983854/140‑1 (MQ=255)
                                                                        tGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAAcc                                                                              >  2:1021095/1‑150 (MQ=255)
                                                                         gTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCt                                                                             >  2:1393153/1‑150 (MQ=255)
                                                                          tCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTg                                                                            <  1:1583267/150‑1 (MQ=255)
                                                                           ccATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGc                                                                           <  1:916532/150‑1 (MQ=255)
                                                                                  tGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAgg                                                                    <  2:1549200/150‑1 (MQ=255)
                                                                                  tGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAgg                                                                    <  1:18409/150‑1 (MQ=255)
                                                                                   ggTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGt                                                                   >  2:736774/1‑150 (MQ=255)
                                                                                    gTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTc                                                                  >  1:1353932/1‑150 (MQ=255)
                                                                                     tGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCa                                                                 <  1:675158/150‑1 (MQ=255)
                                                                                     tGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGGCGCTGATCAATCCTTCCCGTGTGGATGCCAGTGAAGAAAACCTCCAGATGTCCAACCTGCACGCAGGTCa                                                                 >  2:272736/1‑150 (MQ=255)
                                                                                      gCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCag                                                                >  2:200713/1‑150 (MQ=255)
                                                                                       cATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCaga                                                               <  2:1495237/150‑1 (MQ=255)
                                                                                          aCATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGct                                                            >  2:729221/1‑150 (MQ=255)
                                                                                          aCATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGct                                                            <  1:1231950/150‑1 (MQ=255)
                                                                                               aCCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCttcttc                                                       <  2:1117027/150‑1 (MQ=255)
                                                                                                ccATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCg                                                      >  2:773454/1‑150 (MQ=255)
                                                                                                ccATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCg                                                      <  2:831363/150‑1 (MQ=255)
                                                                                                 cATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACCTCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGa                                                     >  1:266163/1‑150 (MQ=255)
                                                                                                      tGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGc                                                >  1:569066/1‑150 (MQ=255)
                                                                                                      tGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGc                                                <  1:1376790/150‑1 (MQ=255)
                                                                                                        gAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTg                                              >  1:1312939/1‑150 (MQ=255)
                                                                                                        gAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTg                                              <  1:1420933/150‑1 (MQ=255)
                                                                                                         aaCCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTgg                                             <  2:1477953/150‑1 (MQ=255)
                                                                                                          aCCGAAATATGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTggg                                            <  1:1696994/150‑1 (MQ=255)
                                                                                                          aCCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTggg                                            <  1:1715748/150‑1 (MQ=255)
                                                                                                            cGAAATACGTCGTTGCTGCGCTTGTACTTAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTg                                          <  1:569707/150‑1 (MQ=255)
                                                                                                            cGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTg                                          >  2:1364031/1‑150 (MQ=255)
                                                                                                            cGAAATACGGCGTTGCTGCGCTGGTACTGACCATGTTCAGCAGCTTTATCGGGCTGTCGCTGATCTATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCTGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTg                                          <  2:818745/150‑1 (MQ=255)
                                                                                                             gAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGCTGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGa                                         >  1:1082724/1‑150 (MQ=255)
                                                                                                             gAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACCAAGGTCAGAGCTTCTTCGAAATGATGGGTGa                                         >  1:1064822/1‑150 (MQ=255)
                                                                                                                     tcgtTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTc                                 >  1:879572/1‑150 (MQ=255)
                                                                                                                      cgtTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCt                                <  1:1139340/150‑1 (MQ=255)
                                                                                                                      cgtTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTAACCGTGTTGATGCCAGTGAAGAAAACATCCAGAAGTCCAACCTGCACGAAGGTCAGACCTTCTACGAAATGCAGGGTGAATACAATca                                >  2:1655767/1‑149 (MQ=255)
                                                                                                                         tgctgcGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGc                             <  1:736774/150‑1 (MQ=255)
                                                                                                                          gctgcGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTTATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCa                            <  1:897867/150‑1 (MQ=255)
                                                                                                                          gctgcGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCa                            >  2:1361323/1‑150 (MQ=255)
                                                                                                                          gctgcGCTGGTACGGAACATGTTCAGCACCTTTATCGTGCTGTCGCGGATCAATCCTTACCGTGTTGATGCCAGTGAATAAAACATCCAGATGTCCAACCTGCCCGAAGGTAAGAGCTTATTCGAAATGCGGGGTGAATAAATTCTGGCa                            >  2:1286898/1‑150 (MQ=255)
                                                                                                                           cgataaaggTGCTGAACATGTTCAGTACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCACCCTGCACGAAGGTCAGAGCTTCTTAGAAATGCTGGGTGAATACTTTCTGGCAg                           >  2:394788/8‑150 (MQ=255)
                                                                                                                           cgataaaggTGCTGAACATGTTCAGTACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAg                           >  2:1004729/8‑150 (MQ=255)
                                                                                                                           cgataaaggTGCTGAACATGTTCAGTACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGCGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAg                           >  1:802535/8‑150 (MQ=255)
                                                                                                                                cTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTc                      <  1:644049/150‑1 (MQ=255)
                                                                                                                                 aggTGCTGAACATGTTCAGTACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAAACTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCa                     >  2:436715/2‑150 (MQ=255)
                                                                                                                                   gTACTGACCATGTTCAGCACCTTTATCGTGCTGGCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCaaa                   <  2:98223/150‑1 (MQ=255)
                                                                                                                                      cTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGtt                >  2:893617/1‑150 (MQ=255)
                                                                                                                                      cTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACCAAGGTCAGAGCTTCTTCCAAATGCTGGGTGAATACATACTGGCAGGTTTCAAAGtt                >  1:346152/1‑150 (MQ=255)
                                                                                                                                       tGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTg               <  1:1393153/150‑1 (MQ=255)
                                                                                                                                        gAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGc              >  2:1280099/1‑150 (MQ=255)
                                                                                                                                        gAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGc              <  2:1078419/150‑1 (MQ=255)
                                                                                                                                        gAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGc              >  1:1711015/1‑150 (MQ=255)
                                                                                                                                               ttCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGCCATTATc       >  1:896371/1‑150 (MQ=255)
                                                                                                                                               ggccgCACCGTTAGCGGGCTGTCGCTGCTCAATCCTTACCGTGGTGATGCCAGTGAAGACAACATCCAGATGTCCAACCTGCACGACGGTCAGAGCTTCTTCGACCTGCTGGGTGAATACATTCTGGCAGGTTTCACAGTTGCCATTagc       <  1:1153915/146‑3 (MQ=255)
                                                                                                                                                tCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGCCATTATCg      >  2:1445810/1‑150 (MQ=255)
                                                                                                                                                tCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGCCATTATCg      >  2:1395796/1‑150 (MQ=255)
                                                                                                                                                   gCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGCCATTATCGTTg   >  2:1216326/1‑150 (MQ=255)
                                                                                                                                                    cACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGCCATTATCGTTGc  <  2:1526238/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGATATCCTCGGCAAAATCTCCCGTAATCGTATGTACACCATGGCTGCCACGGCAATGTCCACCGTGTCGATGTCCATCGTTGGTGCATACATGACCATGCTGGAACCGAAATACGTCGTTGCTGCGCTGGTACTGAACATGTTCAGCACCTTTATCGTGCTGTCGCTGATCAATCCTTACCGTGTTGATGCCAGTGAAGAAAACATCCAGATGTCCAACCTGCACGAAGGTCAGAGCTTCTTCGAAATGCTGGGTGAATACATTCTGGCAGGTTTCAAAGTTGCCATTATCGTTGC  >  NC_000913/2513505‑2513802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: