Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,416,082 Δ1 bp 9.2% intergenic (+24/+95) PP_3901 → / ← PP_3902 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,416,0820C. 9.2% 200.7 / 14.4 65intergenic (+24/+95)PP_3901/PP_3902hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base . (3/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

TCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAG  >  NC_002947/4415942‑4416224
                                                                                                                                            |                                                                                                                                                
tCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCg                                                                                                                                        >  5:80095/1‑149 (MQ=255)
     ggAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCG‑C‑TCGAGGCGGGcaa                                                                                                                                   <  7:163323/149‑3 (MQ=255)
        aGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCttttt                                                                                                                                <  6:209473/149‑1 (MQ=255)
        aGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCttttt                                                                                                                                <  1:43505/149‑1 (MQ=255)
          aCTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCG‑C‑TCGAGGCGGGcaaggtgg                                                                                                                              <  1:227761/149‑8 (MQ=255)
              gggCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTgcgc                                                                                                                          <  3:114705/149‑1 (MQ=255)
              gggCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTgcgc                                                                                                                          <  8:330964/149‑1 (MQ=255)
                  gcgcTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATa                                                                                                                      <  1:182454/149‑1 (MQ=255)
                      tcgtcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCG‑C‑TCGA‑ggcgg                                                                                                                                      >  1:266221/1‑124 (MQ=255)
                      tcgtcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCG‑C‑TCGA‑ggcgg                                                                                                                                      <  2:266221/128‑5 (MQ=255)
                        gtcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGTGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGaaaa                                                                                                                 <  2:199797/148‑1 (MQ=255)
                            aGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGaa                                                                                                            <  6:80095/149‑1 (MQ=255)
                              gCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACGAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAAtt                                                                                                          <  1:2527/149‑1 (MQ=255)
                                cATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTAt                                                                                                        >  8:52212/1‑149 (MQ=255)
                                    gAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATggg                                                                                                     >  8:299702/1‑148 (MQ=255)
                                         cACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTacc                                                                                                 >  2:191164/1‑147 (MQ=255)
                                         cACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTacc                                                                                                 <  1:191164/147‑1 (MQ=255)
                                          aCCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTaccac                                                                                               >  5:354731/1‑148 (MQ=255)
                                            cGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGTGGCGTTCTTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCAt                                                                                            <  3:268106/149‑1 (MQ=255)
                                                ggCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACa                                                                                        <  1:248075/149‑1 (MQ=255)
                                                 gCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACAt                                                                                       >  1:70806/1‑149 (MQ=255)
                                                      aGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTAt                                                                                   <  2:148527/148‑1 (MQ=255)
                                                       gCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTAt                                                                                   >  6:4559/1‑147 (MQ=255)
                                                         cTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAg                                                                                >  8:71117/1‑148 (MQ=255)
                                                                 tGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGaaaa                                                                                                                 >  8:7826/1‑107 (MQ=255)
                                                                 tGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGaaaa                                                                                                                 <  7:7826/107‑1 (MQ=255)
                                                                   aaGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                      <  7:299702/148‑1 (MQ=255)
                                                                          tgGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGt                                                                                                               >  6:43974/1‑100 (MQ=255)
                                                                          tgGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGt                                                                                                               <  5:43974/100‑1 (MQ=255)
                                                                               gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAatat                                                          >  4:149068/1‑148 (MQ=255)
                                                                                cgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                        <  7:52212/149‑1 (MQ=255)
                                                                                cgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                        <  2:70806/149‑1 (MQ=255)
                                                                                     gtgtGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTTGGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGca                                                     >  1:102599/1‑148 (MQ=255)
                                                                                        tgACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCAGCC‑TCGA‑GCGGGCTTTTTTAGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactc                                                <  8:238285/149‑1 (MQ=255)
                                                                                        tgAACTTGGAAAACCCGCTAGAAGTCTTTTTTCGACAAACCCACCTTGCCCGCC‑TCCA‑GCGGGCTTTTTTTGCGGCATAGAAAATTGTATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactc                                                <  5:125332/149‑1 (MQ=255)
                                                                                         gACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactca                                               >  3:227636/1‑149 (MQ=255)
                                                                                          aCCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCAACCTTGCACGCC‑TAGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAAACCATTGACACAAACTATGAGcagactcc                                               >  5:98733/1‑147 (MQ=255)
                                                                                           gcTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACATTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatac                                                  <  8:163538/143‑1 (MQ=255)
                                                                                           ccTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactaga                                              >  7:163538/1‑145 (MQ=255)
                                                                                             tGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAat                                                            >  3:210404/1‑132 (MQ=255)
                                                                                             tGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAat                                                            <  4:210404/132‑1 (MQ=255)
                                                                                                    aGCCGCTAGAAGACTGATTTCTACTAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTTAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGc                                                      <  3:190988/131‑1 (MQ=255)
                                                                                                     gCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCa                                    >  8:35394/1‑148 (MQ=255)
                                                                                                        gCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCAc                                >  2:78682/1‑149 (MQ=255)
                                                                                                                aCTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                        >  4:120961/1‑117 (MQ=255)
                                                                                                                aCTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                        <  3:120961/117‑1 (MQ=255)
                                                                                                                 cTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGa                        <  5:4559/148‑1 (MQ=255)
                                                                                                                  tGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGaca                      >  3:333687/1‑149 (MQ=255)
                                                                                                                  tGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGaca                      >  6:114822/1‑149 (MQ=255)
                                                                                                                  tGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGaca                      >  1:27667/1‑149 (MQ=255)
                                                                                                                    aTTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACAcc                    >  7:268537/1‑149 (MQ=255)
                                                                                                                        cGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGc                 <  1:78682/148‑1 (MQ=255)
                                                                                                                        cGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGc                 <  6:354731/148‑1 (MQ=255)
                                                                                                                         gACAACCCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGc                 >  2:304300/1‑147 (MQ=255)
                                                                                                                         gACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCAACATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGa                            >  6:73710/1‑136 (MQ=255)
                                                                                                                                ccACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGcaaca         >  6:334376/1‑148 (MQ=255)
                                                                                                                                ccACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGcaaca         >  8:184974/1‑148 (MQ=255)
                                                                                                                                 cACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa       <  2:27667/149‑1 (MQ=255)
                                                                                                                                 cACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa       <  4:227636/149‑1 (MQ=255)
                                                                                                                                   ccTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCAGACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAgg     >  5:309095/1‑149 (MQ=255)
                                                                                                                                   aaaagCCCG‑C‑TCGAGGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAgg     >  3:94495/5‑149 (MQ=255)
                                                                                                                                   aaaagCCCG‑C‑TCGAGGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAgg     >  3:99471/5‑149 (MQ=255)
                                                                                                                                     ttGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatact                                                 >  6:228686/1‑103 (MQ=255)
                                                                                                                                     ttGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatact                                                 <  5:228686/103‑1 (MQ=255)
                                                                                                                                      tGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAg  >  5:19278/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                
TCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCC‑TCGA‑GCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAG  >  NC_002947/4415942‑4416224

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: ‑