Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,938 T→A 5.2% S207T (TCC→ACC)  ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,9380TA 5.2% 285.3 / 5.1 77S207T (TCC→ACC) ycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base A (2/2);  total (40/37)
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.31e-01

AGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGCCTGGGCGCTGG  >  NC_000913/964792‑965081
                                                                                                                                                  |                                                                                                                                               
aGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATc                                                                                                                                                <  1:766562/148‑1 (MQ=255)
   cGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGACGGATAACCGtt                                                                                                                                            <  4:1023764/148‑2 (MQ=255)
             ccGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTAt                                                                                                                                                 >  2:416183/1‑134 (MQ=255)
             ccGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTAt                                                                                                                                                 <  1:416183/134‑1 (MQ=255)
                     tGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaaaa                                                                                                                           >  4:535497/1‑148 (MQ=255)
                      gCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAat                                                                                                                          >  1:277417/1‑148 (MQ=255)
                         ggCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCaa                                                                                                                                       <  4:375086/132‑1 (MQ=255)
                         ggCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCaa                                                                                                                                       >  3:375086/1‑132 (MQ=255)
                           cAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATc                                                                                                                                                >  2:460299/1‑121 (MQ=255)
                           cAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATc                                                                                                                                                <  1:460299/121‑1 (MQ=255)
                               tgtCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGGACGGATAACCGttcc                                                                                                                                          >  3:1046509/1‑120 (MQ=255)
                               tgtCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGGACGGATAACCGTtcc                                                                                                                                          <  4:1046509/123‑1 (MQ=255)
                                  cATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGaa                                                                                                                                   <  1:1072649/127‑1 (MQ=255)
                                  cATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGaa                                                                                                                                   >  2:1072649/1‑127 (MQ=255)
                                  cATTGAACCGAATTGTAGCCTGCGTACACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCaa                                                                                                                                       >  4:286367/1‑123 (MQ=255)
                                  cATTGAACCGAATTGTAGCCTGCGTACACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCaa                                                                                                                                       <  3:286367/123‑1 (MQ=255)
                                    ttGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATAc                                                                                                            >  1:585938/1‑148 (MQ=255)
                                     tGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTc                                                                                                                                            >  4:1107835/1‑115 (MQ=255)
                                     tGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTc                                                                                                                                            <  3:1107835/115‑1 (MQ=255)
                                     tGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATAc                                                                                                            <  3:198049/147‑1 (MQ=255)
                                     tGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATAc                                                                                                            >  4:198049/1‑147 (MQ=255)
                                        aCCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaaa                                                                                                                            >  2:868175/1‑128 (MQ=255)
                                        aCCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaaa                                                                                                                            <  1:868175/128‑1 (MQ=255)
                                        aCCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGcgcg                                                                                                                <  1:67399/140‑1 (MQ=255)
                                        aCCGAATTGTAGCATGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGcgcg                                                                                                                >  2:67399/1‑140 (MQ=255)
                                          cGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGaa                                                                                                                                   <  3:455131/119‑1 (MQ=255)
                                          cGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGaa                                                                                                                                   >  1:453224/1‑119 (MQ=255)
                                          cGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGaa                                                                                                                                   <  2:453224/119‑1 (MQ=255)
                                          cGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGaa                                                                                                                                   >  4:455131/1‑119 (MQ=255)
                                           gAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaaa                                                                                                                            <  4:22185/125‑1 (MQ=255)
                                           gAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaaa                                                                                                                            >  3:22185/1‑125 (MQ=255)
                                           gAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACg                                                                                                     >  4:967273/1‑148 (MQ=255)
                                               tGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGgcgc                                                                                                 >  4:290215/1‑148 (MQ=255)
                                                 tAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGt                                                                                                                                             >  2:956889/1‑102 (MQ=255)
                                                 tAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGt                                                                                                                                             <  1:956889/102‑1 (MQ=255)
                                                    ccTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTa                                                                                            <  3:967273/148‑1 (MQ=255)
                                                      tGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAAt                                                                                          <  2:277417/148‑1 (MQ=255)
                                                       gCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATg                                                                                         <  3:535497/148‑1 (MQ=255)
                                                       gCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATg                                                                                         <  3:290215/148‑1 (MQ=255)
                                                        cGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATgg                                                                                        <  2:585938/148‑1 (MQ=255)
                                                                gTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATc                                                                                                                                >  2:863331/1‑100 (MQ=255)
                                                                gTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATc                                                                                                                                <  1:863331/100‑1 (MQ=255)
                                                                          tCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCAt                                                                                                                     <  1:38397/101‑1 (MQ=255)
                                                                          tCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCAt                                                                                                                     >  2:38397/1‑101 (MQ=255)
                                                                          tCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTg                                                                                                          <  3:318037/112‑1 (MQ=255)
                                                                          tCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTg                                                                                                          >  4:318037/1‑112 (MQ=255)
                                                                              tACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGTAACGGCGCATTTAATGGCGATATCGGGATTGCACGTCGc                                                                  >  4:711284/1‑148 (MQ=255)
                                                                               acaaacaaCGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCt                                                                 >  3:956989/1‑148 (MQ=255)
                                                                                 aaacaaCGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGgcg                                                                                                  <  2:303099/113‑1 (MQ=255)
                                                                                 aaacaaCGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGgcg                                                                                                  >  1:303099/1‑113 (MQ=255)
                                                                                    caacGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTgc                                                            >  3:1075601/1‑148 (MQ=255)
                                                                                    caacGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTgc                                                            >  3:890320/1‑148 (MQ=255)
                                                                                         cTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATc                                                                                                                       <  2:857510/84‑1 (MQ=255)
                                                                                         cTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATc                                                                                                                       >  1:857510/1‑84 (MQ=255)
                                                                                                 cTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATa                                                                                                             >  3:895063/1‑86 (MQ=255)
                                                                                                 cTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATa                                                                                                             <  4:895063/86‑1 (MQ=255)
                                                                                                 cTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGcac                                                                       >  2:327854/1‑124 (MQ=255)
                                                                                                 cTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGcac                                                                       <  1:327854/124‑1 (MQ=255)
                                                                                                   aTCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGa                                                                             >  1:1023404/1‑116 (MQ=255)
                                                                                                   aTCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGa                                                                             <  2:1023404/116‑1 (MQ=255)
                                                                                                       gTGGAATGTGGTTATTCTTGGTTTAGGTATGGAGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGACCGGCGCATTTAATGGCGATATCGGGACTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGa                                         <  3:711284/148‑1 (MQ=255)
                                                                                                          gAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATc                                                                                                                       <  4:1032045/67‑1 (MQ=255)
                                                                                                          gAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATc                                                                                                                       >  3:1032045/1‑67 (MQ=255)
                                                                                                             tGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCAtt                                   <  4:956989/148‑1 (MQ=255)
                                                                                                                         tGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATg                                                                                         >  3:262445/1‑82 (MQ=255)
                                                                                                                         tGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATg                                                                                         <  4:262445/82‑1 (MQ=255)
                                                                                                                            tttAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGt                             >  1:915592/1‑139 (MQ=255)
                                                                                                                            tttAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGt                             <  2:915592/139‑1 (MQ=255)
                                                                                                                                   aTGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGc             <  4:1075601/148‑1 (MQ=255)
                                                                                                                                   aTGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGc             <  4:890320/148‑1 (MQ=255)
                                                                                                                                     ttggggACGGATAACCGTTCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGCCt           >  4:22919/3‑148 (MQ=255)
                                                                                                                                       gggAACGGTTATCCGTCCCCACAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTCATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGCCTgg         >  1:222342/1‑148 (MQ=255)
                                                                                                                                       gggAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGCCTgg         >  1:450586/1‑148 (MQ=255)
                                                                                                                                         gAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTg                                                                          <  3:1007365/81‑1 (MQ=255)
                                                                                                                                         gAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTg                                                                          >  4:1007365/1‑81 (MQ=255)
                                                                                                                                           aCGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGCCTGGgcgc     >  1:18639/1‑148 (MQ=255)
                                                                                                                                              gTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTgctg                                                  <  4:648174/100‑1 (MQ=255)
                                                                                                                                              gTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTgctg                                                  >  3:648174/1‑100 (MQ=255)
                                                                                                                                              gTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGCCTGGGCGCTgg  >  4:545141/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCTGCCTGGGCGCTGG  >  NC_000913/964792‑965081

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