Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,616,399 T→G 5.1% V73G (GTG→GGG)  nikC → Ni(2(+)) ABC transporter membrane subunit NikC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,616,3990TG 5.1% 241.0 / 4.7 79V73G (GTG→GGG) nikCNi(2(+)) ABC transporter membrane subunit NikC
Reads supporting (aligned to +/- strand):  ref base T (38/37);  new base G (2/2);  total (40/39)
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.74e-01

GATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGT  >  NC_000913/3616250‑3616548
                                                                                                                                                     |                                                                                                                                                     
gATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCgt                                                                                                                                                       >  1:696582/1‑150 (MQ=255)
  ttGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCgtgt                                                                                                                                                     >  2:1269407/1‑150 (MQ=255)
  ttGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCgtgt                                                                                                                                                     <  1:654218/150‑1 (MQ=255)
    gcgcTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCg                                                                                                                                                   >  1:978949/1‑150 (MQ=255)
      gcTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGcg                                                                                                                                                 >  1:662535/1‑149 (MQ=255)
       cgcACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGTTTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtc                                                                                                                                                <  2:1646284/148‑1 (MQ=255)
        tCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGAGGCGCAGCACTGGCTGGGCACCGATCACTTAGGGCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgc                                                                                                                                               <  1:215123/150‑3 (MQ=255)
          aCCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgcca                                                                                                                                             <  1:113987/150‑5 (MQ=255)
            cagccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttc                                                                                                                                           <  2:254222/150‑1 (MQ=255)
               ccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttcggt                                                                                                                                        <  1:1402794/150‑1 (MQ=255)
               ccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttcggt                                                                                                                                        <  2:900157/150‑1 (MQ=255)
                   tggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATg                                                                                                                                    >  1:1100656/1‑150 (MQ=255)
                      tggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcc                                                                                                                                 >  1:895465/1‑150 (MQ=255)
                          tGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctgcc                                                                                                                             <  1:424475/150‑1 (MQ=255)
                             cGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCctgcctgc                                                                                                                          <  2:943312/150‑1 (MQ=255)
                              gTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgct                                                                                                                         <  1:674984/150‑1 (MQ=255)
                                aTGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCCCATGCGCAGCACTGGCAGGGCACCGATAACTTAGGTCGCGATATTTTCTCGAGGCTGATCGCACCGACCCGCGTGTCGCTCGGTTAGCTAATGGCCTGCcagctcc                                                                                                                       >  2:1634436/1‑148 (MQ=255)
                                   aCCCACAGGCGATTGATTTGCCGGCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTgg                                                                                                                    <  2:407679/150‑1 (MQ=255)
                                    cccACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGt                                                                                                                   >  1:1089252/1‑150 (MQ=255)
                                      cacaGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGc                                                                                                                 <  1:77818/150‑1 (MQ=255)
                                        caGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTg                                                                                                               <  1:832873/150‑1 (MQ=255)
                                        caGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGGCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTg                                                                                                               <  1:452665/150‑1 (MQ=255)
                                            cGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACAt                                                                                                           <  2:1190598/150‑1 (MQ=255)
                                              attgattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTa                                                                                                         <  2:1007257/150‑1 (MQ=255)
                                               ttgattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACAtt                                                                                                          >  1:1446800/1‑148 (MQ=255)
                                               ttgattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACAtt                                                                                                          <  2:1446800/148‑1 (MQ=255)
                                                 gattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAggg                                                                                                      >  1:136863/1‑150 (MQ=255)
                                                 gattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAggg                                                                                                      <  1:493173/150‑1 (MQ=255)
                                                  attTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGc                                                                                                     >  2:764698/1‑150 (MQ=255)
                                                     tGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTgg                                                                                                  >  1:1579242/1‑150 (MQ=255)
                                                        cGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTa                                                                                               >  1:394536/1‑150 (MQ=255)
                                                        cGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTa                                                                                               >  1:400542/1‑150 (MQ=255)
                                                          tCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTAtt                                                                                             <  1:1193171/150‑1 (MQ=255)
                                                           cgcgCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTg                                                                                            >  2:53740/1‑150 (MQ=255)
                                                             cgcCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggc                                                                                          >  2:528284/1‑150 (MQ=255)
                                                                    tttCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAgcg                                                                                   >  2:1316760/1‑150 (MQ=255)
                                                                    tttCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAgcg                                                                                   >  1:1059677/1‑150 (MQ=255)
                                                                     ttCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAgcgc                                                                                  <  2:895465/150‑1 (MQ=255)
                                                                      tCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGcc                                                                                 >  1:1222284/1‑150 (MQ=255)
                                                                       cGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGGTGGTTATTGGCGGCAGCGCCg                                                                                >  1:1512032/1‑150 (MQ=255)
                                                                           ggATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGtt                                                                            >  1:143755/1‑150 (MQ=255)
                                                                                cgcAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTg                                                                       >  1:839304/1‑150 (MQ=255)
                                                                                cgcAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTg                                                                       >  1:1114856/1‑150 (MQ=255)
                                                                                cgcAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCAACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTg                                                                       >  1:988233/1‑150 (MQ=255)
                                                                                  cagcaCTGGCTGGGCAGCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggg                                                                     <  1:659685/150‑2 (MQ=255)
                                                                                  cagcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggc                                                                     >  2:782145/1‑150 (MQ=255)
                                                                                    gcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgg                                                                   <  1:252266/150‑1 (MQ=255)
                                                                                       ctggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgggcg                                                                <  2:479032/150‑1 (MQ=255)
                                                                                       ctggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgggcg                                                                >  2:1137500/1‑150 (MQ=255)
                                                                                        tggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGgcgc                                                               <  1:482244/150‑1 (MQ=255)
                                                                                              ggCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGAt                                                         >  1:904833/1‑150 (MQ=255)
                                                                                                cACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCa                                                       >  2:627686/1‑150 (MQ=255)
                                                                                                    gATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGcc                                                   >  2:316610/1‑150 (MQ=255)
                                                                                                       cACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGccacc                                                <  1:1559649/150‑1 (MQ=255)
                                                                                                       cACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGccacc                                                <  2:696582/150‑1 (MQ=255)
                                                                                                        aCTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAg                                                      <  1:1281971/143‑1 (MQ=255)
                                                                                                        aCTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAg                                                      >  2:1281971/1‑143 (MQ=255)
                                                                                                        aCTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGccacca                                               >  1:243063/1‑150 (MQ=255)
                                                                                                        aCTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGccacca                                               <  2:1619104/150‑1 (MQ=255)
                                                                                                         cTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCGTGCCTGATTCTGGTGCTGAAATTAGTGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGGGGGGGCGTTGATCAGGACACCAt                                              <  2:1157423/150‑1 (MQ=255)
                                                                                                         cTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGccaccaa                                              >  2:1061166/1‑149 (MQ=255)
                                                                                                         cTTAGGTCGAGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCAt                                              >  2:1085537/1‑150 (MQ=255)
                                                                                                          ttAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATg                                             <  2:1100656/150‑1 (MQ=255)
                                                                                                           tAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATgc                                            <  1:627686/150‑1 (MQ=255)
                                                                                                             ggTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATgcgc                                          >  2:443113/1‑150 (MQ=255)
                                                                                                                cgcgATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTc                                       >  1:1390794/1‑150 (MQ=255)
                                                                                                                 gcgATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCg                                      <  1:53740/150‑1 (MQ=255)
                                                                                                                  cgATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGc                                     <  1:1137500/150‑1 (MQ=255)
                                                                                                                         ttCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATg                              >  2:980929/1‑150 (MQ=255)
                                                                                                                          tctcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGt                             >  1:626984/1‑150 (MQ=255)
                                                                                                                           ctcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGtt                            >  2:1617576/1‑150 (MQ=255)
                                                                                                                             cgcgGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTa                          >  1:397908/1‑150 (MQ=255)
                                                                                                                              gcgGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTAt                         <  1:980929/150‑1 (MQ=255)
                                                                                                                               cgGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACAATGCGCGTCGCCGATATGTTTATg                        <  2:1222284/150‑1 (MQ=255)
                                                                                                                                ggCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGa                       <  2:52473/150‑1 (MQ=255)
                                                                                                                                      accgaGCGACACGCGGGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTcc                 >  1:620057/5‑150 (MQ=255)
                                                                                                                                       ggCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTccc                >  1:959579/1‑150 (MQ=255)
                                                                                                                                         gaGCGACACGCGGGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGa              >  2:795567/2‑150 (MQ=255)
                                                                                                                                           gCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGAcc            <  1:1477421/150‑1 (MQ=255)
                                                                                                                                           gCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGAcc            <  2:978949/150‑1 (MQ=255)
                                                                                                                                            cGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCt           <  1:528284/150‑1 (MQ=255)
                                                                                                                                            cGACCCGCGTGTCGCTCGGTGCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCt           <  1:797527/150‑1 (MQ=255)
                                                                                                                                                     tgtCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGt  <  2:143755/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGT  >  NC_000913/3616250‑3616548

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: ‑