Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 799,844 T→A 16.7% D211V (GAC→GTC)  ispB ← octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947799,8440TA16.7% 151.6 / 19.8 59D211V (GAC→GTC) ispBoctaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (23/26);  new base A (7/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACT  >  NC_002947/799704‑799993
                                                                                                                                            |                                                                                                                                                     
ccTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCt                                                                                                                                              >  2:1045464/1‑150 (MQ=255)
  tGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGGTCGACGAcc                                                                                                                                                  >  1:1016784/1‑144 (MQ=255)
  tGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGGTCGACGAcc                                                                                                                                                  <  2:1016784/144‑1 (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:169145/150‑1 (MQ=255)
                                                           aGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTg                                                                                   <  1:1058187/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                                                          >  2:366539/1‑144 (MQ=255)
                                                                                          cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                          <  1:366539/144‑1 (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:746426/1‑150 (MQ=255)
                                                                                                                                 ccAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCt             >  1:153763/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   >  2:1016549/1‑150 (MQ=255)
                                                                                                                                           gtcgACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGAc   >  1:1063185/1‑150 (MQ=255)
                                                                                                                                           gtcgACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGAc   >  1:1063077/1‑150 (MQ=255)
                                                                                                                                            tcgACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACt  >  1:311436/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCATGTAGACT  >  NC_002947/799704‑799993

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: