Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 222,958 Δ1 bp 100% coding (126/576 nt) gmhB → D,D‑heptose 1,7‑bisphosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913222,9580A.98.1% 222.4 / ‑1.7 52coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (22/29);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTA  >  NC_000913/222813‑223101
                                                                                                                                                 |                                                                                                                                               
ttatcaCGGAGCAATAAACGGCGGCGACGAGCGCACCCGCAATTTTTCTTGACCGTGATGCCCCCCTTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGGTATTGACGCCATGCGCGAGCTaaaa                                                                                                                                              <  2:489070/149‑1 (MQ=255)
 tattaAGGTGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAAAGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAATTATCGACGGTGATATTGACGCCATGCGCGAGCTaaaaa                                                                                                                                             >  1:455189/1‑149 (MQ=255)
  attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                            >  1:274512/1‑149 (MQ=255)
  attaAGGAGCTATAAAAGGGCGACAAGAGCGTACCCGCAATTTTTCTTGACCGAGAAGGCACCATTAATGTCGAACACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaa                                                                                                                                            <  2:849105/149‑1 (MQ=255)
    taAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                          <  2:48837/149‑2 (MQ=255)
    taAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                          <  1:431575/149‑2 (MQ=255)
     aaGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatg                                                                                                                                         <  2:88152/149‑3 (MQ=255)
     aaGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatg                                                                                                                                         <  1:885194/149‑3 (MQ=255)
     aaGCAGCTATAAAAGGGGGCGAAGAGCGTACCCACAATTTTTCTTGCCCGGGATGGCACCATTAATGTCGAACACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatg                                                                                                                                         <  2:1001188/149‑3 (MQ=255)
       ggcgcaataAAAGGCGCCGCAGAGCGTACCCGCAATTTTTCTTGACCGTCATGCCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGCCGCTGTTATTGACGCCATGCGCGAGCTAAAAAAAAggg                                                                                                                                       <  2:847752/143‑1 (MQ=255)
       ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGCCCGTGATGGCACCATTAATGTCGATCACGGCTAAGTCCATCAGATCGACAACTTTGATTTTATCGACGGTGTTATTGAAGCCATGCCCGAGCTaaaaaaatgtg                                                                                                                                       >  2:666133/1‑145 (MQ=255)
             ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGCTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgc                                                                                                                                >  1:826925/1‑149 (MQ=255)
             ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgc                                                                                                                                <  2:950205/149‑1 (MQ=255)
               aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcgc                                                                                                                              <  2:274512/149‑1 (MQ=255)
               aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACAGTGTGATTGACGCCATGCGCGAGCT‑AAAAAAATGGCCTTTgcgc                                                                                                                              >  1:607385/1‑149 (MQ=255)
                   ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGAACGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGATATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggt                                                                                                                          <  1:627257/149‑1 (MQ=255)
                     tGGCGAAGTCCGTACCCGCAATTTTTCTATACCGGGATTGCACCATAAATGTCGATCACGGCTATGACGATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGACCT‑AAAAAATTGGGCTTTGCGCtggtgg                                                                                                                        <  2:607385/149‑1 (MQ=255)
                          aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCCCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGGTATTGCCGACATGCGCGAGCT‑AAAAAATTGGGCTTTGCGCTGGTGgtagta                                                                                                                   >  2:1005307/1‑149 (MQ=255)
                           agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                  <  2:369146/149‑1 (MQ=255)
                               cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                              >  2:463672/1‑149 (MQ=255)
                                    ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGt                                                                                                         <  1:24162/149‑1 (MQ=255)
                                     cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                        >  1:425032/1‑149 (MQ=255)
                                       cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                      <  2:2901/149‑1 (MQ=255)
                                           ttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAt                                                                                                  >  1:1014570/1‑149 (MQ=255)
                                                 tGACCGTGATGGCACCATTAATGTCGAGCACGGCTATGTCCAAGAGATCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCACCCAGGCTGGCATTGCTcg                                                                                            >  1:583196/1‑149 (MQ=255)
                                                    ccGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTGGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGAGTTGCGCTGGTGGTATTAACCAACCAGTCTTGCATTGCTCGCgg                                                                                         >  2:601795/1‑149 (MQ=255)
                                                      gTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTa                                                                                       >  2:287022/1‑149 (MQ=255)
                                                            gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTa                                                                                 <  2:365903/149‑1 (MQ=255)
                                                            gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCGGGTGGTAGTAACCAACCAGTCTGGAATTGCTCGCGGTAAATTTa                                                                                 <  2:370059/149‑1 (MQ=255)
                                                             cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                <  1:640768/149‑1 (MQ=255)
                                                              aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATAGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTATCCAACCAGTCTGGCATAGATCGCGGTAAATTTAcc                                                                               >  1:167084/1‑149 (MQ=255)
                                                                cATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGCCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGa                                                                             >  2:212279/1‑149 (MQ=255)
                                                                 aTTAATGTCGATCACGGCTATGTCCATGATATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AATAAAATGGGCTTTGCTCTGGTGGTAGTAACCATCCTGTCTGGCATTGCTCGCGGTAAATATAcatag                                                                            >  1:614517/1‑145 (MQ=255)
                                                                 aTTAATCTCGATCACGGCCATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATCGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                            <  2:707783/149‑1 (MQ=255)
                                                                   tAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGc                                                                          <  1:197051/149‑1 (MQ=255)
                                                                     aTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                        >  2:595657/1‑149 (MQ=255)
                                                                         cGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGtt                                                                    <  1:595657/149‑1 (MQ=255)
                                                                          gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGttt                                                                   >  1:985707/1‑149 (MQ=255)
                                                                            tCACGGCTATGTCCATGAGATCGACAACTTTGAAATTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                                 <  2:69584/149‑1 (MQ=255)
                                                                             cACGGATATGTCCTTGAGATCGACTACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTCCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaa                                                                <  2:737172/149‑1 (MQ=255)
                                                                              ccGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                               <  1:42723/148‑1 (MQ=255)
                                                                              aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                               >  2:969628/1‑149 (MQ=255)
                                                                              aCGGATATGTCAATGAGAGCGACAACTGTGAATTTATAGACGTTGTTATTGAAGCAATTCGCGAGCT‑AAAAAAATGGGCTTTGCGATGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaca                                                               <  1:640360/149‑3 (MQ=255)
                                                                                tgatATGTCCATGAGATCGAAAAATTTGAATTTATAGACTGTGTTATTGACGCCATGCGCGAGCT‑AAAAAATTGGGCTTTGCGCTTGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                             <  1:1017519/146‑1 (MQ=255)
                                                                                             agaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgga                                                <  2:425032/149‑1 (MQ=255)
                                                                                                 cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                            >  1:840526/1‑149 (MQ=255)
                                                                                                  gACAACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                           <  1:300373/149‑1 (MQ=255)
                                                                                                  gACAACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTAACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                           <  2:910513/149‑1 (MQ=255)
                                                                                                             aTTTATCGACGGTGTTAATGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCCCAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGc                                >  1:512559/1‑149 (MQ=255)
                                                                                                                     aCGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                        >  1:12716/1‑149 (MQ=255)
                                                                                                                       ggTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGAt                      <  2:1014570/149‑1 (MQ=255)
                                                                                                                             aTTGACGCCATGCGCGAGCT‑AAAATAATGGGCTTTGCGCTGGTGGTAGTAACCCCCCAGTGTGGCATTGGTCGCGGTATATTTACCGAAGGTCAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGTTGTCGAt                >  1:418511/1‑149 (MQ=255)
                                                                                                                              ttGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc               <  1:463672/149‑1 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCt              <  2:109961/147‑1 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCt              >  1:109961/1‑147 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg            >  2:795029/1‑149 (MQ=255)
                                                                                                                                   gCCATGTGCGAGCT‑TAATAAATGGTCTTTTCGCTGGTGGTAGTAACCAACAAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAGCTGGAt          <  1:660421/149‑1 (MQ=255)
                                                                                                                                    ccATGCGCGAGCT‑AAAAAAATGGGCGTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg         <  2:88901/149‑1 (MQ=255)
                                                                                                                                     cATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCTGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg        <  1:474564/149‑1 (MQ=255)
                                                                                                                                        gcgcGAGCT‑AAAAAAATCGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTAGAAACGATGACCGAGGGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTAt     >  2:338913/1‑149 (MQ=255)
                                                                                                                                          gcgAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTATTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCt   <  1:6902/149‑1 (MQ=255)
                                                                                                                                           cgAGCT‑AAAAACCTGGGCTTGGCGCTGGTGGTAGTAACCAACCAGTCTGTCATTGCTCGCGGAAAATTTACCGAAGCACAGTTTGATACGCTGACCTAGTGCATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCta  <  2:189814/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTA  >  NC_000913/222813‑223101

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