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)
             ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGCTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgc                                                                                                                                >  1:826925/1‑149 (MQ=255)
             ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgc                                                                                                                                <  2:950205/149‑1 (MQ=255)
               aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTgcgc                                                                                                                              <  2:274512/149‑1 (MQ=255)
               aaaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACAGTGTGATTGACGCCATGCGCGAGCTAAAAAAATGGCCTTTgcgc                                                                                                                              >  1:607385/1‑149 (MQ=255)
                   ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGAACGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGATATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggt                                                                                                                          <  1:627257/149‑1 (MQ=255)
                     tGGCGAAGTCCGTACCCGCAATTTTTCTATACCGGGATTGCACCATAAATGTCGATCACGGCTATGACGATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGACCTAAAAAATTGGGCTTTGCGCtggtgg                                                                                                                        <  2:607385/149‑1 (MQ=255)
                          aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCCCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGGTATTGCCGACATGCGCGAGCTAAAAAATTGGGCTTTGCGCTGGTGgtagta                                                                                                                   >  2:1005307/1‑149 (MQ=255)
                           agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                  <  2:369146/149‑1 (MQ=255)
                               cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                              >  2:463672/1‑149 (MQ=255)
                                    ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGt                                                                                                         <  1:24162/149‑1 (MQ=255)
                                     cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                        >  1:425032/1‑149 (MQ=255)
                                       cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                      <  2:2901/149‑1 (MQ=255)
                                           ttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAt                                                                                                  >  1:1014570/1‑149 (MQ=255)
                                                 tGACCGTGATGGCACCATTAATGTCGAGCACGGCTATGTCCAAGAGATCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCACCCAGGCTGGCATTGCTcg                                                                                            >  1:583196/1‑149 (MQ=255)
                                                    ccGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTGGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGAGTTGCGCTGGTGGTATTAACCAACCAGTCTTGCATTGCTCGCgg                                                                                         >  2:601795/1‑149 (MQ=255)
                                                      gTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTa                                                                                       >  2:287022/1‑149 (MQ=255)
                                                            gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTa                                                                                 <  2:365903/149‑1 (MQ=255)
                                                            gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCGGGTGGTAGTAACCAACCAGTCTGGAATTGCTCGCGGTAAATTTa                                                                                 <  2:370059/149‑1 (MQ=255)
                                                             cACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                <  1:640768/149‑1 (MQ=255)
                                                              aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATAGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTATCCAACCAGTCTGGCATAGATCGCGGTAAATTTAcc                                                                               >  1:167084/1‑149 (MQ=255)
                                                                cATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGCCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGa                                                                             >  2:212279/1‑149 (MQ=255)
                                                                 aTTAATGTCGATCACGGCTATGTCCATGATATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAATAAAATGGGCTTTGCTCTGGTGGTAGTAACCATCCTGTCTGGCATTGCTCGCGGTAAATATAcatag                                                                            >  1:614517/1‑145 (MQ=255)
                                                                 aTTAATCTCGATCACGGCCATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATCGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                            <  2:707783/149‑1 (MQ=255)
                                                                   tAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGc                                                                          <  1:197051/149‑1 (MQ=255)
                                                                     aTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                        >  2:595657/1‑149 (MQ=255)
                                                                         cGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGtt                                                                    <  1:595657/149‑1 (MQ=255)
                                                                          gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGttt                                                                   >  1:985707/1‑149 (MQ=255)
                                                                            tCACGGCTATGTCCATGAGATCGACAACTTTGAAATTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                                 <  2:69584/149‑1 (MQ=255)
                                                                             cACGGATATGTCCTTGAGATCGACTACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTCCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaa                                                                <  2:737172/149‑1 (MQ=255)
                                                                              ccGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                               <  1:42723/148‑1 (MQ=255)
                                                                              aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                               >  2:969628/1‑149 (MQ=255)
                                                                              aCGGATATGTCAATGAGAGCGACAACTGTGAATTTATAGACGTTGTTATTGAAGCAATTCGCGAGCTAAAAAAATGGGCTTTGCGATGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaca                                                               <  1:640360/149‑3 (MQ=255)
                                                                                tgatATGTCCATGAGATCGAAAAATTTGAATTTATAGACTGTGTTATTGACGCCATGCGCGAGCTAAAAAATTGGGCTTTGCGCTTGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                             <  1:1017519/146‑1 (MQ=255)
                                                                                             agaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgga                                                <  2:425032/149‑1 (MQ=255)
                                                                                                 cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                            >  1:840526/1‑149 (MQ=255)
                                                                                                  gACAACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                           <  1:300373/149‑1 (MQ=255)
                                                                                                  gACAACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTAACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                           <  2:910513/149‑1 (MQ=255)
                                                                                                             aTTTATCGACGGTGTTAATGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCCCAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGc                                >  1:512559/1‑149 (MQ=255)
                                                                                                                     aCGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                        >  1:12716/1‑149 (MQ=255)
                                                                                                                       ggTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGAt                      <  2:1014570/149‑1 (MQ=255)
                                                                                                                             aTTGACGCCATGCGCGAGCTAAAATAATGGGCTTTGCGCTGGTGGTAGTAACCCCCCAGTGTGGCATTGGTCGCGGTATATTTACCGAAGGTCAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGTTGTCGAt                >  1:418511/1‑149 (MQ=255)
                                                                                                                              ttGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc               <  1:463672/149‑1 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCt              <  2:109961/147‑1 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCt              >  1:109961/1‑147 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg            >  2:795029/1‑149 (MQ=255)
                                                                                                                                   gCCATGTGCGAGCTTAATAAATGGTCTTTTCGCTGGTGGTAGTAACCAACAAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAGCTGGAt          <  1:660421/149‑1 (MQ=255)
                                                                                                                                    ccATGCGCGAGCTAAAAAAATGGGCGTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg         <  2:88901/149‑1 (MQ=255)
                                                                                                                                     cATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCTGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg        <  1:474564/149‑1 (MQ=255)
                                                                                                                                        gcgcGAGCTAAAAAAATCGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTAGAAACGATGACCGAGGGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTAt     >  2:338913/1‑149 (MQ=255)
                                                                                                                                          gcgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTATTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCt   <  1:6902/149‑1 (MQ=255)
                                                                                                                                           cgAGCTAAAAACCTGGGCTTGGCGCTGGTGGTAGTAACCAACCAGTCTGTCATTGCTCGCGGAAAATTTACCGAAGCACAGTTTGATACGCTGACCTAGTGCATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCta  <  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: