Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,337,394 A→G S522G (AGC→GGC)  acnA → aconitate hydratase 1; aconitase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,337,3940AG96.7% 1287.9 / 15.7 363S522G (AGC→GGC) acnAaconitate hydratase 1; aconitase A
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base G (181/167);  total (187/176)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 380 total aligned reads displayed.

TTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTG  >  NC_000913/1337246‑1337541
                                                                                                                                                    |                                                                                                                                                   
ttCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGc                                                                                                                                                   >  2:1680669/1‑151 (MQ=255)
  cGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCg                                                                                                                                                  >  8:180821/1‑150 (MQ=255)
   gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAt                                                                                                                                                >  4:922091/1‑151 (MQ=255)
   gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAACGCGAt                                                                                                                                                >  5:689243/1‑151 (MQ=255)
   gAAAGTCGTTTCTGATTATATGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGAt                                                                                                                                                >  8:926344/1‑151 (MQ=255)
     aaGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAttt                                                                                                                                              <  7:1350092/151‑1 (MQ=255)
      aGTCGTTTCTGATTATCTGGCAAAAGCGAAACTTACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTa                                                                                                                                             <  4:1043931/151‑1 (MQ=255)
      aGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGAGACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTCTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTa                                                                                                                                             <  2:848434/151‑1 (MQ=255)
         cGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAcc                                                                                                                                          <  1:354551/151‑1 (MQ=255)
          gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCg                                                                                                                                         <  5:1386743/151‑1 (MQ=255)
            ttCTGATTATCTGGCATATGCGAAACTGACTGCGTATCTCGACGAACTGGGGTTTATCGTTGTGGGATACGGTTGTACCACGTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAAGGGCTATCAAATTAGGCGATTTAACCGt                                                                                                                                        <  6:493951/150‑1 (MQ=255)
            ttCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                       >  8:311137/1‑151 (MQ=255)
             tCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGCAACGGCAATCAAAAAAGGCCATTTAACCGTcg                                                                                                                                      >  2:1462680/1‑151 (MQ=255)
                gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcgg                                                                                                                                     <  5:207953/149‑1 (MQ=255)
                   tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcgg                                                                                                                                >  1:657950/1‑151 (MQ=255)
                   tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcg                                                                                                                                 <  8:1161357/150‑1 (MQ=255)
                   tATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcg                                                                                                                                 <  6:15587/150‑1 (MQ=255)
                    aTCTGGCCAACGCGCAACTGACACCGTAGCTCGCCGAACTGGGGTGTACCCTGGGGGGCTACGGTTGTACCCCCTGTATTGGTACCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAACGGCGATTTAACCGGcggtgcgg                                                                                                                                <  1:1004726/150‑1 (MQ=255)
                      cTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCggtgcggtgc                                                                                                                             >  2:313776/1‑151 (MQ=255)
                         gCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGt                                                                                                                          <  3:560560/151‑1 (MQ=255)
                               ccGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCa                                                                                                                    <  5:989512/150‑1 (MQ=255)
                                 gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGCTGCGGTGCTGTCCGGCAAc                                                                                                                  >  5:1666054/1‑151 (MQ=255)
                                 gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACAGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAc                                                                                                                  >  1:819566/1‑151 (MQ=255)
                                  aaaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAcc                                                                                                                 <  4:1233626/151‑1 (MQ=255)
                                        acacCGTATCTCGCCGAACTGGGGCTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACt                                                                                                           <  6:830221/151‑1 (MQ=255)
                                        acacCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACt                                                                                                           >  7:792588/1‑151 (MQ=255)
                                            cGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTTGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCCACCGTAACTTTg                                                                                                        >  6:166521/1‑150 (MQ=255)
                                            cGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGa                                                                                                       <  8:623782/151‑1 (MQ=255)
                                              tATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGa                                                                                                       >  8:458072/1‑149 (MQ=255)
                                                 ctcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACAGTAACTTTGAAggac                                                                                                  >  4:1569667/1‑149 (MQ=255)
                                                   cgacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGc                                                                                                   <  4:511545/148‑1 (MQ=255)
                                                   cgacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCg                                                                                                 <  5:1519967/150‑1 (MQ=255)
                                                    gacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTa                                                                                               <  7:1519105/151‑1 (MQ=255)
                                                      cgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                             >  6:1392338/1‑151 (MQ=255)
                                                       gaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatcc                                                                                            >  3:851212/1‑151 (MQ=255)
                                                        aaCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatcca                                                                                           <  4:1429849/151‑1 (MQ=255)
                                                        aaCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCACGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTACGGCAACCGTAACTTTGAAGGCCGTatc                                                                                             <  1:1062422/149‑1 (MQ=255)
                                                           tGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccatcc                                                                                        >  2:1530876/1‑151 (MQ=255)
                                                            ggggTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccatcc                                                                                        <  5:355043/150‑1 (MQ=255)
                                                               gTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTg                                                                                    <  5:260154/151‑1 (MQ=255)
                                                               gTTTAACCTTGTGGGATACGGTTGTACCACCTGTATGGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAATGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTg                                                                                    <  7:719092/151‑1 (MQ=255)
                                                                  tAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGtt                                                                                 >  3:773650/1‑151 (MQ=255)
                                                                   aaCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTa                                                                                <  7:583668/151‑1 (MQ=255)
                                                                     ccTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTaaa                                                                              <  2:597517/151‑1 (MQ=255)
                                                                     ccTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTaaa                                                                              <  3:1007556/151‑1 (MQ=255)
                                                                        tgtgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTaaaa                                                                             <  2:1135144/149‑1 (MQ=255)
                                                                          tgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaacta                                                                          >  8:821755/1‑150 (MQ=255)
                                                                           gggATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaactaac                                                                        <  3:236578/151‑1 (MQ=255)
                                                                            ggATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaactaact                                                                       >  1:478959/1‑151 (MQ=255)
                                                                             gATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActg                                                                      >  3:340802/1‑151 (MQ=255)
                                                                             gATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActg                                                                      <  3:669723/151‑1 (MQ=255)
                                                                              aTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActgg                                                                     <  7:991149/151‑1 (MQ=255)
                                                                                         ccaccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgcc                                                          <  7:677281/151‑1 (MQ=255)
                                                                                         acaccTGTATTTTTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgcc                                                          <  2:1314820/150‑1 (MQ=255)
                                                                                          caccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgccg                                                         <  8:399119/151‑1 (MQ=255)
                                                                                             cTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgccgccg                                                      >  8:1586236/1‑151 (MQ=255)
                                                                                                 aTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtg                                                   <  5:947869/150‑1 (MQ=255)
                                                                                                   tGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggt                                                 <  5:49584/150‑1 (MQ=255)
                                                                                                    ggTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCGACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtgg                                               >  2:354231/1‑151 (MQ=255)
                                                                                                    ggTAACTCTGGGCCCCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGCTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCAGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtcgtgg                                               >  6:1291642/1‑151 (MQ=255)
                                                                                                       aaCTCTGGGCCGGTGCCAGATCCTCTCGAAACGCCAATCAAAAAAGGCGATTTAACCTTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAATGCCGTATCCATCGGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGtt                                             <  4:1183642/150‑1 (MQ=255)
                                                                                                       aaCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTACCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTg                                            <  6:1198076/151‑1 (MQ=255)
                                                                                                             gggCCTCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCt                                         <  3:1416297/148‑1 (MQ=255)
                                                                                                             gggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCt                                         <  7:1435504/148‑1 (MQ=255)
                                                                                                             gggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTAt                                       <  5:810336/150‑1 (MQ=255)
                                                                                                               gCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgcg                                    >  4:673979/1‑151 (MQ=255)
                                                                                                               gCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgcg                                    <  6:686463/151‑1 (MQ=255)
                                                                                                                 cGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCt                                  <  7:163974/151‑1 (MQ=255)
                                                                                                                    tGCCCGATCCTATCGAAACGGCAATCAAATAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGc                               >  3:1313632/1‑151 (MQ=255)
                                                                                                                     gCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCGGCTGGTGGTTGCCTATGCGCTGGCg                              >  2:877663/1‑151 (MQ=255)
                                                                                                                      cccGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCgg                             >  6:1226908/1‑151 (MQ=255)
                                                                                                                      cccGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCgg                             >  8:1475205/1‑151 (MQ=255)
                                                                                                                      cccGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCgg                             >  1:1141588/1‑151 (MQ=255)
                                                                                                                       ccGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCggg                            >  5:1070271/1‑151 (MQ=255)
                                                                                                                        cGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGa                           <  7:927476/151‑1 (MQ=255)
                                                                                                                          aTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTGAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGaa                          >  8:1405559/1‑150 (MQ=255)
                                                                                                                           tCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGatat                        >  4:1446539/1‑151 (MQ=255)
                                                                                                                           tCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAat                        >  1:961814/1‑151 (MQ=255)
                                                                                                                           tCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAat                        >  7:1210308/1‑151 (MQ=255)
                                                                                                                            ccTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAata                       <  7:415776/151‑1 (MQ=255)
                                                                                                                            ccTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAat                        >  1:55496/1‑150 (MQ=255)
                                                                                                                                tCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGCTTGCCTATGCGCTGGCGGGAAATATGa                    >  1:1544555/1‑150 (MQ=255)
                                                                                                                                 cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGGTGGTGGTTGCCTATGCGCTGGGGGGAAATATGAat                  >  6:1453573/1‑151 (MQ=255)
                                                                                                                                 cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGaa                   <  3:168955/150‑1 (MQ=255)
                                                                                                                                 cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAat                  <  3:1117671/151‑1 (MQ=255)
                                                                                                                                  gAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAata                 <  6:1049392/151‑1 (MQ=255)
                                                                                                                                  gAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGACTCGCCGCCGCTGGTGGTTGCCTAGGCGCTGGCGGGAAATATGAat                  >  6:983671/1‑150 (MQ=255)
                                                                                                                                      cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCaa             <  1:1333333/151‑1 (MQ=255)
                                                                                                                                      cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCaa             <  5:849858/151‑1 (MQ=255)
                                                                                                                                         cAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCt          >  8:1051274/1‑151 (MQ=255)
                                                                                                                                         cAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCt          >  6:366195/1‑151 (MQ=255)
                                                                                                                                          aaTCAAAAAAGGCGATCTAACCGTCGGTGCGGTGCTGTCCGGCAACGGTAAGTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATTTCAACCtc         <  1:1208831/151‑2 (MQ=255)
                                                                                                                                           aTCAAAATAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGTCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTgg        <  2:530001/151‑1 (MQ=255)
                                                                                                                                           aTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCACCCTgg        >  3:1630113/1‑151 (MQ=255)
                                                                                                                                            tCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGc       >  3:1190837/1‑151 (MQ=255)
                                                                                                                                              aaaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCtt     >  1:171129/1‑151 (MQ=255)
                                                                                                                                              aaaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCtt     <  8:1211852/151‑1 (MQ=255)
                                                                                                                                              aaaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCtt     <  5:128638/151‑1 (MQ=255)
                                                                                                                                                aaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCACCCTGGCTTCt   >  4:406834/1‑151 (MQ=255)
                                                                                                                                                aaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCt   <  4:1673480/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTG  >  NC_000913/1337246‑1337541

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