Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,135,357 G→A 11.2% Q339* (CAG→TAG) ‡ pitB ← metal phosphate:H(+) symporter PitB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,135,3570GA11.2% 274.3 / 19.3 98Q339* (CAG→TAG) ‡pitBmetal phosphate:H(+) symporter PitB
Reads supporting (aligned to +/- strand):  ref base G (30/57);  new base A (5/6);  total (35/63)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 106 total aligned reads displayed.

TGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGG  >  NC_000913/3135210‑3135504
                                                                                                                                                   |                                                                                                                                                   
tgatgaTCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                                                                   <  1:1100843/150‑1 (MQ=255)
tgatgaTCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                                                                   >  1:201199/1‑150 (MQ=255)
    gatCCACACCGGAGCGTACTCAATGGTGCTTAACCTATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTgcg                                                                                                                                               <  2:1257956/150‑1 (MQ=255)
      tCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgcc                                                                                                                                             >  2:1449252/1‑149 (MQ=255)
      tCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgcc                                                                                                                                             <  1:72271/150‑2 (MQ=255)
      tCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgcc                                                                                                                                             <  1:1449252/150‑2 (MQ=255)
       ccACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTg                                                                                                                                            <  2:145986/150‑1 (MQ=255)
       ccACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccg                                                                                                                                            <  1:757518/150‑3 (MQ=255)
        cacaCCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGa                                                                                                                                           >  2:962872/1‑150 (MQ=255)
        cacaCCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCAGCGCAGCCAGCTgcgccgc                                                                                                                                           <  1:256166/150‑4 (MQ=255)
         acacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAc                                                                                                                                          <  1:949100/150‑1 (MQ=255)
         acacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAc                                                                                                                                          <  2:139786/150‑1 (MQ=255)
         acacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAc                                                                                                                                          <  2:1116464/150‑1 (MQ=255)
         acacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAc                                                                                                                                          <  1:826622/150‑1 (MQ=255)
         acacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAc                                                                                                                                          <  1:353412/150‑1 (MQ=255)
          cacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACt                                                                                                                                         <  1:1500108/150‑1 (MQ=255)
          cacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccgcat                                                                                                                                         <  1:80933/150‑6 (MQ=255)
            ccGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCa                                                                                                                                       >  2:949020/1‑150 (MQ=255)
              ggAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACg                                                                                                                                     <  2:1121375/150‑1 (MQ=255)
                     aCTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACg                                                                                                                              <  1:1144740/150‑1 (MQ=255)
                        cAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCt                                                                                                                           >  2:1139205/1‑150 (MQ=255)
                         aaTGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTc                                                                                                                          <  1:1526066/150‑1 (MQ=255)
                           tGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGt                                                                                                                        >  2:257478/1‑150 (MQ=255)
                              tGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                                                                     <  1:1325234/150‑1 (MQ=255)
                               gCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACt                                                                                                                    >  1:1291258/1‑150 (MQ=255)
                                 ttAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCAGACTCACGCTTAACCGCTCGTAACttt                                                                                                                  >  1:1441498/1‑150 (MQ=255)
                                       tatCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCatat                                                                                                            <  2:864998/150‑1 (MQ=255)
                                        atCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATAtt                                                                                                           >  1:601970/1‑150 (MQ=255)
                                           gCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGcc                                                                                                        <  1:890541/150‑1 (MQ=255)
                                           gCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGcc                                                                                                        <  2:1256969/150‑1 (MQ=255)
                                            cTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCt                                                                                                       >  2:560202/1‑150 (MQ=255)
                                             tGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTg                                                                                                      <  2:633196/150‑1 (MQ=255)
                                              gagAAGTTTTTTCAGCAGGTTCTGGTATTATTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACCCTTAACGGCTCATAACTTTCCATATTTCCTgg                                                                                                     >  2:653896/3‑150 (MQ=255)
                                              gCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTgg                                                                                                     >  2:1362797/1‑150 (MQ=255)
                                              gCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTgg                                                                                                     >  1:804366/1‑150 (MQ=255)
                                               cGAAGTTTTTTCAGAAGGTACTGGTCTTCTTTACTGACGCCTGGCAGTTTCCCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCCGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTgcc                                                                                                    >  2:617867/1‑148 (MQ=255)
                                                 aaGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgcag                                                                                                  <  1:492629/150‑1 (MQ=255)
                                                         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGttt                                                                                          <  1:889164/150‑1 (MQ=255)
                                                         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGttt                                                                                          <  1:1461347/150‑1 (MQ=255)
                                                          cagcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGtttt                                                                                         <  1:140160/150‑1 (MQ=255)
                                                          cagcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGtttt                                                                                         <  2:1201284/150‑1 (MQ=255)
                                                          cagcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGtttt                                                                                         <  1:302401/150‑1 (MQ=255)
                                                             cagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAc                                                                                      >  1:1351715/1‑150 (MQ=255)
                                                                gTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgcg                                                                                   <  1:506915/150‑1 (MQ=255)
                                                                 ttCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgcgc                                                                                  <  1:848925/150‑1 (MQ=255)
                                                                  tCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgcgcg                                                                                 <  1:279663/150‑1 (MQ=255)
                                                                   cTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCGCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACgcgcgc                                                                                >  2:378900/1‑150 (MQ=255)
                                                                    tGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCat                                                                                                              <  1:97533/119‑1 (MQ=255)
                                                                     ggTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCAGCATCGTTTTAACGCGCGCaa                                                                              >  2:725863/1‑150 (MQ=255)
                                                                        cttcttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAg                                                                           <  2:141559/150‑1 (MQ=255)
                                                                          tcttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCa                                                                         <  1:1135063/150‑1 (MQ=255)
                                                                          tcttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCa                                                                         <  1:761640/150‑1 (MQ=255)
                                                                          tcttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCa                                                                         <  1:344729/150‑1 (MQ=255)
                                                                           cttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCAt                                                                        <  1:560202/150‑1 (MQ=255)
                                                                           cttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATc                                                                       >  2:72500/1‑150 (MQ=255)
                                                                               aCTGCCGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCaaa                                                                    >  2:1160568/1‑150 (MQ=255)
                                                                                cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAg                                                                   >  2:1492766/1‑150 (MQ=255)
                                                                                  gACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGt                                                                 >  1:980253/1‑150 (MQ=255)
                                                                                    cGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTAt                                                               >  2:65588/1‑150 (MQ=255)
                                                                                           cAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCgg                                                        <  2:1099091/150‑1 (MQ=255)
                                                                                            aGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGa                                                       <  2:662291/150‑1 (MQ=255)
                                                                                                 cGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACa                                                  >  2:1158760/1‑150 (MQ=255)
                                                                                                  gCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAg                                                 >  1:831418/1‑150 (MQ=255)
                                                                                                  gCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAg                                                 >  1:4659/1‑150 (MQ=255)
                                                                                                  gCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAg                                                 <  2:1216278/150‑1 (MQ=255)
                                                                                                  gCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAg                                                 <  2:1291258/150‑1 (MQ=255)
                                                                                                    tAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTg                                               <  2:86551/150‑1 (MQ=255)
                                                                                                     aGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGa                                              <  2:201199/150‑1 (MQ=255)
                                                                                                     aGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGa                                              >  1:161565/1‑150 (MQ=255)
                                                                                                          cgcgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAActc                                         <  1:962872/150‑1 (MQ=255)
                                                                                                          cgcgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAActc                                         <  2:601970/150‑1 (MQ=255)
                                                                                                            cgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTg                                       <  1:1492766/150‑1 (MQ=255)
                                                                                                              gAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACCTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGtt                                     <  1:904714/150‑1 (MQ=255)
                                                                                                                   aTCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTg                                <  1:65588/150‑1 (MQ=255)
                                                                                                                     cAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTACCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCg                              <  1:50911/150‑1 (MQ=255)
                                                                                                                       gagaTGCACATCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTg                            <  1:402563/150‑1 (MQ=255)
                                                                                                                       gagaTGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTg                            <  2:161565/150‑1 (MQ=255)
                                                                                                                         gaTGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTATCGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGGGAAACTCTGTTACTTGCGTGcc                          <  1:1350900/150‑1 (MQ=255)
                                                                                                                           tGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCAt                        <  1:1158760/150‑1 (MQ=255)
                                                                                                                            gCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCCAAGGTATTTGCCGGATGACAGTGCAACTCTGTTACTTGCGTGCCATc                       <  2:1165077/150‑1 (MQ=255)
                                                                                                                            gCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATc                       >  1:953854/1‑150 (MQ=255)
                                                                                                                            gCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATc                       <  1:1139205/150‑1 (MQ=255)
                                                                                                                            gCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATc                       <  2:1359610/150‑1 (MQ=255)
                                                                                                                               caGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGt                    >  1:825047/1‑150 (MQ=255)
                                                                                                                                aGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTc                   <  2:825047/150‑1 (MQ=255)
                                                                                                                                 gCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCg                  >  2:814723/1‑150 (MQ=255)
                                                                                                                                    tAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATc                       >  2:23945/1‑142 (MQ=255)
                                                                                                                                    tAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATc                       <  1:23945/142‑1 (MQ=255)
                                                                                                                                      agtcagcgcAGCCAGCTGCGCTGACTCACGCTTACCGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAtgct             >  2:888424/6‑150 (MQ=255)
                                                                                                                                       tGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAtgctg            >  1:890311/1‑150 (MQ=255)
                                                                                                                                       gtcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAtgctg            >  2:242076/5‑150 (MQ=255)
                                                                                                                                         cagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCa          >  1:1247703/3‑150 (MQ=255)
                                                                                                                                         cagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCa          >  1:782440/3‑150 (MQ=255)
                                                                                                                                            cgcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGc       <  1:451723/150‑1 (MQ=255)
                                                                                                                                                 ctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAAt                                                                             >  1:1457495/1‑75 (MQ=255)
                                                                                                                                                 ctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAAt                                                                             <  2:1457495/75‑1 (MQ=255)
                                                                                                                                                 ctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATgg  <  2:1351715/150‑1 (MQ=255)
                                                                                                                                                 ctggctgCGCTGACTCACGCTTAACGGCGCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATgg  <  2:176535/150‑1 (MQ=255)
                                                                                                                                                   ggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGt                    <  1:163502/130‑1 (MQ=255)
                                                                                                                                                   ggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGACTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGt                    >  2:163502/1‑130 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGG  >  NC_000913/3135210‑3135504

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: