Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,617,375 G→T 6.1% A121S (GCC→TCC) ‡ 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,3750GT 6.1% 276.7 / 5.8 82A121S (GCC→TCC) ‡nikDNi(2(+)) ABC transporter ATP binding subunit NikD
Reads supporting (aligned to +/- strand):  ref base G (39/38);  new base T (3/2);  total (42/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 = 9.25e-01

CCTTGCGCCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACC  >  NC_000913/3617231‑3617518
                                                                                                                                                |                                                                                                                                               
ccTTGCGCCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATa                                                                                                                                            >  1:715095/1‑150 (MQ=255)
 cTTGCGCCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAg                                                                                                                                           >  1:1305846/1‑150 (MQ=255)
  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                                                                                                      >  2:453121/1‑150 (MQ=255)
                                      gCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                                      >  2:295805/1‑150 (MQ=255)
                                      gCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                                      <  1:1214730/150‑1 (MQ=255)
                                        aGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGc                                                                                                    <  1:304706/150‑1 (MQ=255)
                                        aGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGc                                                                                                    <  1:985103/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                                                                                    >  2:655428/1‑150 (MQ=255)
                                                        cTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGAt                                                                                    >  1:1342867/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                                                                                                     >  1:1160297/1‑121 (MQ=255)
                                                                    gCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAg                                                                                                     <  2:1160297/121‑1 (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:1152474/150‑1 (MQ=255)
                                                                                                          cGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCtgtgtg                                  >  2:1265348/1‑150 (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)
                                                                                                                                                |                                                                                                                                               
CCTTGCGCCCTGCGCGGCATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCTGGCGTTAGGGAAACCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCGATGAACCGACC  >  NC_000913/3617231‑3617518

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: