Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 799,850 A→T 10.8% L209Q (CTG→CAG)  ispB ← octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947799,8500AT10.8% 188.0 / 13.8 65L209Q (CTG→CAG) ispBoctaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base T (5/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

CCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACTTCCTCG  >  NC_002947/799704‑799999
                                                                                                                                                  |                                                                                                                                                     
ccTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCt                                                                                                                                                    >  2:1045464/1‑150 (MQ=255)
    gCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGGTCGACGACCTGCTGGac                                                                                                                                                <  1:174740/150‑2 (MQ=255)
      gTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                              >  2:160013/1‑150 (MQ=255)
      gTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                              >  1:1146002/1‑150 (MQ=255)
          cAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCAc                                                                                                                                          <  1:1077609/150‑1 (MQ=255)
               ccgccTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCa                                                                                                                                     <  2:702484/150‑1 (MQ=255)
                 gccTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGa                                                                                                                                   >  2:169145/1‑150 (MQ=255)
                      cTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTc                                                                                                                              <  1:1211224/150‑1 (MQ=255)
                        cGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGc                                                                                                                            <  1:1164142/150‑1 (MQ=255)
                             gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGTAGTCCACGCCCAGATGGTCGCCAAAg                                                                                                                       <  1:814029/150‑1 (MQ=255)
                             gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                                       <  1:160013/150‑1 (MQ=255)
                                  ttCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACg                                                                                                                  >  1:109386/1‑150 (MQ=255)
                                  ttCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACg                                                                                                                  >  1:1149606/1‑150 (MQ=255)
                                     gcgcATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAg                                                                                                               <  2:109386/150‑1 (MQ=255)
                                          tGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCtt                                                                                                          >  1:756401/1‑150 (MQ=255)
                                               tAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTg                                                                                                                                                   >  2:643182/1‑104 (MQ=255)
                                               tAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTg                                                                                                                                                   <  1:643182/104‑1 (MQ=255)
                                                aGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCgct                                                                                                    >  2:1087365/1‑150 (MQ=255)
                                                 gATCAGCGGCAGTGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTGCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCgctg                                                                                                   >  2:1013480/1‑150 (MQ=255)
                                                 gATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCGCAGGGTTTGCGAGTCCCCCTGGTAGTCCAGCAGATCGACAACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGACTTCGCgctg                                                                                                   >  2:494250/1‑150 (MQ=255)
                                                  aTCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCgctgc                                                                                                  <  2:683117/150‑1 (MQ=255)
                                                  aTCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGCTCTGGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCgctgc                                                                                                  >  2:501244/1‑150 (MQ=255)
                                                   tCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCgctgct                                                                                                 >  2:401385/1‑150 (MQ=255)
                                                         gCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCgg                                                                                           <  1:1087365/150‑1 (MQ=255)
                                                          cAGGGGAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGt                                                                                          >  1:1025922/1‑150 (MQ=255)
                                                           aGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTg                                                                                         <  1:1058187/150‑1 (MQ=255)
                                                           aGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTg                                                                                         <  1:169145/150‑1 (MQ=255)
                                                               tAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCt                                                                                     >  1:1135184/1‑150 (MQ=255)
                                                                 ggTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTgg                                                                                   <  2:1193788/150‑1 (MQ=255)
                                                                  gTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTggg                                                                                  <  1:1013480/150‑1 (MQ=255)
                                                                        ccTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAgcg                                                                            <  2:1140216/150‑1 (MQ=255)
                                                                                      tcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                              >  1:473011/1‑70 (MQ=255)
                                                                                      tcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                              <  2:473011/70‑1 (MQ=255)
                                                                                          cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                                <  1:366539/144‑1 (MQ=255)
                                                                                          cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                                >  2:366539/1‑144 (MQ=255)
                                                                                            aCGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGaa                                                        <  1:1045464/150‑1 (MQ=255)
                                                                                              gTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGc                                                      >  2:314058/1‑150 (MQ=255)
                                                                                                tCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCtt                                                    >  2:542610/1‑150 (MQ=255)
                                                                                                      ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                                   <  2:900927/145‑1 (MQ=255)
                                                                                                      ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                                   >  1:900927/1‑145 (MQ=255)
                                                                                                      ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACa                                              >  1:155092/1‑150 (MQ=255)
                                                                                                      ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACa                                              <  1:314058/150‑1 (MQ=255)
                                                                                                      ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACa                                              <  2:155092/150‑1 (MQ=255)
                                                                                                        aGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGc                                            <  2:1149606/150‑1 (MQ=255)
                                                                                                        aGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGc                                            <  2:117586/150‑1 (MQ=255)
                                                                                                            tCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCg                                        >  1:683692/1‑150 (MQ=255)
                                                                                                              tGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGcc                                      >  1:823066/1‑150 (MQ=255)
                                                                                                                 gAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTc                                   <  2:1135184/150‑1 (MQ=255)
                                                                                                                  aGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCt                                  <  1:401385/150‑1 (MQ=255)
                                                                                                                         cTTGTAGTCCAGCAGGTCGTCGACCAGCGGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCcgcg                           >  1:236109/1‑150 (MQ=255)
                                                                                                                          tggccttCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGACCTGGGCCAGCGCCGCAGCACAGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGa                          >  1:219349/7‑150 (MQ=255)
                                                                                                                                 ccAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCt                   >  1:153763/1‑150 (MQ=255)
                                                                                                                                 ccAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCt                   >  1:746426/1‑150 (MQ=255)
                                                                                                                                  cagctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATTACCTc                  >  2:598855/1‑150 (MQ=255)
                                                                                                                                  cagcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTc                  <  2:756401/150‑1 (MQ=255)
                                                                                                                                    gctggTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                                   >  2:365044/4‑115 (MQ=255)
                                                                                                                                    gctggTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                                   <  1:365044/112‑1 (MQ=255)
                                                                                                                                           gtcgACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGAc         >  1:1063077/1‑150 (MQ=255)
                                                                                                                                           gtcgACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGAc         >  2:1016549/1‑150 (MQ=255)
                                                                                                                                           gtcgACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGAc         >  1:1063185/1‑150 (MQ=255)
                                                                                                                                            tcgACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACt        >  1:311436/1‑150 (MQ=255)
                                                                                                                                                 cAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATg              >  2:445684/1‑139 (MQ=255)
                                                                                                                                                 cAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATg              <  1:445684/139‑1 (MQ=255)
                                                                                                                                                 cAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACTTCCTc   >  1:1150094/1‑150 (MQ=255)
                                                                                                                                                 cAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACTTCCTc   <  2:823066/150‑1 (MQ=255)
                                                                                                                                                  aGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACTTCCTCg  >  1:933404/1‑150 (MQ=255)
                                                                                                                                                  aGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACTTCCTCg  <  2:311436/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACTTCCTCG  >  NC_002947/799704‑799999

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