Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,240,380 T→C V19A (GTA→GCA)  psiE → phosphate starvation inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,240,3800TC98.2% 397.7 / ‑1.6 108V19A (GTA→GCA) psiEphosphate starvation inducible protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (62/44);  total (63/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 108 total aligned reads displayed.

TCCAGGTTGAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTATTT  >  NC_000913/4240233‑4240527
                                                                                                                                                   |                                                                                                                                                   
tCCAGGTTGAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACt                                                                                                                                                  >  2:1690482/1‑151 (MQ=255)
  cAGGTTGAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCa                                                                                                                                                <  2:1047714/151‑1 (MQ=255)
   aGGTTGAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCaa                                                                                                                                               <  1:1706376/151‑1 (MQ=255)
    ggTTGAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAAt                                                                                                                                              <  2:1520785/151‑1 (MQ=255)
    ggTTGAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAAt                                                                                                                                              >  1:1208270/1‑151 (MQ=255)
      ttGAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCt                                                                                                                                            >  1:49705/1‑151 (MQ=255)
        gAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCtt                                                                                                                                           <  1:495581/150‑1 (MQ=255)
        gAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTg                                                                                                                                          <  1:123800/151‑1 (MQ=255)
         aaCAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTgg                                                                                                                                         <  2:603995/151‑1 (MQ=255)
          aCAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGc                                                                                                                                        >  1:341038/1‑151 (MQ=255)
           cAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGcc                                                                                                                                       >  1:306373/1‑151 (MQ=255)
           cAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGcc                                                                                                                                       <  2:1343888/151‑1 (MQ=255)
             aaaCATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCctg                                                                                                                                     >  1:77793/1‑151 (MQ=255)
               aCATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCctgct                                                                                                                                   <  2:1357885/151‑1 (MQ=255)
               aCATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCctgct                                                                                                                                   >  2:844187/1‑151 (MQ=255)
                   acacAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGc                                                                                                                               <  1:1456132/151‑1 (MQ=255)
                     acaAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCt                                                                                                                             <  1:937369/151‑1 (MQ=255)
                     acaAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCt                                                                                                                             <  2:1552323/151‑1 (MQ=255)
                      caAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTc                                                                                                                            <  1:167763/151‑1 (MQ=255)
                       aaaaaaTATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCg                                                                                                                           <  2:365600/151‑1 (MQ=255)
                          aaaTATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGcc                                                                                                                        >  1:706904/1‑151 (MQ=255)
                            atatAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTg                                                                                                                      <  1:851771/151‑1 (MQ=255)
                              ataGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGAt                                                                                                                    >  1:1555235/1‑151 (MQ=255)
                               taGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGAtt                                                                                                                   >  2:846265/1‑151 (MQ=255)
                                aGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGAttt                                                                                                                  >  1:1347810/1‑151 (MQ=255)
                                  aTCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTg                                                                                                                >  1:1493780/1‑151 (MQ=255)
                                   tctcCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTgg                                                                                                               >  1:1474137/1‑151 (MQ=255)
                                    ctcCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGt                                                                                                              >  2:618361/1‑151 (MQ=255)
                                      ccGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTg                                                                                                            <  2:561944/151‑1 (MQ=255)
                                       cGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGt                                                                                                           >  1:1404067/1‑151 (MQ=255)
                                        gTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTc                                                                                                          <  1:771862/151‑1 (MQ=255)
                                         tCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCGATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCt                                                                                                         >  2:840604/1‑151 (MQ=255)
                                          cacaTTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCtt                                                                                                        <  2:864477/151‑1 (MQ=255)
                                           acaTTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTc                                                                                                       >  1:268528/1‑151 (MQ=255)
                                             aTTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCt                                                                                                     >  1:287973/1‑151 (MQ=255)
                                             aTTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCt                                                                                                     <  2:842499/151‑1 (MQ=255)
                                             aTTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCt                                                                                                     >  2:499027/1‑151 (MQ=255)
                                              tttttGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTc                                                                                                    >  2:11696/1‑151 (MQ=255)
                                               ttttGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTc                                                                                                    <  2:969551/150‑1 (MQ=255)
                                                tttGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCgg                                                                                                  >  2:942664/1‑151 (MQ=255)
                                                 ttGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCgg                                                                                                  <  2:968610/150‑1 (MQ=255)
                                                  tGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGGTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGc                                                                                                 >  2:217304/1‑150 (MQ=255)
                                                     gTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAg                                                                                             <  1:276101/151‑1 (MQ=255)
                                                      ttATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCaa                                                                                               >  1:1134443/1‑148 (MQ=255)
                                                      ttATACAGGAAGCTCGCCACTGTGAAGGAGGGACTGCTATGACGTCACTCTCTCGTCCCCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGGTTTTGGTTGTCTTCCTCGGCaa                                                                                               >  1:1263418/1‑148 (MQ=255)
                                                       tataCAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGa                                                                                            >  1:1361883/1‑150 (MQ=255)
                                                       tataCAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGa                                                                                            >  1:581465/1‑150 (MQ=255)
                                                        ataCAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGaa                                                                                           >  2:233279/1‑150 (MQ=255)
                                                           cAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACgg                                                                                       <  1:1279987/151‑1 (MQ=255)
                                                            aGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGt                                                                                      >  2:724362/1‑151 (MQ=255)
                                                            aGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGt                                                                                      >  1:79675/1‑151 (MQ=255)
                                                             ggAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTg                                                                                     <  1:1442884/151‑1 (MQ=255)
                                                                    cGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTgg                                                                              <  2:803604/151‑1 (MQ=255)
                                                                      ccACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCt                                                                            >  1:1516436/1‑151 (MQ=255)
                                                                      ccACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCt                                                                            >  2:1220879/1‑151 (MQ=255)
                                                                         cTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGAt                                                                         >  2:1562882/1‑151 (MQ=255)
                                                                           gtgAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGAtgt                                                                       <  1:1193485/151‑1 (MQ=255)
                                                                              aaGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCt                                                                    >  2:307944/1‑151 (MQ=255)
                                                                              aaGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCt                                                                    >  1:538292/1‑151 (MQ=255)
                                                                                 gaggTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGtt                                                                 <  1:28916/151‑1 (MQ=255)
                                                                                  aggTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTc                                                                <  1:285501/151‑1 (MQ=255)
                                                                                    gTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTcgc                                                              >  2:1221683/1‑151 (MQ=255)
                                                                                       cTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGcc                                                           >  1:1444031/1‑151 (MQ=255)
                                                                                         gCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGcc                                                           >  1:1135207/1‑149 (MQ=255)
                                                                                            aTGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGa                                                        <  1:1196780/149‑1 (MQ=255)
                                                                                            aTGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAAc                                                      <  1:1068182/151‑1 (MQ=255)
                                                                                            aTGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAAc                                                      >  2:1547582/1‑151 (MQ=255)
                                                                                              gACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACaa                                                    <  2:1686749/151‑1 (MQ=255)
                                                                                               aCGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACaa                                                    >  1:148160/1‑150 (MQ=255)
                                                                                               aCGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACa                                                     <  2:242783/149‑1 (MQ=255)
                                                                                                 gTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAAcc                                                 <  2:770007/151‑1 (MQ=255)
                                                                                                   cACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCa                                                <  2:105477/150‑1 (MQ=255)
                                                                                                   cACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCa                                                <  2:1694411/150‑1 (MQ=255)
                                                                                                    aCTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGc                                              >  2:1038940/1‑151 (MQ=255)
                                                                                                      tctctcGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCa                                             >  1:540341/1‑150 (MQ=255)
                                                                                                          tcGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAatat                                        >  1:148700/1‑151 (MQ=255)
                                                                                                            gTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATg                                       >  1:333898/1‑150 (MQ=255)
                                                                                                              ccGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGc                                    >  1:1510430/1‑151 (MQ=255)
                                                                                                               cgcgcgTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCt                                   <  1:918841/151‑1 (MQ=255)
                                                                                                                gcgcgTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTg                                  >  2:1460051/1‑151 (MQ=255)
                                                                                                                  gcgTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGt                                >  2:1473041/1‑151 (MQ=255)
                                                                                                                   cgTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGt                                >  2:764847/1‑150 (MQ=255)
                                                                                                                   cgTGGAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGt                                >  1:128131/1‑150 (MQ=255)
                                                                                                                      ggAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGaa                            <  2:185235/151‑1 (MQ=255)
                                                                                                                       gAGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGa                             >  1:1446173/1‑149 (MQ=255)
                                                                                                                        aGTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAgg                          <  1:846200/151‑1 (MQ=255)
                                                                                                                         gTTTATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAgg                          >  2:1580010/1‑150 (MQ=255)
                                                                                                                            tATCTCCACCATTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACt                       <  1:992517/150‑1 (MQ=255)
                                                                                                                                 ccaccaTTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACtggtgg                  <  2:1146575/150‑1 (MQ=255)
                                                                                                                                 ccaccaTTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACtggtg                   <  1:667948/149‑1 (MQ=255)
                                                                                                                                   accaTTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGttt               <  1:618307/151‑1 (MQ=255)
                                                                                                                                    ccaTTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGttt               <  2:1219774/150‑1 (MQ=255)
                                                                                                                                     caTTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTAc             <  2:306401/151‑1 (MQ=255)
                                                                                                                                     caTTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTAc             >  2:278785/1‑151 (MQ=255)
                                                                                                                                      aTTTTGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACt            >  2:910638/1‑151 (MQ=255)
                                                                                                                                       ttttGCAGACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACtt           >  2:332558/1‑151 (MQ=255)
                                                                                                                                          tGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTtc         <  2:1616489/150‑1 (MQ=255)
                                                                                                                                              gACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTAt    <  1:942591/151‑1 (MQ=255)
                                                                                                                                              gACCGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTAt    >  1:732161/1‑151 (MQ=255)
                                                                                                                                                ccGCACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTAttt  >  1:570190/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCCAGGTTGAACAAAACATACACAAAAAATATAGATCTCCGTCACATTTTTGCGTTATACAGGAAGCTCGCCACTGTGAAGGAGGTACTGCTATGACGTCACTCTCTCGTCCGCGCGTGGAGTTTATCTCCACCATTTTGCAGACCGTACTCAATCTTGGCCTGCTGTGCCTCGGCCTGATTTTGGTTGTCTTCCTCGGCAAAGAAACGGTGCATCTGGCTGATGTGCTGTTCGCGCCAGAACAAACCAGCAAATATGAGCTGGTAGAAGGACTGGTGGTTTACTTTCTCTATTT  >  NC_000913/4240233‑4240527

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