Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,480,757 A→T 5.3% W169R (TGG→AGG)  PP_3090 ← OmpA domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,480,7570AT 5.3% 239.0 / 4.1 75W169R (TGG→AGG) PP_3090OmpA domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base T (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GCTGACGCTGACGATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAGGTACCAGGGCAGGGCATACAGGGC  >  NC_002947/3480608‑3480904
                                                                                                                                                     |                                                                                                                                                     
gctgacgctgacGATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCa                                                                                                                                                       >  2:124673/1‑150 (MQ=255)
  tgacgctgacgATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTACTTGGGGGCGCTGAAGGAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGt                                                                                                                                                    >  1:525902/1‑150 (MQ=255)
  caacgctgacgATCAGGCCGTTAAGGGGACGGCGTGGGCGGTGCGGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCAGGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGt                                                                                                                                                    <  1:689367/148‑1 (MQ=255)
   gacgctgacgATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTc                                                                                                                                                   >  2:1063826/1‑150 (MQ=255)
      gctgacgATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCa                                                                                                                                                >  1:1318529/1‑150 (MQ=255)
       ctgacgATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAg                                                                                                                                               <  2:1205030/150‑1 (MQ=255)
         gacgATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGtt                                                                                                                                             <  2:1163076/150‑1 (MQ=255)
            gATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTgcg                                                                                                                                         >  1:961760/1‑150 (MQ=255)
                         aGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAgg                                                                                                                            >  1:689118/1‑150 (MQ=255)
                          ggggACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGc                                                                                                                           >  2:566782/1‑150 (MQ=255)
                           gggACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGcc                                                                                                                          >  1:54674/1‑150 (MQ=255)
                            ggccGGCGCTGGCTGTGCTGGCGCAGCAGTTGCAAAAACCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGccc                                                                                                                         <  1:152549/147‑1 (MQ=255)
                                    cTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCg                                                                                                                 <  2:961760/150‑1 (MQ=255)
                                    cTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCg                                                                                                                 >  2:1160827/1‑150 (MQ=255)
                                     tGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCTGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGc                                                                                                                >  1:845116/1‑150 (MQ=255)
                                     tGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCTGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGc                                                                                                                >  1:845513/1‑150 (MQ=255)
                                             cTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTc                                                                                                        <  1:1020490/150‑1 (MQ=255)
                                                    agcagTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCg                                                                                                 <  2:525902/150‑1 (MQ=255)
                                                         ttGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATa                                                                                            >  1:243541/1‑150 (MQ=255)
                                                          tGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAg                                                                                           <  2:1076823/150‑1 (MQ=255)
                                                              aaaaGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCcg                                                                                                                        <  2:681902/117‑1 (MQ=255)
                                                              aaaaGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCcg                                                                                                                        >  1:681902/1‑117 (MQ=255)
                                                              aaaaGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAggggg                                                                                       >  2:679164/1‑149 (MQ=255)
                                                               aaaGCCGCGCCACTTGCCGGCTTCCGCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGaa                                                                                      <  2:54674/150‑1 (MQ=255)
                                                                      cgcCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAgg                                                                               <  1:1160827/150‑1 (MQ=255)
                                                                        ccACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCCGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGcc                                                                             >  2:130113/1‑150 (MQ=255)
                                                                          aCTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCg                                                                           >  1:917375/1‑150 (MQ=255)
                                                                            ttGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCg                                                                                                 >  1:734793/1‑126 (MQ=255)
                                                                            ttGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCg                                                                                                 >  1:735050/1‑126 (MQ=255)
                                                                            ttGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCg                                                                                                 <  2:734793/126‑1 (MQ=255)
                                                                            ttGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCg                                                                                                 <  2:735050/126‑1 (MQ=255)
                                                                            ttGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGAc                                                                         >  1:1235192/1‑150 (MQ=255)
                                                                            ctGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAg                                                                                     >  2:1411962/2‑138 (MQ=255)
                                                                            ctGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAg                                                                                     <  1:1411962/137‑1 (MQ=255)
                                                                                cGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTGGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCa                                                                     <  2:451388/150‑1 (MQ=255)
                                                                                  ccTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCCATTCGTGATGTGTTCCACGGAACGGCCGTCACCTGCCGCCCATGCACGAGACCCGCCCGCCTTCGCGCTCTGCATAGGGGAAGTTCAGGCCCGACTGCAGt                                                                   <  1:679164/149‑1 (MQ=255)
                                                                                    ttCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTat                                                                 >  1:960380/1‑150 (MQ=255)
                                                                                     tCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTatc                                                                <  1:566782/150‑1 (MQ=255)
                                                                                      ccTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTatca                                                               <  2:689118/150‑1 (MQ=255)
                                                                                       cTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTatcat                                                              <  2:1318529/150‑1 (MQ=255)
                                                                                        tCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCGGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGTGGAAGTTCAGGCCCGACTGCAGTatcatc                                                             <  1:235428/150‑1 (MQ=255)
                                                                                        tCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCcgc                                                                                                                       <  2:646376/92‑1 (MQ=255)
                                                                                        tCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCcgc                                                                                                                       >  1:646376/1‑92 (MQ=255)
                                                                                        tCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCAc                                                                                                              >  1:806270/1‑101 (MQ=255)
                                                                                        tCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCAc                                                                                                              <  2:806270/101‑1 (MQ=255)
                                                                                          gtggCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAACTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCACACGCCTTCGCGCCCGGCAAAGGGGAAGACCAGTCACGACTGCAGTATCATCGt                                                           >  2:726045/3‑150 (MQ=255)
                                                                                               tGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTatcatc                                                             <  1:717658/143‑1 (MQ=255)
                                                                                               tGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTatcatc                                                             >  2:717658/1‑143 (MQ=255)
                                                                                                 ttCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTgcgc                                                                                                                                        <  1:828157/66‑1 (MQ=255)
                                                                                                 ttCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTgcgc                                                                                                                                        >  2:828157/1‑66 (MQ=255)
                                                                                                 ttCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTg                                                    >  2:657425/1‑150 (MQ=255)
                                                                                                     tGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTc                                                                                                        <  2:874565/94‑1 (MQ=255)
                                                                                                     tGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTc                                                                                                        >  1:874565/1‑94 (MQ=255)
                                                                                                              cGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAg                                                                                >  1:1266675/1‑109 (MQ=255)
                                                                                                              cGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAg                                                                                <  2:1266675/109‑1 (MQ=255)
                                                                                                              cGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTTCATTATTATCTTGCTCTTGCCGTCGGCAGGCt                                       <  2:1361021/150‑1 (MQ=255)
                                                                                                               ggCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAgg                                                                                                                            >  2:771234/1‑64 (MQ=255)
                                                                                                               ggCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAgg                                                                                                                            <  1:771234/64‑1 (MQ=255)
                                                                                                                gCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTgg                                     <  2:243541/150‑1 (MQ=255)
                                                                                                                       ccAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGa                              <  1:130113/150‑1 (MQ=255)
                                                                                                                        cagcagTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGAc                             <  1:1063826/150‑1 (MQ=255)
                                                                                                                          gcagTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGc                                                   <  2:913455/126‑1 (MQ=255)
                                                                                                                          gcagTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGc                                                   >  1:913455/1‑126 (MQ=255)
                                                                                                                           cagTACCGCTTCAGAGCTGAAGAACCTGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAg                          <  2:845116/150‑1 (MQ=255)
                                                                                                                           cagTACCGCTTCAGAGCTGAAGAACCTGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAg                          <  2:845513/150‑1 (MQ=255)
                                                                                                                                ccGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATACGGGGAGGTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAGGTAcc                     >  1:449256/1‑150 (MQ=255)
                                                                                                                                       agagCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAg                                                                                           >  1:1353663/1‑73 (MQ=255)
                                                                                                                                       agagCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAg                                                                                           <  2:1353663/73‑1 (MQ=255)
                                                                                                                                         agCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAGGTACcagggcaggg            >  1:1271805/1‑150 (MQ=255)
                                                                                                                                           cTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAGGTACCAgggcagggca          >  2:1004968/1‑150 (MQ=255)
                                                                                                                                             gaagaaCCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAGGTACCAGGGCAGGGCATa        >  1:241131/1‑150 (MQ=255)
                                                                                                                                              aagaaCCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCa                                                                     >  1:445481/1‑88 (MQ=255)
                                                                                                                                              aagaaCCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCa                                                                     <  2:445481/88‑1 (MQ=255)
                                                                                                                                                  aCCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAGGTACCAGGGCAGGGCATACAggg   >  1:614143/1‑150 (MQ=255)
                                                                                                                                                   ccAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCCCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCCGCGGTAGGCTGGCCAATGACCAGGTACCAGGGCAGGGCATACAGGGc  >  1:536850/1‑150 (MQ=255)
                                                                                                                                                   ccAGTCGCAGTTGCGCGTGCCGCCCAGGCACGGGCAGCGACTGCGCTAGCGCCACTGAGATGGGAAGTGCACGCGCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTTGCCAAGTGCCAGGTACCAGGGCAGGGCATACAggga  >  1:248411/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGACGCTGACGATCAAGCCGTTAAGGGGACGGCGCTGGCGGTGCTGGCGCAGCAGTTGCAAAAAGCCGCGCCACTTGCCGGCTTCCTCGGGGCTGTTCATGTAACGCCCGGCAGTGTCCAGCAGTACCGCTTCAGAGCTGAAGAACCAGTCGCAGTTGCGCGTGCCGCCCAGGCCCGCGACCCGCACGCCTTCGCGCTCGGCATAGGGGAAGTTCAGGCCCGACTGCAGTATCATCGTGCTCTTGCCGGCGGCAGGCTGGCCAATGACCAGGTACCAGGGCAGGGCATACAGGGC  >  NC_002947/3480608‑3480904

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