Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,236,466 G→T 5.6% A69S (GCT→TCT)  yjbF → extracellular polysaccharide production lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,236,4660GT 5.6% 253.6 / 3.8 72A69S (GCT→TCT) yjbFextracellular polysaccharide production lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base T (2/2);  total (38/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 = 5.20e-01

CTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCAC  >  NC_000913/4236328‑4236610
                                                                                                                                          |                                                                                                                                                
cTGGGCAATTCACTGTGGGACAGGCTGTTCGGCAGGCCAGGCGTACAGCTGACGGATGATGATATTCATAATATGCCCTACGCCAGCCAGTACATGCAGCTTACTGGCGGGCCGCAGTTATTTGTGGTGCTGGCGTTAGCTGAAGAc                                                                                                                                          <  8:327945/147‑1 (MQ=255)
 tGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACg                                                                                                                                         >  7:291899/1‑147 (MQ=255)
   ggCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGa                                                                                                                                       <  2:216566/147‑1 (MQ=255)
       aTTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGacaac                                                                                                                                   >  4:253604/1‑147 (MQ=255)
           aCTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACaaaa                                                                                                                               <  3:43226/147‑1 (MQ=255)
           aCTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACaaaa                                                                                                                               <  3:253604/147‑1 (MQ=255)
              gtgGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCTCTGAAGACGGACAACAAAAATgg                                                                                                                           >  2:81558/1‑148 (MQ=255)
              gtgGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATgg                                                                                                                           >  2:18492/1‑148 (MQ=255)
               tgGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATggg                                                                                                                          <  8:23592/148‑1 (MQ=255)
                    cAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCTCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACt                                                                                                                     <  4:103403/148‑1 (MQ=255)
                       tCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGa                                                                                                                                           <  3:83198/123‑1 (MQ=255)
                       tCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGa                                                                                                                                           >  4:83198/1‑123 (MQ=255)
                       tCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAg                                                                                                                  >  6:117394/1‑148 (MQ=255)
                        cTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAg                                                                                                                                            <  8:4856/121‑1 (MQ=255)
                        cTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAg                                                                                                                                            >  7:4856/1‑121 (MQ=255)
                          gTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACaaaa                                                                                                                               >  5:107763/1‑132 (MQ=255)
                          gTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACaaaa                                                                                                                               <  6:107763/132‑1 (MQ=255)
                                     cAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTActct                                                                                                    <  3:29509/148‑1 (MQ=255)
                                            aCAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACaaaa                                                                                                                               >  2:349411/1‑114 (MQ=255)
                                            aCAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACaaaa                                                                                                                               <  1:349411/114‑1 (MQ=255)
                                            aCAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTAcc                                                                                             >  7:167252/1‑148 (MQ=255)
                                            aCAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTAcc                                                                                             >  5:162829/1‑148 (MQ=255)
                                            aCAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTAcc                                                                                             >  7:161838/1‑148 (MQ=255)
                                            aCAGCTGACGGATGATGATAATCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTAcc                                                                                             >  7:173414/1‑148 (MQ=255)
                                               gCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCaa                                                                                          >  7:334964/1‑148 (MQ=255)
                                                cTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAAc                                                                                         <  6:162829/148‑1 (MQ=255)
                                                cTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAAc                                                                                         <  1:18492/148‑1 (MQ=255)
                                                cTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAAc                                                                                         <  8:161838/148‑1 (MQ=255)
                                                cTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAAc                                                                                         <  8:167252/148‑1 (MQ=255)
                                                cTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAAc                                                                                         <  8:334964/148‑1 (MQ=255)
                                                cTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGACGTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAAc                                                                                         <  1:382733/148‑1 (MQ=255)
                                                     ggATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGc                                                                                    <  8:291899/148‑1 (MQ=255)
                                                       atgatgatATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCg                                                                                  >  3:85990/1‑148 (MQ=255)
                                                                 tCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCg                                                                                  <  3:94090/138‑1 (MQ=255)
                                                                 tCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCg                                                                                  >  4:94090/1‑138 (MQ=255)
                                                                 tCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTg                                                                           >  1:66075/1‑145 (MQ=255)
                                                                  cAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAg                                                                        <  4:85990/147‑1 (MQ=255)
                                                                  cAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGa                                                                       >  8:36562/1‑148 (MQ=255)
                                                                                     gccagTACATGCAGCTTAATGGCGGTACGCAGTTATTTGTGGTGCTGGACTTCGCTGAAGACTGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATTGCCGTCTGGTGAAGACACTGCTCGGCGGCGACaa                                                    >  8:138342/1‑148 (MQ=255)
                                                                                          tACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCTCTGAAGACGGACAACAAAAATGGGTCACTCAGGATc                                                                                                              <  6:308963/85‑1 (MQ=255)
                                                                                          tACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCTCTGAAGACGGACAACAAAAATGGGTCACTCAGGATc                                                                                                              >  5:308963/1‑85 (MQ=255)
                                                                                           aCATGCAGCTTAATGGCGGGCCGCTGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGAt                                              >  5:140694/1‑148 (MQ=255)
                                                                                           aCATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAgg                                                                                                           <  8:388526/87‑1 (MQ=255)
                                                                                           aCATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAgg                                                                                                           >  7:388526/1‑87 (MQ=255)
                                                                                           aCATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCt                                                                               <  8:130070/115‑1 (MQ=255)
                                                                                           aCATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCt                                                                               >  7:130070/1‑115 (MQ=255)
                                                                                           aCATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGAt                                              >  4:81713/1‑148 (MQ=255)
                                                                                                cAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAg                                         >  3:184826/1‑148 (MQ=255)
                                                                                                  gCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCTCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTg                                       <  1:81558/148‑1 (MQ=255)
                                                                                                  gCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTg                                       <  3:81713/148‑1 (MQ=255)
                                                                                                  gCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTg                                       <  8:173414/148‑1 (MQ=255)
                                                                                                   cTTAATGGCGGGCCGCAGTTATTTTTGGTGCGGTCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCTACAACCTGattgaattga                                      <  7:138342/148‑1 (MQ=255)
                                                                                                      aaTGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAAt                                    >  6:126862/1‑147 (MQ=255)
                                                                                                          gCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCTGCGACAACCTGATTGAAGTGAATAACCt                               <  2:66075/148‑1 (MQ=255)
                                                                                                                  cAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGcc                        >  5:336005/1‑147 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAg                                                                        >  2:138478/1‑95 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAg                                                                        <  1:138478/95‑1 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACaa                                                    >  3:265926/1‑115 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACaa                                                    <  4:265926/115‑1 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGaa                                          >  3:240203/1‑125 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGaa                                          <  4:240203/125‑1 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGAccc                   >  1:185570/1‑148 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGAccc                   >  7:176103/1‑148 (MQ=255)
                                                                                                                      tATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGAccc                   >  3:155064/1‑148 (MQ=255)
                                                                                                                         ttGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCt                >  8:40785/1‑148 (MQ=255)
                                                                                                                                gCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAAc                                 >  6:87381/1‑124 (MQ=255)
                                                                                                                                gCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTaaa         >  8:185273/1‑148 (MQ=255)
                                                                                                                                gCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACAGGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAAc                                 <  5:87381/124‑1 (MQ=255)
                                                                                                                                gCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACAAGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTaaa         <  5:117394/148‑1 (MQ=255)
                                                                                                                                  tGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTa           >  8:96323/1‑144 (MQ=255)
                                                                                                                                  tGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTa           <  7:96323/144‑1 (MQ=255)
                                                                                                                                        tCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGcac  >  7:115663/1‑147 (MQ=255)
                                                                                                                                          |                                                                                                                                                
CTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCAC  >  NC_000913/4236328‑4236610

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