Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 620,790 G→A 5.0% V133V (GTC→GTT fepG ← iron‑enterobactin ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913620,7900GA 5.0% 259.0 / 2.3 80V133V (GTC→GTTfepGiron‑enterobactin ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (42/33);  new base A (2/2);  total (44/36)
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.67e-01

GCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAGGGCTGCCGAGCGGGTTACG  >  NC_000913/620644‑620936
                                                                                                                                                  |                                                                                                                                                  
gCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGAc                                                                                                                                                  <  1:230149/149‑1 (MQ=255)
   ttCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACcag                                                                                                                                               >  1:393045/1‑149 (MQ=255)
       acaAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCagcagc                                                                                                                                           >  1:263155/1‑149 (MQ=255)
           aCCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAg                                                                                                                                       <  2:138724/149‑1 (MQ=255)
            ccTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGt                                                                                                                                      <  1:21250/149‑1 (MQ=255)
             ctgctgTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACca                                                                                                                                                >  6:351237/1‑138 (MQ=255)
             ctgctgTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACca                                                                                                                                                <  5:351237/138‑1 (MQ=255)
                ctgTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGa                                                                                                                                    <  6:295178/147‑1 (MQ=255)
                     aGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACcagca                                                                                                                                             >  3:183735/1‑133 (MQ=255)
                     aGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACcagca                                                                                                                                             <  4:183735/133‑1 (MQ=255)
                        gcCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATgccgcc                                                                                                                          <  3:391269/149‑1 (MQ=255)
                         ccGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAACTGACAATgccgcc                                                                                                                          >  1:371065/1‑148 (MQ=255)
                            tttCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCAAGACCAGCAGCGAAGTGACAATGCCGCccac                                                                                                                       >  3:297111/1‑148 (MQ=255)
                                tAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGACCGCCGCTACCGATAATAATCAACCGAAAGGTGTCGCTGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCCCCATTg                                                                                                                  >  8:434601/1‑149 (MQ=255)
                                tAAAGACGCTTTCAGCAACAGCCAGGTACTAAAGGCCACCAGCATGGCGCGAACGCCGAAACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTg                                                                                                                  >  2:108208/1‑149 (MQ=255)
                                 aaaGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGc                                                                                                                 >  3:139665/1‑149 (MQ=255)
                                     aCGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCCCCATTGCTGAc                                                                                                             >  4:381035/1‑149 (MQ=255)
                                       gCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAg                                                                                                           >  3:161603/1‑149 (MQ=255)
                                              gCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGa                                                                                                      >  5:154361/1‑147 (MQ=255)
                                              gCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATa                                                                                                    <  1:347523/149‑1 (MQ=255)
                                              gCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATa                                                                                                    <  4:161603/149‑1 (MQ=255)
                                                 aCAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAAACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGcc                                                                                                 >  5:132638/1‑149 (MQ=255)
                                                 aCAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAAACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGcc                                                                                                 >  5:126904/1‑149 (MQ=255)
                                                    gCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGt                                                                                               >  8:303831/1‑148 (MQ=255)
                                                     cgaGGTATTAAAGGCCAGGAGCTTGGCGCGAACGCCGATACCGATTATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCa                                                                                             <  2:393045/147‑1 (MQ=255)
                                                          tATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGCCAGCGCGATAGCCGTCAGGTcc                                                                                        >  5:69521/1‑149 (MQ=255)
                                                          tATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTcc                                                                                        >  5:193577/1‑149 (MQ=255)
                                                             tAAAGGCCACCAGCATGGCGCGAACGCTGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGa                                                                                     >  3:324397/1‑149 (MQ=255)
                                                              aaaGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGa                                                                                     <  2:263155/148‑1 (MQ=255)
                                                              aaaGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGa                                                                                     <  3:19503/148‑1 (MQ=255)
                                                                aGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCa                                                                                  <  7:303831/149‑1 (MQ=255)
                                                                   ccaccaGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAAc                                                                               <  4:162159/149‑1 (MQ=255)
                                                                    caccaGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGATGCCAGCGCGATACCCGACAGGTCCTGACCaaccc                                                                              >  7:90885/1‑146 (MQ=255)
                                                                     accaGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAg                                                                             >  1:50756/1‑149 (MQ=255)
                                                                      ccaGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGc                                                                            <  4:139665/149‑1 (MQ=255)
                                                                         gCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCAc                                                                          >  8:133200/1‑148 (MQ=255)
                                                                               cgcgAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGcc                                                                   >  4:362557/1‑149 (MQ=255)
                                                                               cgcgAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAAACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGcc                                                                   <  6:126904/149‑1 (MQ=255)
                                                                               cgcgAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAAACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGcc                                                                   <  6:132638/149‑1 (MQ=255)
                                                                                      gCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCa                                                            >  7:191015/1‑149 (MQ=255)
                                                                                        cgataccgataATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCAcg                                                          >  7:12804/1‑149 (MQ=255)
                                                                                        cgataccgataATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCAcg                                                          >  2:138359/1‑149 (MQ=255)
                                                                                               gataATAATCAACCGAAAGGTGTCGCTGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTc                                                    >  5:47155/1‑148 (MQ=255)
                                                                                               gataATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTcc                                                   >  1:234621/1‑149 (MQ=255)
                                                                                                     aatCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGccc                                             >  1:123018/1‑149 (MQ=255)
                                                                                                      atCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCGCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGcccc                                            >  5:335606/1‑149 (MQ=255)
                                                                                                       tCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCg                                           <  1:108208/149‑1 (MQ=255)
                                                                                                       tCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCg                                           >  7:365471/1‑149 (MQ=255)
                                                                                                         aaCCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCt                                                                                       >  8:389753/1‑103 (MQ=255)
                                                                                                         aaCCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCt                                                                                       <  7:389753/103‑1 (MQ=255)
                                                                                                         aaCCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGgt                                         >  2:139817/1‑149 (MQ=255)
                                                                                                         aaCCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGgt                                         >  2:116407/1‑149 (MQ=255)
                                                                                                         aaCCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGgt                                         <  4:324397/149‑1 (MQ=255)
                                                                                                           ccGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGcca                                                                  >  4:158291/1‑122 (MQ=255)
                                                                                                           ccGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGcca                                                                  <  3:158291/122‑1 (MQ=255)
                                                                                                              aaaGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTaa                                    <  6:193577/149‑1 (MQ=255)
                                                                                                              aaaGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTa                                     >  5:380204/1‑148 (MQ=255)
                                                                                                              aaaGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTa                                     >  8:321483/1‑148 (MQ=255)
                                                                                                               aaGGTGTCGATGCCGTTGCGCCAGGCGCGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTa                                     >  8:341979/1‑147 (MQ=255)
                                                                                                               aaGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCCCCAGCACGCCGCTCCACGCCCCGGTGTTaaa                                   <  1:116407/149‑1 (MQ=255)
                                                                                                                  gtgtCGATGCCTTTGCGCCAGGCGAGCAGACATAACAGCAGCGAAGTGACAATGACGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGc                                 <  2:371065/148‑1 (MQ=255)
                                                                                                                  gtgtCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGcc                                >  2:254481/1‑149 (MQ=255)
                                                                                                                  gtgtCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGcc                                >  6:147351/1‑149 (MQ=255)
                                                                                                                   tgtCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGccc                               <  8:191015/149‑1 (MQ=255)
                                                                                                                    gtCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCa                              >  4:399554/1‑149 (MQ=255)
                                                                                                                     tCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATTGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCAt                             <  3:74037/149‑1 (MQ=255)
                                                                                                                     tCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCa                              <  3:167563/148‑1 (MQ=255)
                                                                                                                            cGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCTATTACGTCa                      <  2:50756/149‑1 (MQ=255)
                                                                                                                              ttGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATgcc                                                                                                                             >  6:253970/1‑44 (MQ=255)
                                                                                                                              ttGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATgcc                                                                                                                             <  5:253970/44‑1 (MQ=255)
                                                                                                                              ttGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAgg                    <  6:47155/149‑1 (MQ=255)
                                                                                                                              ttGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAgg                    <  2:234621/149‑1 (MQ=255)
                                                                                                                                         gAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAGGGCTGCCGAGCg         <  1:138359/149‑1 (MQ=255)
                                                                                                                                          agcagcCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAGGGCTGCCGAGCgg        <  3:362557/149‑1 (MQ=255)
                                                                                                                                          agcagcCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAGGGCTGCCGAGCgg        <  1:139817/149‑1 (MQ=255)
                                                                                                                                               ccAGTCCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGt                                                                                          >  3:145530/1‑62 (MQ=255)
                                                                                                                                               ccAGTCCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGt                                                                                          <  4:145530/62‑1 (MQ=255)
                                                                                                                                               ccAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAGGGCTGCCGAGCggg       >  8:15570/1‑145 (MQ=255)
                                                                                                                                               ccAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAGGGCTGCCGAGCggg       <  7:15570/145‑1 (MQ=255)
                                                                                                                                                cAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAGGGCTGCCGAGCGGGTTACg  >  6:112565/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCATTCCACAAACCTGCTGTTAGCGCCGTTTCTAAAGACGCTTTCAGCAACAGCCAGGTATTAAAGGCCACCAGCATGGCGCGAACGCCGATACCGATAATAATCAACCGAAAGGTGTCGATGCCGTTGCGCCAGGCGAGCAGCCAGACCAGCAGCGAAGTGACAATGCCGCCCACCATTGCTGACAGCGCGATAGCCGTCAGGTCCTGACCAAACAGCACCATCGCCACCAGCACGCCGCTCCACGCCCCGGTGTTAAAGCCCATTACGTCAGGGCTGCCGAGCGGGTTACG  >  NC_000913/620644‑620936

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