Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,794,441 C→T 6.0% A238T (GCG→ACG)  btuC ← vitamin B12 ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,794,4410CT 6.0% 214.7 / 2.7 67A238T (GCG→ACG) btuCvitamin B12 ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base T (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGATCAATGCCAGCATCAGCCAGCTTTGC  >  NC_000913/1794302‑1794586
                                                                                                                                           |                                                                                                                                                 
gCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCCCCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAg                                                                                                                                          >  6:62767/1‑149 (MQ=255)
gCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGccacca                                                                                                                                           <  3:117718/148‑1 (MQ=255)
gCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGccacca                                                                                                                                           <  7:362860/148‑1 (MQ=255)
gCCCCTGCCAGCGAGTAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACATCCGGAGAATATGTGGTATCACCAGATCAATATAGCCGATAGAAACCGGCAGCGCCACACTTAAGCCAACCATCCAGCCGGTCGCTGcctcct                                                                                                                                           <  6:216058/148‑2 (MQ=255)
  cccTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGca                                                                                                                                        >  5:158996/1‑149 (MQ=255)
  cccTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGca                                                                                                                                        >  2:166234/1‑149 (MQ=255)
        cAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTgc                                                                                                                                  <  2:303686/149‑1 (MQ=255)
                aGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAAc                                                                                                                          <  3:96470/149‑1 (MQ=255)
                aGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCAACCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAAc                                                                                                                          >  7:289363/1‑149 (MQ=255)
                      gAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGCCCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAgg                                                                                                                    >  1:221895/1‑149 (MQ=255)
                      gAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAgg                                                                                                                    <  1:192380/149‑1 (MQ=255)
                                  gatgatCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAggg                                                                                                                   >  5:247152/1‑138 (MQ=255)
                                  gatgatCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAggg                                                                                                                   <  6:247152/138‑1 (MQ=255)
                                  gatgatCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGt                                                                                                        >  1:291981/1‑149 (MQ=255)
                                      atCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGcc                                                                                                    <  6:158996/149‑1 (MQ=255)
                                       tCGGTTAAACCACACAACCGGAGAATATGGGGGATCCCCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCcg                                                                                                   >  1:184183/1‑149 (MQ=255)
                                         ggTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGt                                                                                                        >  8:121246/1‑142 (MQ=255)
                                         ggTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGt                                                                                                        <  7:121246/142‑1 (MQ=255)
                                           ttAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGcc                                                                                               <  6:166371/149‑1 (MQ=255)
                                                   cacaACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGAtctc                                                                                       <  2:212978/149‑1 (MQ=255)
                                                      aaCCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCa                                                                                                                              >  4:13242/1‑107 (MQ=255)
                                                      aaCCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCa                                                                                                                              <  3:13242/107‑1 (MQ=255)
                                                        ccGGAGAATATGGGGGATCACCAGCCCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGcc                                                                                    >  2:140537/1‑147 (MQ=255)
                                                           gagaATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCTCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCTCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTg                                                                               <  2:158095/149‑1 (MQ=255)
                                                            agaATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGc                                                                              <  5:62767/149‑1 (MQ=255)
                                                             gaATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGc                                                                              >  4:279456/1‑148 (MQ=255)
                                                               atatGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCTCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTTCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTaa                                                                           >  7:247377/1‑149 (MQ=255)
                                                                   gggggATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGAGCTAACata                                                                       >  5:76556/1‑149 (MQ=255)
                                                                    ggggATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACata                                                                       <  3:279456/148‑1 (MQ=255)
                                                                      tgATCACCAGACCAATAAAGCCGTTAGCACCCGCCAGCGCCACACTGACTCCAACCATCCAGCCGGTCGCTGCCACCAGGACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCATGTGCTAAGATATTc                                                                    <  8:247377/148‑1 (MQ=255)
                                                                      ggATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCCGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATAtt                                                                     >  2:122174/1‑148 (MQ=255)
                                                                           aCCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCGTCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGCACCACAGGGGTACACGCAGTTGCCGCGCCGAGATCTCGCCAAGTGGTAACATATTCATCgg                                                               >  4:79446/1‑149 (MQ=255)
                                                                             cAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGGCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGcc                                                             <  2:291981/149‑1 (MQ=255)
                                                                                aCCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGcc                                                                                    >  6:216980/1‑123 (MQ=255)
                                                                                aCCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGcc                                                                                    <  5:216980/123‑1 (MQ=255)
                                                                                 ccAATAAAGCCGATAGCATCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGa                                                         <  1:178187/149‑1 (MQ=255)
                                                                                  cAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGTTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGAc                                                        >  2:375522/1‑149 (MQ=255)
                                                                                  cAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGTTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGAc                                                        >  6:128073/1‑149 (MQ=255)
                                                                                    aTAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTg                                                      >  7:76648/1‑149 (MQ=255)
                                                                                     tAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGAtctc                                                                                       <  5:110456/115‑1 (MQ=255)
                                                                                     tAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGAtctc                                                                                       >  6:110456/1‑115 (MQ=255)
                                                                                     tAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGa                                                     <  1:166234/149‑1 (MQ=255)
                                                                                       aaGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCCCCAGCACATTGCGCCAGAACCACAGGGGTAACCCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGac                                                    >  1:48133/1‑148 (MQ=255)
                                                                                          ccGATAGCCCCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGCACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTACCCTATTCATCGGCCTGGACTGacaacc                                                >  5:158586/1‑148 (MQ=255)
                                                                                          ccGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCCGCACCACAGGGGTAAACGCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCCTCGGCCTGGACTGacca                                                  >  6:220671/1‑145 (MQ=255)
                                                                                          ccGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGacaa                                                  >  5:164209/1‑147 (MQ=255)
                                                                                          ccGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTCGACTGacaaca                                                <  8:289363/149‑1 (MQ=255)
                                                                                                 cACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCcac                                         <  2:184183/149‑1 (MQ=255)
                                                                                                   cccGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGTTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCacaa                                       <  5:128073/149‑1 (MQ=255)
                                                                                                   cccGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGTTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCacaa                                       <  1:375522/149‑1 (MQ=255)
                                                                                                   cccGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCacaa                                       <  1:140537/149‑1 (MQ=255)
                                                                                                       ccAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCacaacaa                                    >  3:58865/1‑148 (MQ=255)
                                                                                                            gcCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACggg                              <  8:76648/149‑1 (MQ=255)
                                                                                                            gcCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACggg                              >  7:258640/1‑149 (MQ=255)
                                                                                                             ccACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACggg                              >  8:92801/1‑148 (MQ=255)
                                                                                                              cacaCTGACGCCAACCATCCAGCCGGGCGCTGCCACCAGCACATTGCGCCAGAACCCCAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGa                            <  4:298148/149‑1 (MQ=255)
                                                                                                               acacTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCATAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGAGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGAt                           <  6:76556/149‑1 (MQ=255)
                                                                                                                cacTGACGCCAACCATCCAGCCGGTCGCTGCAACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGAt                           >  5:189582/1‑148 (MQ=255)
                                                                                                                 acTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGATc                          >  3:133209/1‑148 (MQ=255)
                                                                                                                     aCGCCAACCATCCAGCCGGTCGCTGCCACCAGCACCTTGCGCCAGAACCACAGGGGTAACCCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACCGGGCTCACTGc                     >  4:19590/1‑149 (MQ=255)
                                                                                                                     aCGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGCTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGATCAAt                       <  2:221895/147‑1 (MQ=255)
                                                                                                                          aaCGATCCAGCCGGTCGCTGCCACCAGGACATTGCGCCAGTACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGTCAACAGATCCACAACAACACGGGGATCAATGCCAGCa                <  4:133209/149‑1 (MQ=255)
                                                                                                                          aaCCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGATCAATGCCAGCa                <  6:164209/149‑1 (MQ=255)
                                                                                                                          aaCCATCAAGCCGGTCGCTGACACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGACTGGACTGACAACAGATCCACAACAACACGGGGATCAATGCCAGCa                <  6:158586/149‑1 (MQ=255)
                                                                                                                                       gTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGATCAATGc                     >  2:45564/1‑131 (MQ=255)
                                                                                                                                       gTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGATCAATGc                     <  1:45564/131‑1 (MQ=255)
                                                                                                                                        tCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGATCAATGCCAGCATCAGCCAGCTTTgc  <  7:92801/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GCCCCTGCCAGCGCGCAGCCGGGAAGTAATACGCGATGATCGGTTAAACCACACAACCGGAGAATATGGGGGATCACCAGACCAATAAAGCCGATAGCACCCGCCAGCGCCACACTGACGCCAACCATCCAGCCGGTCGCTGCCACCAGCACATTGCGCCAGAACCACAGGGGTAAACCCAGTTGCCGCGCCGAGATCTCGCCAAGTGCTAACATATTCATCGGCCTGGACTGACAACAGATCCACAACAACACGGGGATCAATGCCAGCATCAGCCAGCTTTGC  >  NC_000913/1794302‑1794586

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