Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,647,612 G→A D384D (GAC→GAT oprB‑II ← carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,647,6120GA98.4% 191.8 / ‑4.8 64D384D (GAC→GAToprB‑IIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (30/33);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCCTGCACATAGTT  >  NC_002947/1647469‑1647758
                                                                                                                                               |                                                                                                                                                    
ccATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATag                                                                                                                                                <  2:801925/150‑1 (MQ=255)
      aCCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGGGTGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTGGTAATCATagcgcgcg                                                                                                                                          >  2:932273/1‑144 (MQ=255)
          gCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCcggcgg                                                                                                                                      >  1:875990/1‑150 (MQ=255)
            tacacCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGcc                                                                                                                                    >  1:71531/2‑150 (MQ=255)
                ccccGCCTGGGTGGCGGCAATACTGCAGGTTGGGGCGTCCAGAGCGCCAGTCGGCAAAGCGACTGCCGTAATACATTTATCCGGTTTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGt                                                                                                                                <  2:953216/150‑1 (MQ=255)
                 cccgcccgGGTGGCGGATGTACTTCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCACGTGAATGCCGTAATACATCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGGTGTCGTAATGATAGAGGCCGGCGGGCTCGtt                                                                                                                               >  2:267234/1‑150 (MQ=255)
                         ggTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTCAACCAggc                                                                                                                       <  2:1185918/150‑1 (MQ=255)
                           tGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAggcgg                                                                                                                     <  1:978736/150‑1 (MQ=255)
                            ggCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAggcggg                                                                                                                    <  2:708945/150‑1 (MQ=255)
                              cGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAggcgggcg                                                                                                                  >  1:512890/1‑150 (MQ=255)
                                  tGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCCGACTGGTTGACTAGGTGGGTGTTCt                                                                                                              <  2:1090603/150‑1 (MQ=255)
                                   gTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCtt                                                                                                             >  2:1028180/1‑150 (MQ=255)
                                     aCTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑c                                                                                                         <  2:875990/150‑1 (MQ=255)
                                      cTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCGTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑Cg                                                                                                        <  2:324307/150‑1 (MQ=255)
                                       tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGa                                                                                                       >  1:321761/1‑150 (MQ=255)
                                       tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGa                                                                                                       <  2:612237/150‑1 (MQ=255)
                                            gTTGGGGCGTACCGTGAGCCAGTCGACCAAGTGAATGCCGTAATACCGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGc                                                                                                  >  1:880257/1‑150 (MQ=255)
                                                gggCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAggg                                                                                              <  1:671211/150‑1 (MQ=255)
                                                     gaCCGTGAGCCAGTCGGCCAAGGGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGGAGGGTTGAc                                                                                         <  2:126262/149‑1 (MQ=255)
                                                      aCCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACg                                                                                        >  2:1167494/1‑150 (MQ=255)
                                                        cGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCCGGGTTGACGTg                                                                                      >  1:392072/1‑150 (MQ=255)
                                                           gAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCAc                                                                                   <  2:1016983/150‑1 (MQ=255)
                                                              ccAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCAcgcg                                                                                >  1:781248/1‑150 (MQ=255)
                                                               cAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCgg                                                                               <  2:1169878/150‑1 (MQ=255)
                                                                 gTCGACCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGc                                                                             <  2:880257/150‑1 (MQ=255)
                                                                     gCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAgg                                                                         >  2:710577/1‑150 (MQ=255)
                                                                     gCCAAGTGAATGCCGTAATACAGCGCGGCGCTGTACTCGGTGTCCTGCAATGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGTTAGGCAGGGTTGACGTGCACGCGGGCCAgg                                                                         <  2:619946/150‑1 (MQ=255)
                                                                      ccAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAggg                                                                        >  1:285324/1‑150 (MQ=255)
                                                                       cAAGTGAATGCCGTAATACTGCTCGGGGCTGTACTCGGTGTCCGCCACTGGCAGGAAGCCCGGGTGGTCGTAATCATAGAGGCCGGCGGCGCGGTTGACCAGGCGGGCGGTATTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGc                                                                       <  2:519988/150‑1 (MQ=255)
                                                                       cAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGc                                                                       <  2:512890/150‑1 (MQ=255)
                                                                          gTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGaa                                                                    <  2:892334/150‑1 (MQ=255)
                                                                           tGAATGCCGTAATACAGCTCGGCGCTGTACTCGTTGTCCTGCACTGGCAGGAAGCCCGTTCTGGCCCCGCTTTATAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGaaa                                                                   <  1:1034213/150‑1 (MQ=255)
                                                                                gctgTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGAt                                                              <  2:119667/147‑1 (MQ=255)
                                                                                 ccGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACGAGGCGGGCGTTCTTG‑‑CGATAGGCAGGATTGACGTGCACGCGGGCCAGGGCGAAATCGAtt                                                             >  1:901762/1‑149 (MQ=255)
                                                                                 ccGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCCGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCCCCACCGTCATGGT‑‑ACGTGCACGCGGGCCAGGGCGAAACCGATg                                                             >  2:452280/1‑150 (MQ=255)
                                                                                 ccGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTAGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATg                                                             <  2:40752/150‑1 (MQ=255)
                                                                                  cGTAATACAGCTCCGCGCTGTACTTGGTGTCCTGCACTGTCAGGAAGCCGTGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATgt                                                            <  2:254339/150‑1 (MQ=255)
                                                                                       tACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CCATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTcc                                                       >  1:82344/1‑150 (MQ=255)
                                                                                          aGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTg                                                    <  1:688985/150‑1 (MQ=255)
                                                                                          aGCTCGGCGCTGTACTCGGTGTCCTGCACTGGAAGGAAGCCAGGGCTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTg                                                    >  1:909274/1‑150 (MQ=255)
                                                                                            cTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGgc                                                  <  2:486014/150‑1 (MQ=255)
                                                                                             tCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGgcg                                                 >  1:1197826/1‑150 (MQ=255)
                                                                                                gcgcTGTACTCGGTGTCCTGCACTGGCAGGAAGGCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGAGGTCGTCCTTGGCGCgg                                              <  1:630758/150‑1 (MQ=255)
                                                                                                    tGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCAt                                          >  2:68106/1‑150 (MQ=255)
                                                                                                           ggTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAgg                                   >  2:203185/1‑150 (MQ=255)
                                                                                                           ggTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATACGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAgg                                   >  1:161402/1‑150 (MQ=255)
                                                                                                                 cTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGa                                                               <  2:654526/116‑1 (MQ=255)
                                                                                                                 cTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGa                                                               >  1:654526/1‑116 (MQ=255)
                                                                                                                 cTGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGCAAGGCCCttt                             >  2:67098/1‑150 (MQ=255)
                                                                                                                  tGCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTg                            >  1:1015292/1‑150 (MQ=255)
                                                                                                                  cgCACTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGTCAGGGTTGACGTGCACGCGGGCTAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTg                            <  2:48999/149‑1 (MQ=255)
                                                                                                                     aCTGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAAt                         <  1:832197/150‑1 (MQ=255)
                                                                                                                       tGGCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGATCCGTTGGTGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATAc                       <  1:83142/150‑1 (MQ=255)
                                                                                                                         gCAGGAAGCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCa                     <  2:351329/150‑1 (MQ=255)
                                                                                                                                gCCCGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgc              >  2:373904/1‑150 (MQ=255)
                                                                                                                                 cccGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgcc             <  1:267234/150‑1 (MQ=255)
                                                                                                                                 cccGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgcc             >  2:1109704/1‑150 (MQ=255)
                                                                                                                                  ccGGGATGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGAGGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATGCCAGTcctgcct            >  2:1058550/1‑150 (MQ=255)
                                                                                                                                   cGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGa                                                                                                       >  1:484048/1‑58 (MQ=255)
                                                                                                                                   cGGGTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGa                                                                                                       <  2:484048/58‑1 (MQ=255)
                                                                                                                                     ggTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCCTGca         <  1:383589/150‑1 (MQ=255)
                                                                                                                                     ggTTGTCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGATGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCCTGca         >  2:968118/1‑150 (MQ=255)
                                                                                                                                          tCGTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCCTGCACATAg    <  2:71531/150‑1 (MQ=255)
                                                                                                                                            gTAATCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTAGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCCTGCACATAGtt  >  1:897458/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCATCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTG‑‑CGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCCTGCACATAGTT  >  NC_002947/1647469‑1647758

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