Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,616,399 T→G 5.4% 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.4% 278.8 / 7.6 93V73G (GTG→GGG) nikCNi(2(+)) ABC transporter membrane subunit NikC
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base G (3/2);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

GATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGT  >  NC_000913/3616250‑3616548
                                                                                                                                                     |                                                                                                                                                     
gATTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCgt                                                                                                                                                       >  1:679309/1‑150 (MQ=255)
 aTTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCgtg                                                                                                                                                      <  2:935440/150‑1 (MQ=255)
 aTTGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCgtg                                                                                                                                                      >  1:1432179/1‑150 (MQ=255)
  ttGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCgtgt                                                                                                                                                     >  2:988415/1‑150 (MQ=255)
   tGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTc                                                                                                                                                    <  2:1262342/150‑1 (MQ=255)
   tGCGCTCACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTc                                                                                                                                                    >  2:1712996/1‑150 (MQ=255)
         cACCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgcc                                                                                                                                              <  1:1585655/150‑4 (MQ=255)
          aCCAGCCAGTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggt                                                                                                                                             <  1:755523/150‑1 (MQ=255)
           ccagccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggtt                                                                                                                                            <  1:460955/150‑1 (MQ=255)
           ccagccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACACGCGGGTCGCtgccat                                                                                                                                            >  2:1317350/1‑145 (MQ=255)
            cagccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttg                                                                                                                                           <  1:788166/150‑2 (MQ=255)
              gccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttcgg                                                                                                                                         >  2:10892/1‑150 (MQ=255)
              gccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttcgg                                                                                                                                         <  1:1947985/150‑1 (MQ=255)
               ccagTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCATGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCtcggttcggt                                                                                                                                        >  2:1840515/1‑150 (MQ=255)
                 agTGGTGGCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTaa                                                                                                                                      <  2:274888/150‑1 (MQ=255)
                    ggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATgg                                                                                                                                   >  2:353734/1‑150 (MQ=255)
                    ggtggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCGGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATgg                                                                                                                                   <  2:1783271/150‑1 (MQ=255)
                      tggCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcc                                                                                                                                 <  2:359820/150‑1 (MQ=255)
                        gCTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctg                                                                                                                               <  2:1425064/150‑1 (MQ=255)
                         cTGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctgc                                                                                                                              >  1:126833/1‑150 (MQ=255)
                          tGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctgcc                                                                                                                             <  1:391160/150‑1 (MQ=255)
                          tGCCGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGgcctgcc                                                                                                                             <  1:1963993/150‑1 (MQ=255)
                            ccGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGcctgcctg                                                                                                                           >  1:387152/1‑150 (MQ=255)
                             cGTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGTTCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCctgcctgc                                                                                                                          <  1:1688916/150‑1 (MQ=255)
                              gTATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgct                                                                                                                         <  2:1863357/150‑1 (MQ=255)
                               tATGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctg                                                                                                                        >  1:1402230/1‑150 (MQ=255)
                                 tGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgct                                                                                                                      <  1:1359176/150‑1 (MQ=255)
                                 tGACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgct                                                                                                                      <  1:1366048/150‑1 (MQ=255)
                                  gACCCACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCctgctgctg                                                                                                                     >  1:1106939/1‑150 (MQ=255)
                                    cccACAGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGt                                                                                                                   >  1:1547472/1‑150 (MQ=255)
                                       acaGGCGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCt                                                                                                                <  1:659189/150‑1 (MQ=255)
                                         aGGCGATTTATTTTCCGTCGCGCCTGCTTTCGCGGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGa                                                                                                              <  2:387152/150‑1 (MQ=255)
                                            cGATTGATTTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATAACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTAGGTTAGGTAATGGCCTGCATGCTGCTGGTGCTGACAt                                                                                                           >  2:870939/1‑150 (MQ=255)
                                                 gattTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGGCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAggg                                                                                                      <  1:1673903/150‑1 (MQ=255)
                                                  attTGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGc                                                                                                     >  2:968191/1‑150 (MQ=255)
                                                     tGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTgg                                                                                                  <  2:1246305/150‑1 (MQ=255)
                                                     tGCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTgg                                                                                                  >  1:1598911/1‑150 (MQ=255)
                                                      gCCGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGt                                                                                                 >  1:53674/1‑150 (MQ=255)
                                                       ccGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGtt                                                                                                <  1:540953/150‑1 (MQ=255)
                                                       ccGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGtt                                                                                                >  1:1856419/1‑150 (MQ=255)
                                                       ccGTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGtt                                                                                                >  2:1642823/1‑150 (MQ=255)
                                                         gTCGCGCCTGCTTTCGCCGGATGCGCAGCACTGTCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTAt                                                                                              <  1:870939/150‑1 (MQ=255)
                                                           cgcgCCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTg                                                                                            >  1:198489/1‑150 (MQ=255)
                                                             cgcCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggc                                                                                          <  1:871785/150‑1 (MQ=255)
                                                             cgcCTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggc                                                                                          <  1:353734/150‑1 (MQ=255)
                                                                cTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggcggc                                                                                       <  2:209851/150‑1 (MQ=255)
                                                                cTGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggcggc                                                                                       <  2:1332021/150‑1 (MQ=255)
                                                                 tGCTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCa                                                                                      >  2:941768/1‑150 (MQ=255)
                                                                   cTTTCGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAgc                                                                                    >  2:1470562/1‑150 (MQ=255)
                                                                       cGCCGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCg                                                                                >  1:1012254/1‑150 (MQ=255)
                                                                         ccGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCggg                                                                              <  1:831803/150‑1 (MQ=255)
                                                                         ccGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCggg                                                                              >  2:1103305/1‑150 (MQ=255)
                                                                          cGGATGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGt                                                                             <  2:355233/150‑1 (MQ=255)
                                                                          cGGATGCGCAGCACTCGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTGTTGGCGGCAGCGCCGGGt                                                                             >  2:954813/1‑150 (MQ=255)
                                                                             aTGCGCAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGa                                                                          <  2:1598911/150‑1 (MQ=255)
                                                                               gcgcAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGAtt                                                                        <  2:1012254/150‑1 (MQ=255)
                                                                               gcgcAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGAtt                                                                        <  2:642552/150‑1 (MQ=255)
                                                                                cgcAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTg                                                                       >  1:464357/1‑150 (MQ=255)
                                                                                cgcAGCACTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTg                                                                       >  2:907113/1‑150 (MQ=255)
                                                                                 gcagcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTgg                                                                      >  2:1069048/1‑150 (MQ=255)
                                                                                   agcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGGGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGTGCTGTTTATTGGCGGCAGCGCCGGGTTTATTggcg                                                                    <  1:526381/150‑1 (MQ=255)
                                                                                   agcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcg                                                                    <  1:941768/150‑1 (MQ=255)
                                                                                    gcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGTTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgg                                                                   <  2:1547472/150‑1 (MQ=255)
                                                                                    gcaCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgg                                                                   >  2:1204259/1‑150 (MQ=255)
                                                                                      aCTGGCTGGGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgggc                                                                 <  2:679309/150‑1 (MQ=255)
                                                                                       ctggctgcGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTggcgggcg                                                                >  2:1547457/1‑150 (MQ=255)
                                                                                         ggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCCGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGcgcg                                                              >  1:1248073/1‑150 (MQ=255)
                                                                                         ggctggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGcgcg                                                              >  2:1121876/1‑150 (MQ=255)
                                                                                            tggGCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTg                                                           <  2:888967/150‑1 (MQ=255)
                                                                                             gggCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGa                                                          >  1:1079512/1‑150 (MQ=255)
                                                                                               gCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggcg                                                                                         <  1:326569/117‑1 (MQ=255)
                                                                                               gCACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTggcg                                                                                         >  2:326569/1‑117 (MQ=255)
                                                                                                cACCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCa                                                       >  1:1843038/1‑150 (MQ=255)
                                                                                                 aCCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAg                                                      <  1:968191/150‑1 (MQ=255)
                                                                                                 aCCGATCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAg                                                      <  1:435367/150‑1 (MQ=255)
                                                                                                      tCACTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGccac                                                 <  1:907113/150‑1 (MQ=255)
                                                                                                         cTTAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCAt                                              >  2:1366047/1‑150 (MQ=255)
                                                                                                           tAGGTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATgc                                            >  2:133931/1‑150 (MQ=255)
                                                                                                             ggTCGCGATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATgcgc                                          >  2:1303870/1‑150 (MQ=255)
                                                                                                                  cgATATTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGc                                     >  1:1399050/1‑150 (MQ=255)
                                                                                                                      atTTTCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGat                                 >  1:531186/1‑150 (MQ=255)
                                                                                                                        tttCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGatat                               <  1:1470562/150‑1 (MQ=255)
                                                                                                                        tttCTCGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGatat                               >  1:64091/1‑150 (MQ=255)
                                                                                                                          tctcGCGGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGt                             >  1:102209/1‑150 (MQ=255)
                                                                                                                             cgcgGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTa                          <  2:1106939/150‑1 (MQ=255)
                                                                                                                              gcgGCTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTAt                         >  2:1759327/1‑150 (MQ=255)
                                                                                                                                  cTGATGGCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGAcc                     >  2:903628/1‑150 (MQ=255)
                                                                                                                                    gaaccgaGCGACACGCGGGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCtt                   >  2:1520785/7‑150 (MQ=255)
                                                                                                                                      accgaGCGACACGCGGGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTcc                 >  1:1534896/5‑150 (MQ=255)
                                                                                                                                        gCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCg               >  2:49889/1‑150 (MQ=255)
                                                                                                                                        gCAGCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCg               <  1:1759327/150‑1 (MQ=255)
                                                                                                                                           gCGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGAcc            <  2:53674/150‑1 (MQ=255)
                                                                                                                                            cGACCCGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCt           <  2:126833/150‑1 (MQ=255)
                                                                                                                                               cccGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGa        >  1:879626/1‑150 (MQ=255)
                                                                                                                                               cccGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGa        >  2:1453190/1‑150 (MQ=255)
                                                                                                                                                ccGCGTGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGAt       <  2:1432179/150‑1 (MQ=255)
                                                                                                                                                acgcgGGTCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTATTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGAt       <  2:1534896/149‑1 (MQ=255)
                                                                                                                                                     tgtCGCTCGGTTCGGTAATGGCCTGCCTGCTGCTGGTGCTGACATTAGGGCTGGTTCTTGGCGGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGt  <  1:954813/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: