Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 154,908 G→C 6.8% L130V (CTG→GTG) ‡ yjbB ← Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947154,9080GC 6.8% 224.4 / 4.0 72L130V (CTG→GTG) ‡yjbBNa+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base C (3/2);  total (36/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 = 1.00e+00

CCCGACCAGAGCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATCGCCAACGCCGGG  >  NC_002947/154762‑155056
                                                                                                                                                  |                                                                                                                                                    
cccGACCAGAGCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCACACCCGGCCCAGtt                                                                                                                                                    >  7:286809/1‑149 (MQ=255)
   gACCAGAGCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTggc                                                                                                                                                 <  1:252610/149‑1 (MQ=255)
        gagCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTggc                                                                                                                                                 >  4:245117/1‑144 (MQ=255)
        gagCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTggc                                                                                                                                                 <  3:245117/144‑1 (MQ=255)
         agCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGccc                                                                                                                                           <  4:15245/149‑1 (MQ=255)
         agCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGccc                                                                                                                                           <  5:351077/149‑1 (MQ=255)
         agCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGccc                                                                                                                                           <  6:290525/149‑1 (MQ=255)
             tCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCgtgt                                                                                                                                       >  7:99996/1‑149 (MQ=255)
               gAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGTCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTc                                                                                                                                      <  2:4923/148‑1 (MQ=255)
               gAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTAGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCt                                                                                                                                     >  1:267954/1‑149 (MQ=255)
                 gCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGc                                                                                                                                   <  7:71737/149‑1 (MQ=255)
                  cATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCt                                                                                                                                  >  7:296113/1‑149 (MQ=255)
                      atgTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTgc                                                                                                                               <  8:232210/148‑1 (MQ=255)
                           gCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGaaa                                                                                                                         <  5:30104/149‑1 (MQ=255)
                           gCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGa                                                                                                                           >  6:353019/1‑147 (MQ=255)
                               gTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGgaa                                                                                                                     <  1:112700/149‑1 (MQ=255)
                               gTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGgaa                                                                                                                     <  8:296113/149‑1 (MQ=255)
                                tCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGgaag                                                                                                                    <  2:267954/149‑1 (MQ=255)
                                   gTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGaagaag                                                                                                                 <  7:118459/149‑1 (MQ=255)
                                    tGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGgaag                                                                                                                    >  4:194973/1‑145 (MQ=255)
                                      aGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATc                                                                                                              <  3:106156/149‑1 (MQ=255)
                                       ggCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCa                                                                                                             >  5:330460/1‑149 (MQ=255)
                                        gCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCAc                                                                                                            >  6:39728/1‑149 (MQ=255)
                                        gCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCAACTGggccg                                                                                                                                               >  5:151951/1‑114 (MQ=255)
                                        gCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCAACTGggccg                                                                                                                                               <  6:151951/114‑1 (MQ=255)
                                           aaCAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGc                                                                                                          >  4:104344/1‑148 (MQ=255)
                                            aCAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCCCGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTGGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGTTCACGcc                                                                                                         >  8:167091/1‑148 (MQ=255)
                                              aGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGAGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAg                                                                                                      >  1:105116/1‑149 (MQ=255)
                                                 aCTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAgg                                                                                                     >  5:127326/1‑147 (MQ=255)
                                                     tCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCATGGAAGAt                                                                                                <  5:39728/148‑1 (MQ=255)
                                                       aCGCCTTGGGCTTGGGTGATCGGCGCAGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCACACCCCGCCCAGTTGGCCGGACCGTGTCTGCTTGCGCGAAAGGAAGAACATAACGCCAAGGAAGATcag                                                                                             >  6:334109/1‑149 (MQ=255)
                                                           cttgggcttgggTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGgaagaa                                                                                                                  >  7:349702/1‑124 (MQ=255)
                                                           attgggcttgggTGATCGGCTCCGCTGCCTGCACGATAAGTTCGATTTCCAGTATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTTGCCGGCCCGTGTCTGCTTGCGCGAAAGgaagaa                                                                                                                  <  8:349702/123‑1 (MQ=255)
                                                              gggcttgggTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTg                                                                                        <  7:198194/147‑1 (MQ=255)
                                                               ggcttgggTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGc                                                                                      <  5:353019/148‑1 (MQ=255)
                                                                  ttgggTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTGGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGa                                                                                    >  3:294663/1‑147 (MQ=255)
                                                                    gggTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGc                                                                                >  4:171138/1‑149 (MQ=255)
                                                                       tGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCACACCCGGCCCAGTTGGTCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAg                                                                             <  8:286809/149‑1 (MQ=255)
                                                                        gATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAg                                                                             >  5:242281/1‑148 (MQ=255)
                                                                         aTCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGa                                                                           <  7:185186/149‑1 (MQ=255)
                                                                            ggCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAg                                                                        >  4:361369/1‑149 (MQ=255)
                                                                             gcgcCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAgg                                                                       <  3:194973/149‑1 (MQ=255)
                                                                              cgccgcTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGt                                                                      <  8:99996/149‑1 (MQ=255)
                                                                                  gcTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGaa                                                                  <  6:242281/149‑1 (MQ=255)
                                                                                      cctgcACGATAAGTTCGAGTGCCAGGATGATGAGGACAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGt                                                               <  5:334109/148‑1 (MQ=255)
                                                                                       ctgcACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCa                                                             <  3:171138/149‑1 (MQ=255)
                                                                                                  gTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGc                                                   >  6:150118/1‑148 (MQ=255)
                                                                                                  gTTCCAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGc                                                   >  1:211861/1‑148 (MQ=255)
                                                                                                   ttCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTc                                                              <  5:152318/136‑1 (MQ=255)
                                                                                                   ttCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTc                                                              >  6:152318/1‑136 (MQ=255)
                                                                                                            cAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCgg                                        >  5:96062/1‑149 (MQ=255)
                                                                                                            cAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCgg                                        <  3:104344/149‑1 (MQ=255)
                                                                                                            cAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCgg                                        >  7:197616/1‑149 (MQ=255)
                                                                                                             aGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGt                                       <  6:330460/149‑1 (MQ=255)
                                                                                                              ggATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTa                                      <  3:361369/149‑1 (MQ=255)
                                                                                                                  gatCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCGGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCca                                  >  1:94713/1‑149 (MQ=255)
                                                                                                                   atCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCCCGAAAGGAAGAAGCTCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGACGAAGGTCAGTACCCCGGCCATCAGTGCGGTACCca                                  >  2:214052/1‑148 (MQ=255)
                                                                                                                    tCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCcaca                                >  8:6028/1‑149 (MQ=255)
                                                                                                                      aGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATc                              >  6:66796/1‑149 (MQ=255)
                                                                                                                       ggCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCg                             >  3:286872/1‑149 (MQ=255)
                                                                                                                       ggCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCAAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCg                             <  2:234337/149‑1 (MQ=255)
                                                                                                                                cGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGc                    <  2:105116/149‑1 (MQ=255)
                                                                                                                                      ccACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATc              >  6:57728/1‑149 (MQ=255)
                                                                                                                                      ccACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATc              <  2:211861/149‑1 (MQ=255)
                                                                                                                                      ccACCCGGCCCAGTTCGCCGGCCCGTGTCTGATTGAGCGAAAGGAAGAAGATCACGCAAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGATGGTCAGTACCCGGGCCAACAGTGCGGAACCCACATCGGCACCCAGCatcatc              >  2:119804/1‑149 (MQ=255)
                                                                                                                                       gggccGGCCAACTGGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTCCGGTCCCCACATCGGCCCCAAGCATGAACg             >  5:115666/4‑149 (MQ=255)
                                                                                                                                         gccGGCCAACTGGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGc                    >  5:259797/2‑140 (MQ=255)
                                                                                                                                         gccGGCCAACTGGGCAGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGc                    <  6:259797/139‑1 (MQ=255)
                                                                                                                                            ggCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATCGCCaa         >  3:298204/1‑148 (MQ=255)
                                                                                                                                            ggCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATCGCCaa         <  4:298204/148‑1 (MQ=255)
                                                                                                                                            ggCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATCGCCAAc        <  4:294663/149‑1 (MQ=255)
                                                                                                                                               ccAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATCGCCAACGcc     >  7:113825/1‑149 (MQ=255)
                                                                                                                                                  gTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATCGCCAACGCCggg  <  5:150118/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCCGACCAGAGCGTCGAGCATGATGTCGCCGGTCAGTGAGGCAAACAGCACTTTCACGCCTTGGGCTTGGGTGATCGGCGCCGCTGCCTGCACGATAAGTTCGAGTGCCAGGATGATCAGGCCAAGCCCGATGCCCACCCGGCCCAGTTGGCCGGCCCGTGTCTGCTTGCGCGAAAGGAAGAAGATCACGCCAAGGAAGATCAGCAGTGGCGACAGCCAGGACAGGTCGAAGGTCAGTACCCGGGCCATCAGTGCGGTACCCACATCGGCACCCAGCATGATCGCCAACGCCGGG  >  NC_002947/154762‑155056

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