Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,566,645 T→A 5.5% K201I (AAA→ATA) ‡ yjiN ← zinc‑type alcohol dehydrogenase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,566,6450TA 5.5% 315.0 / 7.1 90K201I (AAA→ATA) ‡yjiNzinc‑type alcohol dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (46/39);  new base A (3/2);  total (49/41)
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.54e-01

TAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAA  >  NC_000913/4566499‑4566792
                                                                                                                                                  |                                                                                                                                                   
tAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTa                                                                                                                                                   >  2:179874/1‑149 (MQ=255)
 aaTATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTAt                                                                                                                                                  >  2:47592/1‑149 (MQ=255)
  atatCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATc                                                                                                                                                 >  3:124616/1‑149 (MQ=255)
   tatCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATcg                                                                                                                                                >  1:487019/1‑149 (MQ=255)
        tcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaa                                                                                                                                           <  1:774398/149‑3 (MQ=255)
        tcaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATAAATCgcgcaa                                                                                                                                           <  1:633692/149‑3 (MQ=255)
          aaGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGga                                                                                                                                         <  2:481517/149‑1 (MQ=255)
            gCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGTGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGgaga                                                                                                                                       <  1:237582/149‑1 (MQ=255)
            gCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGgaga                                                                                                                                       >  2:742649/1‑149 (MQ=255)
             cAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCgcgc                                                                                                                                             <  3:303969/142‑1 (MQ=255)
             cAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCgcgc                                                                                                                                             >  4:303969/1‑142 (MQ=255)
             cAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGaa                                                                                                                                      >  4:548738/1‑149 (MQ=255)
              aaaGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAg                                                                                                                                     >  3:358354/1‑149 (MQ=255)
               aaGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAgg                                                                                                                                    >  1:39839/1‑149 (MQ=255)
                aGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAggg                                                                                                                                   >  4:315189/1‑149 (MQ=255)
                 gAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGc                                                                                                                                  <  1:391656/149‑1 (MQ=255)
                 gAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGc                                                                                                                                  >  1:10381/1‑149 (MQ=255)
                    ttCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGAt                                                                                                                               <  4:691216/149‑1 (MQ=255)
                       aCCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAgg                                                                                                                                    <  2:284552/141‑1 (MQ=255)
                       aCCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAgg                                                                                                                                    >  1:284552/1‑141 (MQ=255)
                       aCCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCaa                                                                                                                            >  1:280138/1‑149 (MQ=255)
                       aCCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCaa                                                                                                                            >  3:628947/1‑149 (MQ=255)
                       aCCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGCGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCaa                                                                                                                            <  1:227288/149‑1 (MQ=255)
                         cgcgTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACt                                                                                                                          >  4:198697/1‑149 (MQ=255)
                           cgTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtgt                                                                                                                        <  4:358354/149‑1 (MQ=255)
                             tCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGc                                                                                                                      >  4:22516/1‑149 (MQ=255)
                               aGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGa                                                                                                                    <  4:691019/149‑1 (MQ=255)
                               aGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGa                                                                                                                    >  2:743125/1‑149 (MQ=255)
                                gAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGAt                                                                                                                   >  1:73403/1‑149 (MQ=255)
                                gAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGAt                                                                                                                   >  4:40747/1‑149 (MQ=255)
                                 aaaCCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATc                                                                                                                  <  4:628947/149‑1 (MQ=255)
                                   aCCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCAAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAg                                                                                                                >  1:588411/1‑149 (MQ=255)
                                       aCTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTa                                                                                                            <  3:562358/149‑1 (MQ=255)
                                       aCTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTa                                                                                                            <  2:6117/149‑1 (MQ=255)
                                         tCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATc                                                                                                          <  2:39839/149‑1 (MQ=255)
                                         tCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATc                                                                                                          >  1:159223/1‑149 (MQ=255)
                                            gcgcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcag                                                                                                       <  1:550194/149‑1 (MQ=255)
                                            gcgcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGACAAGGGCGATCAACTGTGCGATCAGCGTATCcag                                                                                                       >  1:10684/1‑149 (MQ=255)
                                              gcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagca                                                                                                     <  2:369765/149‑1 (MQ=255)
                                                   gTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAgg                                                                                                                                    <  2:91533/113‑1 (MQ=255)
                                                   gTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAgg                                                                                                                                    >  1:91533/1‑113 (MQ=255)
                                                         ccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtgt                                                                                                                        <  4:144164/119‑1 (MQ=255)
                                                         ccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACtgt                                                                                                                        >  3:144164/1‑119 (MQ=255)
                                                         ccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatga                                                                                          <  1:540072/149‑1 (MQ=255)
                                                         ccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatga                                                                                          >  2:748389/1‑149 (MQ=255)
                                                         ccAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCtgatga                                                                                          >  3:209890/1‑149 (MQ=255)
                                                          cAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGAc                                                                                         >  4:627474/1‑149 (MQ=255)
                                                           aaCCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACg                                                                                        >  2:494010/1‑149 (MQ=255)
                                                            aCCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGa                                                                                       <  3:12018/149‑1 (MQ=255)
                                                               aTTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAgg                                                                                                                                    <  1:188742/101‑1 (MQ=255)
                                                               aTTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAgg                                                                                                                                    >  2:188742/1‑101 (MQ=255)
                                                                 tCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGtt                                                                                  <  1:645521/149‑1 (MQ=255)
                                                                  cGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGttt                                                                                 <  1:340603/149‑1 (MQ=255)
                                                                  cGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGttt                                                                                 <  3:627474/149‑1 (MQ=255)
                                                                     tGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTg                                                                              <  2:73403/149‑1 (MQ=255)
                                                                     tGGGGAGAATTTTGGCTTTCAGTGGAGGCTCGCTCTCCAGCCAGCGAACAAGTTGCTGGGCGATAAACTTGCGCGATTTATCGCGGTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATGCAGCAGCACCTGATGACGATCGTTTTTg                                                                              <  3:161500/149‑1 (MQ=255)
                                                                       gggAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGt                                                                            <  1:442939/149‑1 (MQ=255)
                                                                         gagaATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCa                                                                          >  4:304782/1‑149 (MQ=255)
                                                                          agaATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCTGCAGCACCTGATGACGATCGTTTTTGGTCAt                                                                         >  3:30840/1‑149 (MQ=255)
                                                                                     tttCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCa                                                              >  4:748256/1‑149 (MQ=255)
                                                                                     tttCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTATGGTCATACTCTCCAGCa                                                              >  4:172311/1‑149 (MQ=255)
                                                                                        cAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCa                                                           >  2:437293/1‑149 (MQ=255)
                                                                                         aGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCaa                                                          <  1:47592/149‑1 (MQ=255)
                                                                                          gTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAAc                                                         >  3:594537/1‑149 (MQ=255)
                                                                                           tGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACg                                                        >  4:543069/1‑149 (MQ=255)
                                                                                                gctcgctcTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTg                                                   <  3:132066/149‑1 (MQ=255)
                                                                                                  tcgctcTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGt                                                 >  4:530387/1‑149 (MQ=255)
                                                                                                   cgctcTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTg                                                >  1:673771/1‑149 (MQ=255)
                                                                                                        tcCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGa                                           >  4:322672/1‑149 (MQ=255)
                                                                                                                gcGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGa                                   <  4:549396/149‑1 (MQ=255)
                                                                                                                  gAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGAcc                                 >  4:498506/1‑149 (MQ=255)
                                                                                                                  gAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGAcc                                 >  1:284092/1‑149 (MQ=255)
                                                                                                                  gAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGAcc                                 >  3:695232/1‑149 (MQ=255)
                                                                                                                   aaCAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCACCGCACTGGTGCCGGAAAGATCGACCt                                <  2:656919/149‑1 (MQ=255)
                                                                                                                   aaCAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCt                                >  1:19039/1‑149 (MQ=255)
                                                                                                                     cAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTATGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTa                              <  3:172311/149‑1 (MQ=255)
                                                                                                                      aaTTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagca                                                                                                     >  4:461361/1‑77 (MQ=255)
                                                                                                                      aaTTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagca                                                                                                     <  3:461361/77‑1 (MQ=255)
                                                                                                                      aaTTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTAt                             <  3:748256/149‑1 (MQ=255)
                                                                                                                      aaTTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTAt                             <  3:22516/149‑1 (MQ=255)
                                                                                                                        tttGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCa                           >  2:93754/1‑149 (MQ=255)
                                                                                                                         ttGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCaa                          <  3:304782/149‑1 (MQ=255)
                                                                                                                            cTGGGCGATAAACTTGCACGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCg                       <  4:640180/149‑1 (MQ=255)
                                                                                                                                 cGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGa                  <  2:673771/149‑1 (MQ=255)
                                                                                                                                  gATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGAt                 <  3:315189/149‑1 (MQ=255)
                                                                                                                                        ccagcgcGATAAATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcg           >  3:321149/4‑149 (MQ=255)
                                                                                                                                        ccagcgcGATAAATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcg           >  1:672719/4‑149 (MQ=255)
                                                                                                                                        cTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGa                  >  4:526761/1‑142 (MQ=255)
                                                                                                                                        cTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGa                  <  3:526761/142‑1 (MQ=255)
                                                                                                                                          agcgcGATAAATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACcgcg         >  4:632871/2‑149 (MQ=255)
                                                                                                                                                 tttATCGCGCTGGAGGAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTaa  >  4:727106/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAA  >  NC_000913/4566499‑4566792

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: