Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,240 A→T 6.8% N122I (AAT→ATT)  slyB → outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,2400AT 6.8% 208.5 / 5.3 59N122I (AAT→ATT) slyBouter membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (21/34);  new base T (2/2);  total (23/36)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

TTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAT  >  NC_000913/1720092‑1720388
                                                                                                                                                    |                                                                                                                                                    
ttCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTaa                                                                                                                                                      <  2:248509/149‑1 (MQ=255)
 tCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAAt                                                                                                                                                     <  3:193765/149‑1 (MQ=255)
  ccTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATa                                                                                                                                                    <  3:58/149‑1 (MQ=255)
  ccTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGCCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATa                                                                                                                                                    <  1:490524/149‑1 (MQ=255)
   cTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATAc                                                                                                                                                   <  4:120815/149‑1 (MQ=255)
        gAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACcatca                                                                                                                                              >  3:80575/1‑149 (MQ=255)
            aCTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtggt                                                                                                                                          <  4:80575/149‑1 (MQ=255)
              tGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtggtgg                                                                                                                                        >  3:467431/1‑149 (MQ=255)
                ttGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTa                                                                                                                                      >  4:576486/1‑149 (MQ=255)
                  ggtggCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCACGGTGGTACa                                                                                                                                    <  1:594286/149‑1 (MQ=255)
                      gCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGaaa                                                                                                                                <  1:379492/149‑1 (MQ=255)
                        ggAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACa                                                                                                                              >  2:644659/1‑149 (MQ=255)
                         gAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACaa                                                                                                                             <  4:98093/149‑1 (MQ=255)
                         gAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACaa                                                                                                                             <  1:414888/149‑1 (MQ=255)
                          aaCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            <  2:506405/149‑1 (MQ=255)
                          aaCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            <  2:67164/149‑1 (MQ=255)
                          aaCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            <  4:23584/149‑1 (MQ=255)
                              gggCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCaa                                                                                                                        >  4:409422/1‑149 (MQ=255)
                                gCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAACCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAaca                                                                                                                      <  1:480776/149‑1 (MQ=255)
                                gCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAaca                                                                                                                      <  2:290634/149‑1 (MQ=255)
                                    tctcTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCg                                                                                                                  <  3:37854/149‑1 (MQ=255)
                                       ctGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGttt                                                                                                               <  1:493596/149‑1 (MQ=255)
                                              cTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCg                                                                                                        >  1:670334/1‑149 (MQ=255)
                                              cTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCg                                                                                                        <  3:576486/149‑1 (MQ=255)
                                               tGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCgg                                                                                                       >  4:652307/1‑149 (MQ=255)
                                                 cagcagGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGc                                                                                                     <  1:623531/149‑1 (MQ=255)
                                                    cagGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCaa                                                                                                  >  1:94084/1‑149 (MQ=255)
                                                      ggCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACg                                                                                                >  1:14145/1‑149 (MQ=255)
                                                                     ggCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAg                                                                                 <  3:652307/149‑1 (MQ=255)
                                                                        gTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCaa                                                                              <  2:14145/149‑1 (MQ=255)
                                                                                tCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTc                                                                      <  2:375477/149‑1 (MQ=255)
                                                                                 cAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCa                                                                     >  3:400338/1‑149 (MQ=255)
                                                                                   gggCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAgg                                                                   <  2:670334/149‑1 (MQ=255)
                                                                                     gCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACgt                                                                                               <  3:627827/119‑1 (MQ=255)
                                                                                     gCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACgt                                                                                               >  4:627827/1‑119 (MQ=255)
                                                                                          cAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAcc                                                              >  3:674482/1‑147 (MQ=255)
                                                                                          cAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAcc                                                              <  4:674482/147‑1 (MQ=255)
                                                                                          cAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGt                                                            >  2:238971/1‑149 (MQ=255)
                                                                                            gagTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGttt                                                          <  1:644659/149‑1 (MQ=255)
                                                                                            gagTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGttt                                                          <  2:94084/149‑1 (MQ=255)
                                                                                                gCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTc                                                                                            >  3:192462/1‑111 (MQ=255)
                                                                                                gCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTc                                                                                            <  4:192462/111‑1 (MQ=255)
                                                                                                  aaTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAc                                                               >  1:693819/1‑138 (MQ=255)
                                                                                                  aaTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAc                                                               <  2:693819/138‑1 (MQ=255)
                                                                                                  aaTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCcgc                                                    >  2:289960/1‑149 (MQ=255)
                                                                                                       aCAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTaa                                               >  1:344823/1‑149 (MQ=255)
                                                                                                           aaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATaa                                           <  3:409422/149‑1 (MQ=255)
                                                                                                           aaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATaa                                           <  2:344823/149‑1 (MQ=255)
                                                                                                               cAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGttt                                       <  1:316585/149‑1 (MQ=255)
                                                                                                                 gggTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGc                                     <  3:403221/149‑1 (MQ=255)
                                                                                                                 gggTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGc                                     <  4:400338/149‑1 (MQ=255)
                                                                                                                           cTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAg                           <  1:84869/149‑1 (MQ=255)
                                                                                                                              gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                        >  3:482485/1‑149 (MQ=255)
                                                                                                                                     gTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAcc                 >  2:169019/1‑149 (MQ=255)
                                                                                                                                        aaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcaca              >  3:390198/1‑149 (MQ=255)
                                                                                                                                         ccatgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCgtgt                                 <  2:504355/125‑1 (MQ=255)
                                                                                                                                         ccatgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCgtgt                                 >  1:504355/5‑129 (MQ=255)
                                                                                                                                          catgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                   <  2:466549/139‑1 (MQ=255)
                                                                                                                                          catgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                   >  1:466549/4‑142 (MQ=255)
                                                                                                                                           aCGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGTGTGTGGTCAGGTTACTGACCACACGc           >  1:102621/1‑149 (MQ=255)
                                                                                                                                                    aTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAt  >  1:627229/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAT  >  NC_000913/1720092‑1720388

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