Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,617,382 A→C 6.2% E123A (GAA→GCA)  nikD → Ni(2(+)) ABC transporter ATP binding subunit NikD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,617,3820AC 6.2% 244.7 / 7.8 81E123A (GAA→GCA) nikDNi(2(+)) ABC transporter ATP binding subunit NikD
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base C (3/2);  total (41/40)
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.50e-01

TTGCGCCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACCACCGACCTTG  >  NC_000913/3617233‑3617528
                                                                                                                                                     |                                                                                                                                                  
ttGCGCCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGa                                                                                                                                                    <  1:777227/150‑1 (MQ=255)
   cgcCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGc                                                                                                                                                 >  1:961960/1‑150 (MQ=255)
     cccTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCgg                                                                                                                                               <  2:897521/150‑1 (MQ=255)
     cccTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCgg                                                                                                                                               >  1:1152474/1‑150 (MQ=255)
      ccTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTTCTATGGCAGCGGt                                                                                                                                              <  2:756101/150‑1 (MQ=255)
       cTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTTCTATGGCAGCGGta                                                                                                                                             <  1:1096934/150‑2 (MQ=255)
       cTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTg                                                                                                                                             <  1:793432/150‑1 (MQ=255)
        tGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTgg                                                                                                                                            <  1:1216087/150‑1 (MQ=255)
           gcgGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGc                                                                                                                                         >  2:979504/1‑150 (MQ=255)
           gcgGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGc                                                                                                                                         <  1:378278/150‑1 (MQ=255)
           gcgGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGc                                                                                                                                         <  2:715095/150‑1 (MQ=255)
                aTCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGaa                                                                                                                                    <  1:100509/150‑1 (MQ=255)
                  cAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGaaaa                                                                                                                                  >  2:242948/1‑150 (MQ=255)
                   aaaaTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAc                                                                                                                                 <  1:119399/150‑1 (MQ=255)
                   aaaaTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAc                                                                                                                                 >  1:1139813/1‑150 (MQ=255)
                    aaaTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAcg                                                                                                                                >  1:765946/1‑150 (MQ=255)
                    aaaTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAcg                                                                                                                                >  1:376702/1‑150 (MQ=255)
                     aaTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAcgc                                                                                                                               <  2:1305028/150‑1 (MQ=255)
                      aTTGCCACCATCATGCAGAACCAGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGAGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAAcgcc                                                                                                                              >  1:1400439/1‑150 (MQ=255)
                         gCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCcgcg                                                                                                                           >  2:118298/1‑150 (MQ=255)
                              catcatGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGt                                                                                                                      >  1:1194517/1‑150 (MQ=255)
                               atcatGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTg                                                                                                                     <  1:92565/150‑1 (MQ=255)
                                  atGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTg                                                                                                                  <  2:1305846/150‑1 (MQ=255)
                                    gCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                                                <  1:1214730/150‑1 (MQ=255)
                                    gCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                                                >  2:453121/1‑150 (MQ=255)
                                    gCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                                                >  2:295805/1‑150 (MQ=255)
                                      aGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGc                                                                                                              <  1:985103/150‑1 (MQ=255)
                                      aGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGc                                                                                                              <  1:304706/150‑1 (MQ=255)
                                      aGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGGTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGc                                                                                                              <  2:355882/150‑1 (MQ=255)
                                         aCCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGt                                                                                                           <  2:765946/150‑1 (MQ=255)
                                               gcAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGt                                                                                                     >  1:884973/1‑150 (MQ=255)
                                                     ccTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCgaga                                                                                               >  1:1026602/1‑150 (MQ=255)
                                                      cTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGAt                                                                                              >  1:1342867/1‑150 (MQ=255)
                                                      cTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGAt                                                                                              >  2:655428/1‑150 (MQ=255)
                                                       tttAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATg                                                                                             <  1:912622/150‑1 (MQ=255)
                                                         tAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAg                                                                                           >  1:1266343/1‑150 (MQ=255)
                                                           aTCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAgcg                                                                                         >  2:33395/1‑150 (MQ=255)
                                                               ccTGCACACCATGCACACCCACGCGCGGGAAACCTGCCGGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGcggcgg                                                                                     <  2:561904/149‑1 (MQ=255)
                                                                  gCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAg                                                                                                               <  2:1160297/121‑1 (MQ=255)
                                                                  gCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAg                                                                                                               >  1:1160297/1‑121 (MQ=255)
                                                                  gCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCg                                                                                                  <  2:1218232/134‑1 (MQ=255)
                                                                  gCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCg                                                                                                  >  1:1218232/1‑134 (MQ=255)
                                                                  gCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCAt                                                                                  >  1:440448/1‑150 (MQ=255)
                                                                     cacCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGtt                                                                               >  1:1187081/1‑150 (MQ=255)
                                                                           gCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAACCCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgcg                                                                         >  2:1202310/1‑150 (MQ=255)
                                                                             acacCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCa                                                                       <  1:893022/150‑1 (MQ=255)
                                                                             acacCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCa                                                                       <  2:1139813/150‑1 (MQ=255)
                                                                                        cgGGAAACCTGCCTGGCGTTAGTGAAACCCGCCGATGACGCTACGCTTACCGATGCCATACAAGCGGTGGGGCTGGAAAACCCCGCGCGCGTGCTGAACCTGTACCCGTTCGACATCAGCGCCGGCATGTTGCAGCGCATCATGATTgag                                                            >  2:204127/1‑148 (MQ=255)
                                                                                              acatgcctgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATAGCg                                                            >  1:973151/4‑144 (MQ=255)
                                                                                              aCCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCg                                                            <  2:973151/144‑1 (MQ=255)
                                                                                                ctgcctgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGt                                                    <  2:961960/150‑1 (MQ=255)
                                                                                                  gcctgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGc                                                  >  2:267444/1‑150 (MQ=255)
                                                                                                   cctgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCt                                                 <  1:1276331/150‑1 (MQ=255)
                                                                                                     tgGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtgt                                               >  2:1088134/1‑150 (MQ=255)
                                                                                                      ggCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtgtg                                              >  1:884597/1‑150 (MQ=255)
                                                                                                        cGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtgtgtg                                            <  2:1194517/150‑1 (MQ=255)
                                                                                                        cGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtgtgtg                                            >  2:1265348/1‑150 (MQ=255)
                                                                                                        cGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtgtgtg                                            <  2:1152474/150‑1 (MQ=255)
                                                                                                         gTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGa                                           <  1:979504/150‑1 (MQ=255)
                                                                                                          ttAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGaa                                          <  1:295805/150‑1 (MQ=255)
                                                                                                          ttAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCACGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGaa                                          <  1:242948/150‑1 (MQ=255)
                                                                                                           tAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAAt                                         >  2:612461/1‑150 (MQ=255)
                                                                                                            aGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACACGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATc                                        >  1:193066/1‑150 (MQ=255)
                                                                                                             gggAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCa                                       >  1:781928/1‑150 (MQ=255)
                                                                                                                aaaCCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCg                                    <  1:172830/150‑1 (MQ=255)
                                                                                                                  aCCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGtt                                  <  2:884973/150‑1 (MQ=255)
                                                                                                                   cccGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGttt                                 >  1:837715/1‑150 (MQ=255)
                                                                                                                   cccGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGttt                                 >  1:944023/1‑150 (MQ=255)
                                                                                                                   cccGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGGGCTGTGTGAATCACCGttt                                 <  2:945344/150‑1 (MQ=255)
                                                                                                                         gATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTatcatc                           >  2:275262/1‑150 (MQ=255)
                                                                                                                           tGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGc                         >  2:4361/1‑150 (MQ=255)
                                                                                                                             acgctacgctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCg                       <  2:1266343/150‑1 (MQ=255)
                                                                                                                               gctacgctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGAt                     <  2:1342867/150‑1 (MQ=255)
                                                                                                                               gctacgctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGAt                     <  1:1142179/150‑1 (MQ=255)
                                                                                                                                  agccccaCCGCTTCTATGGCAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGaa                  >  1:871740/7‑150 (MQ=255)
                                                                                                                                    gctTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAacc                <  2:1026602/150‑1 (MQ=255)
                                                                                                                                    ccccaCCGCTTCTATGGCAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAacc                >  1:850077/5‑150 (MQ=255)
                                                                                                                                     cccaCCGCTTCTATGGCAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccg               >  1:1410900/4‑150 (MQ=255)
                                                                                                                                       tACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccgac             >  1:267325/1‑150 (MQ=255)
                                                                                                                                        aCCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAaccgacc            <  1:118298/150‑1 (MQ=255)
                                                                                                                                                  tAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACCACCGACCTTg  <  1:275262/150‑1 (MQ=255)
                                                                                                                                                  tAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTAGAGATGAGCGGCGGCATCTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGACGATGAACCGACCACCGACCTTg  >  1:460040/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTGCGCCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACCACCGACCTTG  >  NC_000913/3617233‑3617528

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