Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,576,827 C→T 6.3% A265T (GCA→ACA)  yddB ← putative TonB‑dependent outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,576,8270CT 6.3% 206.7 / 3.5 63A265T (GCA→ACA) yddBputative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (2/2);  total (30/33)
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.98e-01

CCCATTTCGCGATCAGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAAT  >  NC_000913/1576679‑1576967
                                                                                                                                                    |                                                                                                                                            
cccATTTCGCGATCAGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGc                                                                                                                                              <  4:681387/149‑1 (MQ=255)
 ccgtttCGCGATCAGACTTCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAAGCAGGTAAATATGCTTAATGCAGTATCGATAACGTTTTTATACTGTGcc                                                                                                                                             >  4:52390/4‑149 (MQ=255)
 ccATTTCGCGATCAGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGcc                                                                                                                                             >  3:600307/1‑149 (MQ=255)
    tttCGCGATCAGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGa                                                                                                                                          <  4:582526/149‑1 (MQ=255)
        gcgATCAGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCgg                                                                                                                                      <  3:143979/149‑1 (MQ=255)
                aCTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATg                                                                                                                              >  3:571315/1‑149 (MQ=255)
                 cTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGTCCGACCGGCGACAATGc                                                                                                                             >  2:558887/1‑149 (MQ=255)
                 cTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGTCCGACCGGCGACAATGc                                                                                                                             >  4:58221/1‑149 (MQ=255)
                  tGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGcc                                                                                                                            <  1:361286/149‑1 (MQ=255)
                            tGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGaa                                                                                                                  <  4:231651/149‑1 (MQ=255)
                             gCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGaaa                                                                                                                 <  4:247087/149‑1 (MQ=255)
                             gCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGaaa                                                                                                                 <  2:18934/149‑1 (MQ=255)
                               tGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACa                                                                                                               >  2:324988/1‑149 (MQ=255)
                                 gTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATa                                                                                                             <  2:423584/149‑1 (MQ=255)
                                    ttatAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACAGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATc                                                                                                          >  2:538932/1‑149 (MQ=255)
                                     tataATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATcc                                                                                                         >  3:203109/1‑149 (MQ=255)
                                          tCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCcgcgcg                                                                                                    >  3:155644/1‑149 (MQ=255)
                                          tCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCcgcgcg                                                                                                    >  2:590128/1‑149 (MQ=255)
                                            aCGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCggg                                                                                                  >  2:94669/1‑149 (MQ=255)
                                                    aGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCa                                                                                          <  4:203109/149‑1 (MQ=255)
                                                      ccGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGa                                                                                        <  3:52390/149‑1 (MQ=255)
                                                        ggTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACt                                                                                      <  1:35486/149‑1 (MQ=255)
                                                          tataTTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTGTCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCAt                                                                                                                       >  4:47121/1‑114 (MQ=255)
                                                          tataTTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCAt                                                                                                                       <  3:47121/114‑1 (MQ=255)
                                                              ttttAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTggcggc                                                                                >  2:115778/1‑149 (MQ=255)
                                                                  aaGGTTAAATCGTGGGTAAAGCGGCCGCTGGCAAACCAGGTAGATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGATCGGCGACAATGCCGGCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGa                                                                            <  4:600307/149‑1 (MQ=255)
                                                                     gTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTg                                                                                                >  3:303015/1‑126 (MQ=255)
                                                                     gTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTg                                                                                                <  4:303015/126‑1 (MQ=255)
                                                                     gTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATaa                                                                         >  1:332365/1‑149 (MQ=255)
                                                                        aaaTCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGc                                                                                                                             >  1:287404/1‑94 (MQ=255)
                                                                        aaaTCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGc                                                                                                                             <  2:287404/94‑1 (MQ=255)
                                                                           tCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACcggc                                                                   >  2:689287/1‑149 (MQ=255)
                                                                                ggTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTaa                                                              <  1:324988/149‑1 (MQ=255)
                                                                                   aaaGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGc                                                           >  4:82415/1‑149 (MQ=255)
                                                                                     aGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCa                                                                                          >  4:254763/1‑116 (MQ=255)
                                                                                     aGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCa                                                                                          <  3:254763/116‑1 (MQ=255)
                                                                                       cGGTTGCTAGCAAACCAGGTAAATTTGCTTAATGCCGTATCGATAACGTTTTTATACTGTGCCCGACCGTCGACAATTCCGTCATTCGAAACATAATCCGCGCGGGTGATACCAGACTGGCGGCGCGATAACCCGGCGTTAccgccggt                                                       <  1:538932/149‑3 (MQ=255)
                                                                                       cGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGaa                                                       <  3:125346/149‑1 (MQ=255)
                                                                                        ggTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGATAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAg                                                      <  1:365110/149‑1 (MQ=255)
                                                                                         gTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGt                                                     >  2:433602/1‑149 (MQ=255)
                                                                                          tCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGa                                                                                               <  2:31125/106‑1 (MQ=255)
                                                                                          tCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGa                                                                                               >  1:31125/1‑106 (MQ=255)
                                                                                          tCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGtt                                                    >  1:380130/1‑149 (MQ=255)
                                                                                                gCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGc                                              <  3:243444/149‑1 (MQ=255)
                                                                                                    aCCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAGCCGGCGGTAACGCCGAAGTTATCAGCGAGt                                          >  1:449854/1‑149 (MQ=255)
                                                                                                        ggTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCt                                      >  3:154339/1‑149 (MQ=255)
                                                                                                         gTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCa                                                >  4:13057/1‑138 (MQ=255)
                                                                                                         gTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCa                                                <  3:13057/138‑1 (MQ=255)
                                                                                                         gTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTg                                     >  1:656048/1‑149 (MQ=255)
                                                                                                          tAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGa                                    <  1:689287/149‑1 (MQ=255)
                                                                                                                tGCTTAATGCAGTATCGATAACGTTTTTATACTGTGTCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAAc                              <  3:58221/149‑1 (MQ=255)
                                                                                                                tGCTTAATGCAGTATCGATAACGTTTTTATACTGTGTCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAAc                              <  1:558887/149‑1 (MQ=255)
                                                                                                                 gCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACcggcg                                                                  >  1:196096/1‑112 (MQ=255)
                                                                                                                 gCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACcggcg                                                                  <  2:196096/112‑1 (MQ=255)
                                                                                                                     aaTGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGa                                    >  3:509310/1‑138 (MQ=255)
                                                                                                                     aaTGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGa                                    <  4:509310/138‑1 (MQ=255)
                                                                                                                         cAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGta                     >  4:519133/1‑149 (MQ=255)
                                                                                                                            tATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTc                                                                                                                     >  4:85345/1‑50 (MQ=255)
                                                                                                                            tATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTc                                                                                                                     <  3:85345/50‑1 (MQ=255)
                                                                                                                              tCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATaaaa                <  4:571315/149‑1 (MQ=255)
                                                                                                                                    aCGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCtt          <  3:519133/149‑1 (MQ=255)
                                                                                                                                     cGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCttt         <  2:380130/149‑1 (MQ=255)
                                                                                                                                           tataCTGTGCCCGACCGGCGACAATGCCGTCTTTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTaaaa   >  4:574591/1‑149 (MQ=255)
                                                                                                                                            ataCTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAAt  <  3:574591/149‑1 (MQ=255)
                                                                                                                                                tgtgCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCt                                      >  1:388052/1‑109 (MQ=255)
                                                                                                                                                tgtgCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTACt                                      <  2:388052/109‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CCCATTTCGCGATCAGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAAT  >  NC_000913/1576679‑1576967

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: ‑