Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,617,382 A→C 5.8% E123A (GAA→GCA)  nikD → nickel ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,617,3820AC 5.8% 294.4 / 6.3 112E123A (GAA→GCA) nikDnickel ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (58/47);  new base C (3/4);  total (61/51)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 112 total aligned reads displayed.

CCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACCACCGACCTTGA  >  NC_000913/3617238‑3617529
                                                                                                                                                |                                                                                                                                                   
cccTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAACCCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCgg                                                                                                                                                >  2:661746/1‑150 (MQ=255)
 ccTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGt                                                                                                                                               <  2:1578167/150‑1 (MQ=255)
 ccTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGt                                                                                                                                               <  1:431330/150‑1 (MQ=255)
   tGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTgg                                                                                                                                             <  2:1212443/150‑1 (MQ=255)
    gcgcGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTggg                                                                                                                                            <  2:1051460/150‑1 (MQ=255)
    gcgcGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTggg                                                                                                                                            <  2:1491595/150‑1 (MQ=255)
    gcgcGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTggg                                                                                                                                            >  2:18210/1‑150 (MQ=255)
    gcgcGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCAGCCGATGACTCTACGCTTACCGCTTCTATGGCAGCGGtaag                                                                                                                                            <  1:989909/150‑4 (MQ=255)
     cgcgGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCTGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTgggg                                                                                                                                           <  2:285407/150‑1 (MQ=255)
      gcgGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTTCTATGGCAGCGGtaagcg                                                                                                                                          <  2:50907/150‑6 (MQ=255)
         gCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAACCCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTgg                                                                                                                                       >  2:1545517/1‑150 (MQ=255)
             cAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGaaaa                                                                                                                                   >  2:625836/1‑150 (MQ=255)
              aaaaTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAc                                                                                                                                  <  1:1032445/150‑1 (MQ=255)
               aaaTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAcg                                                                                                                                 <  1:1247830/150‑1 (MQ=255)
                 aTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAcgcc                                                                                                                               >  2:645777/1‑150 (MQ=255)
                   tGCCACCATCATGAAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAcgccgc                                                                                                                             <  2:217276/150‑1 (MQ=255)
                     ccaccaTCATGCAGAACCAGCGCAGAGCCATTAATCCACTGCACACCATGCACACCCACGCGCGCGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTAAGCTTACCTCTGCCATAGAAGCGGTCGGGCTGGAAAACGAagcgc                                                                                                                           >  1:723471/1‑150 (MQ=255)
                           tcatGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGc                                                                                                                     >  2:1672419/1‑150 (MQ=255)
                               gCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                                                 <  1:18210/150‑1 (MQ=255)
                               gCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                                                 >  1:1139326/1‑150 (MQ=255)
                                cAGAACCCGCGCAGCGCCTTTCATCCACTGCCCACCATGCACACCCCCGCGCGGGAAACCTGCCTGGCGTTAGGGAACCCCGCCGATGACGCTACGC‑TACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGc                                                                                                               <  1:1577549/150‑1 (MQ=255)
                                cAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCCTTAGGGAAACCCGCCGCTGACCCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAg                                                                                                                >  2:298884/1‑150 (MQ=255)
                                 aGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGc                                                                                                               >  2:690299/1‑150 (MQ=255)
                                 aGAACCCGCGCAGCGACTTTAATCCACTGCACACCATGCCAACCCACGCGCGCGAAACCTGACTGGAGATAGGGAAACCAGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGAGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGc                                                                                                               >  1:1439598/1‑150 (MQ=255)
                                    aCCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGt                                                                                                            >  2:595305/1‑150 (MQ=255)
                                     cccGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTa                                                                                                           >  2:1469888/1‑150 (MQ=255)
                                     ccagcgcAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTa                                                                                                           >  1:694510/4‑150 (MQ=255)
                                       cgcgCAGCGCCTTTACTCCACTGCACACCATGCTCACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTAcc                                                                                                         <  1:330412/150‑1 (MQ=255)
                                         cgcAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCg                                                                                                       >  2:54179/1‑150 (MQ=255)
                                         cgcAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGCGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCg                                                                                                       <  1:820731/150‑1 (MQ=255)
                                           cAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACCCTACGCTTACCGCTGCAATAGAAGCGGTGGGGCTGGAAAACGCCGCTCGCGTGCTGAAGCTGTACCCGtt                                                                                                     >  1:1075028/1‑150 (MQ=255)
                                             gcgcCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGTTGTACCCGTTCg                                                                                                   >  2:691493/1‑150 (MQ=255)
                                                 cTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGAt                                                                                               <  2:1648939/150‑1 (MQ=255)
                                                   ttAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCATGGCGGTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCTAGATGa                                                                                             <  2:124991/150‑1 (MQ=255)
                                                      aTCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAgcg                                                                                          >  1:1093245/1‑150 (MQ=255)
                                                        ccACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAgcggc                                                                                        <  2:165702/150‑1 (MQ=255)
                                                         cACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGCAGCTGTACCCGTTCGAGATGAgcggcg                                                                                       >  1:787424/1‑150 (MQ=255)
                                                         cACTGCACAACATCCACACCCACGCGCGGCAAACCTGCCTGCCGTTAGGGAAACCCGCCGATCACGCTACGCATACCGCTCCCATAGAAGCGGAGGGGCTCGAAAACCCCGCGCGCCTGCTGAAGCTGTACCCGTTCGAGATCAccggcg                                                                                       >  1:780852/1‑150 (MQ=255)
                                                          aCTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGACACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGcggcgg                                                                                      >  1:1662236/1‑150 (MQ=255)
                                                          aCTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGcggcgg                                                                                      >  2:1418676/1‑150 (MQ=255)
                                                          aCTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGcggcgg                                                                                      <  2:740273/150‑1 (MQ=255)
                                                            tGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGACGCTGTACCCGTTCGAGATGAGCGGCGGCa                                                                                    >  2:1706830/1‑150 (MQ=255)
                                                              ctcacCATGCACACCCACGCGCGGGAATCCTGCCTGGCGTTAGGGAAACCCGCCGATGTCGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATg                                                                                  <  1:1057690/148‑1 (MQ=255)
                                                                cacCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGtt                                                                                >  1:104267/1‑150 (MQ=255)
                                                                cacCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGtt                                                                                >  1:358502/1‑150 (MQ=255)
                                                                   cATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCa                                                                             <  2:304961/150‑1 (MQ=255)
                                                                    aTGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAg                                                                            <  1:54179/150‑1 (MQ=255)
                                                                    aTGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAg                                                                            >  1:501101/1‑150 (MQ=255)
                                                                       cacaCCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAACCCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgcgc                                                                         >  2:1069052/1‑150 (MQ=255)
                                                                         cacCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCat                                                                       <  2:373607/150‑1 (MQ=255)
                                                                         cacCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCat                                                                       <  2:17814/150‑1 (MQ=255)
                                                                          acCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCatg                                                                      >  1:878861/1‑150 (MQ=255)
                                                                           ccccCGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGTCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCatga                                                                     <  1:797263/150‑1 (MQ=255)
                                                                           cccACGCTCGGGAAACCTGCCTGGCGTTAGGGAAAGCCGCCGATGACGCTACGCTTACCGCTGCCATAGCAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCatga                                                                     <  1:321996/150‑1 (MQ=255)
                                                                              aCGCGCGGGAAACCGGCCTTGCGGTAGGTCAACCCGCCGATGACGCTACGCTGACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCAtgatga                                                                  <  1:96287/150‑1 (MQ=255)
                                                                               cgcgcgGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATgatgat                                                                 >  2:746577/1‑150 (MQ=255)
                                                                                  gcgGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGTTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGc                                                              <  1:691493/150‑1 (MQ=255)
                                                                                   cgGGAACCCTGCCAGGCGGTAGGGAAACCCGCCGCTGACGCTACGCTTCCCGCTGCCATAGCAGCGGTGGGGCTGGAAATCGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGCTGAGCGGCGGGATGTTGCAGCGCATGATGATTGCg                                                             <  1:152875/150‑1 (MQ=255)
                                                                                   cgGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCg                                                             >  1:355888/1‑150 (MQ=255)
                                                                                     gtaaaCCTGCCTGTCGTTAGGGAAACACGCCGATGACGCTACGCTTAACTCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGAt                                                           <  2:1649911/148‑1 (MQ=255)
                                                                                     ggAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGAt                                                           >  2:890284/1‑150 (MQ=255)
                                                                                     ggAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGGTCGAGATGAGCGGCGGCATGTTGCAGCGCATGCTGATTGCGAt                                                           >  2:326579/1‑150 (MQ=255)
                                                                                        aaCCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGc                                                        >  1:1442415/1‑150 (MQ=255)
                                                                                         aCCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCg                                                       <  2:381745/150‑1 (MQ=255)
                                                                                         aCCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCg                                                       >  2:1191669/1‑150 (MQ=255)
                                                                                          cctgcctgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCgg                                                      >  1:842694/1‑150 (MQ=255)
                                                                                            tgcctgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTg                                                    >  1:1426384/1‑150 (MQ=255)
                                                                                             gcctgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGc                                                   <  1:1469888/150‑1 (MQ=255)
                                                                                               ctgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtg                                                 >  2:877939/1‑150 (MQ=255)
                                                                                                gggCGGTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtgt                                                <  1:438854/149‑1 (MQ=255)
                                                                                                   cGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtgtgtg                                             <  2:1139326/150‑1 (MQ=255)
                                                                                                   cGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCagtgtg                                             >  1:1694534/1‑150 (MQ=255)
                                                                                                       aGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAAta                                         >  2:1655134/1‑149 (MQ=255)
                                                                                                        gggAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGTCTGGAAAACGCCGCGCGCGAGCTGAAGATGTACCCGTTCGACATGAGCGGCGGCAAGTTGCAGCGCATGATGATTGCGATGGAGGTCCTGTGTGAATCa                                        >  1:736965/1‑150 (MQ=255)
                                                                                                        gggAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCa                                        >  1:184076/1‑150 (MQ=255)
                                                                                                          gAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCAcc                                      >  2:1247971/1‑150 (MQ=255)
                                                                                                               ccgccgATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgcg                                                                          >  2:796279/1‑109 (MQ=255)
                                                                                                                cgccgATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTat                                <  1:1672419/150‑1 (MQ=255)
                                                                                                                cgccgATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTat                                <  2:1442415/150‑1 (MQ=255)
                                                                                                                cgccgATGACACTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTat                                <  2:851115/150‑1 (MQ=255)
                                                                                                                  ccgATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTatca                              <  1:890284/150‑1 (MQ=255)
                                                                                                                   cgATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTatcat                             >  2:1607664/1‑150 (MQ=255)
                                                                                                                       gACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGcc                         <  1:690299/150‑1 (MQ=255)
                                                                                                                         cgctacgctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGa                       <  1:645777/150‑1 (MQ=255)
                                                                                                                           ctacgctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATg                     <  2:22303/150‑1 (MQ=255)
                                                                                                                            cagccccaCCGCTTCTATGGCAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGa                    >  1:450004/8‑150 (MQ=255)
                                                                                                                             agccccaCCGCTTCTATGGCAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCGCCGTTTATCATCTCCGATGac                   >  1:877635/7‑149 (MQ=255)
                                                                                                                              gccccaCCGCTTCTATGGCAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAac                  >  1:170467/6‑150 (MQ=255)
                                                                                                                              cgctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAac                  <  1:1191669/150‑1 (MQ=255)
                                                                                                                               gctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAacc                 <  2:104267/150‑1 (MQ=255)
                                                                                                                                ctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccg                >  2:735240/1‑150 (MQ=255)
                                                                                                                                ctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccg                >  1:560099/1‑150 (MQ=255)
                                                                                                                                     cGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccgaccac           <  1:1530388/150‑1 (MQ=255)
                                                                                                                                      gCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCCGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGCaccgcccacc          >  2:1505274/1‑150 (MQ=255)
                                                                                                                                       cTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccgaccaccg         >  2:867137/1‑150 (MQ=255)
                                                                                                                                         gCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccgaccaccgac       >  1:544683/1‑150 (MQ=255)
                                                                                                                                          ccATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccgaccaccgacc      <  2:358502/150‑1 (MQ=255)
                                                                                                                                           cATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACCACCGACCt     >  1:816940/1‑150 (MQ=255)
                                                                                                                                             tAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACCACCGACCTTg   >  1:107384/1‑150 (MQ=255)
                                                                                                                                              aGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGATTATCATCGCCGATGAACCGACCACCGACCTTGa  >  1:1150575/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACCACCGACCTTGA  >  NC_000913/3617238‑3617529

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