Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,512,350:1 +G 5.8% coding (65/957 nt) fecD ← ferric citrate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,3501.G 5.8% 584.5 / 32.1 138coding (65/957 nt)fecDferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (67/63);  new base G (6/2);  total (73/65)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

Only 100 of 150 total aligned reads displayed.

CAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGAC  >  NC_000913/4512076‑4512619
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                             
cggcggcgGAACAATCCCGGGCAGCGCAGGCAAAAGAAGTGGAGCCCCCGCAGAGGCCAGGCTGGCGGCATGGTTAACACCGGGAATAGCCGGTGATGCCAGAGGGTTGCGCACAATCCCCGGTATCAGCACGCCGGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCatat                                                                                                                                                                                                                                                                                <  4:479920/273‑1 (MQ=255)
                                        gggccagccACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtga                                                                                                                                                                                                                                        <  2:606562/267‑1 (MQ=255)
                                                     agGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAg                                                                                                                                                                                                                           >  2:376764/1‑274 (MQ=255)
                                                          aGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCCGCACCGACAAACAATGCCAGCAGCAAGCCCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCATTCGGTCAGCAGCGCGCGCCACGGCGCTGGGATCCCGCCCATAAGGGGTGATAACCGCGCACAGCCCGCCAGGCCAAGGGTGATGaaaagacacctgtccc                                                                                                                                                                                                                      >  4:606560/1‑262 (MQ=255)
                                                                                  agaATATCCGGTGATGCCAGAGGGTTCCGGCCAATCCCCTGTATCCGCTCTCCCGCGACTGCCAGGTCTTCCCCGTCTAACAGTTCCAGCATCACGCGCTGCAGACGGTACTCCATCAATGTATTATAGCGCTCGCGTTCGTCCTCCTCTTCCGCTAGCAGCGCGTGCAACGGCACAGGGCTCACTCCCATATGGAATGTTACCAGCTCAACACCCGCCAGGGCCAGGGTGATGAAAATCACCAGCGCCATTTTCATCACCCACTCCTCACaaac                                                                                                                                                                                             >  4:664084/1‑273 (MQ=255)
                                                                                      tatCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCatat                                                                                                                                                                                                                                                                                <  4:646640/189‑1 (MQ=255)
                                                                                       atCCGGCGTGGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                                                                                                                                                         <  2:80442/274‑1 (MQ=255)
                                                                                             tGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCCCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATCTGGAGGGATAACAGCGCACAACCCGCCAGGGCAAGGGTGCTGAAAATACCCAGCGCAAGTTTCCACCTCGCCCACTCACAAGccccacacacc                                                                                                                                                                                   >  4:477939/1‑265 (MQ=255)
                                                                                             agattcgacggggTGGCGCCAAACCCCCGGTATGAGCCGGCCCTCCACCGCGAGGGCCGAACCGACAAAATGTGCCAGGAGGAAGCGTGGCAGTGGCTACTCCAACAAAACATAATAATGCTCGGGTCCGGCCTGCAAGTGAGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCCGCGCAATTTTCATCCTCGCCTTCTCACAAGCCAGACAAAGCa                                                                                                                                                                                  <  4:376762/265‑1 (MQ=255)
                                                                                                                                gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATg                                                                                                                                                                                                                                                                              >  6:637730/1‑148 (MQ=255)
                                                                                                                                gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGc                                                                                                                                                                                             >  2:188702/1‑229 (MQ=255)
                                                                                                                                   cgcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCATATgggagt                                                                                                                                                                                                                                                                         <  7:328182/149‑5 (MQ=255)
                                                                                                                                    gcccgccACGGCGGGGGATGCACCGACAAACAGTGCCAGCAGAAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGGCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCatat                                                                                                                                                                                                                                                                                <  4:785646/143‑1 (MQ=255)
                                                                                                                                    gcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGTCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGt                                                                                                                                                                                                                                                                          <  5:665339/148‑1 (MQ=255)
                                                                                                                                    gcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                                                                                                                                                         <  4:188701/229‑1 (MQ=255)
                                                                                                                                        gccACGGCGAGGGCTGCACCGACAAACAGTGCCAGTAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTGAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATaa                                                                                                                                                                                                                                                                     <  7:616942/149‑1 (MQ=255)
                                                                                                                                         ccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAAc                                                                                                                                                                                                                                                                    <  7:160871/149‑1 (MQ=255)
                                                                                                                                          cACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACa                                                                                                                                                                                                                                                                   <  8:864290/149‑1 (MQ=255)
                                                                                                                                               cGAGGGCTGCACCGACAAACAGTGCAAGCAGAAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCAGATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGa                                                                                                                                <  2:532341/275‑1 (MQ=255)
                                                                                                                                                 aGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGcaca                                                                                                                                                                                                                                                            >  8:804840/1‑149 (MQ=255)
                                                                                                                                                         aCCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGcc                                                                                                                                                                                                                                                     <  5:152024/148‑1 (MQ=255)
                                                                                                                                                          ccGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGccc                                                                                                                                                                                                                                                        >  8:103540/1‑144 (MQ=255)
                                                                                                                                                           cGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAgg                                                                                                                                                                                                                                                  <  8:528398/149‑1 (MQ=255)
                                                                                                                                                             aCAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGc                                                                                                                                                                                                                                                >  7:841883/1‑149 (MQ=255)
                                                                                                                                                               aaaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                                                                                                                                              >  8:44086/1‑149 (MQ=255)
                                                                                                                                                                   aGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAggg                                                                                                                                                                                                                                           <  5:669527/148‑1 (MQ=255)
                                                                                                                                                                    gTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGACCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAgg                                                                                                                                                                                                                                            >  7:24703/1‑146 (MQ=255)
                                                                                                                                                                    gTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGACCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAgg                                                                                                                                                                                                                                            <  8:24703/146‑1 (MQ=255)
                                                                                                                                                                     tGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtg                                                                                                                                                                                                                                         <  5:112333/148‑1 (MQ=255)
                                                                                                                                                                     tGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtg                                                                                                                                                                                                                                         >  6:395962/1‑148 (MQ=255)
                                                                                                                                                                       ccAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTTCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGcc                                                                                                                                                                                                                                                     >  3:27732/1‑134 (MQ=255)
                                                                                                                                                                       ccAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaaa                                                                                                                                                                                                                                  >  1:206471/1‑153 (MQ=255)
                                                                                                                                                                            agcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAGCAGCGCACAGccc                                                                                                                                                                                                                                                        >  8:482867/1‑126 (MQ=255)
                                                                                                                                                                            agcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAGCAGCGCACAGccc                                                                                                                                                                                                                                                        <  7:482867/126‑1 (MQ=255)
                                                                                                                                                                               aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAGAATAAc                                                                                                                                                                                                                              <  8:737472/149‑1 (MQ=255)
                                                                                                                                                                               aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATaa                                                                                                                                                                                                                               >  6:310908/1‑148 (MQ=255)
                                                                                                                                                                                 gcgcGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAcc                                                                                                                                                                                                                             >  6:885925/1‑148 (MQ=255)
                                                                                                                                                                                    cgGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATc                                                                                                                                                                     <  2:29751/201‑1 (MQ=255)
                                                                                                                                                                                                    aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTc                                                                                                                                                                                                                 <  7:429402/141‑1 (MQ=255)
                                                                                                                                                                                                    aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTc                                                                                                                                                                                                                 >  8:429402/1‑141 (MQ=255)
                                                                                                                                                                                                    aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCg                                                                                                                                                                                                          >  6:1432/1‑148 (MQ=255)
                                                                                                                                                                                                    aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                                                                                                                                                         <  7:152741/149‑1 (MQ=255)
                                                                                                                                                                                                     tCAAGACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGCTGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                                                                                                                                                         >  5:259423/1‑148 (MQ=255)
                                                                                                                                                                                                      cAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCg                                                                                                                                                >  2:333111/1‑204 (MQ=255)
                                                                                                                                                                                                       aaTACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctc                                                                                                                                                                                                      <  7:188824/149‑1 (MQ=255)
                                                                                                                                                                                                       aaTACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcct                                                                                                                                                                                                       <  6:901706/148‑1 (MQ=255)
                                                                                                                                                                                                         tACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctcc                                                                                                                                                                                                     >  5:353266/1‑148 (MQ=255)
                                                                                                                                                                                                          aCATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctcct                                                                                                                                                                                                    <  6:477845/148‑1 (MQ=255)
                                                                                                                                                                                                          aCATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCa                                                                                                                                            <  4:333109/204‑1 (MQ=255)
                                                                                                                                                                                                             taataaTGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcaca                                                                                                                                                                                                >  7:767547/1‑149 (MQ=255)
                                                                                                                                                                                                               ataatGCTCGTGGCAGGCGTGCCCGTCGGTCAGAAGCGCGCGCCACGGCACGTGGGTCACTCCCGTATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGAGATGAAAATAACCAGCGCAATTTTCATCCTAGCCTCCTCACAAGCCAGAAAAAGCAAGGGCTGCCAAGCAGCGCCAGCACTGCGCATGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCa                                                                 <  4:146692/274‑1 (MQ=255)
                                                                                                                                                                                                                 aatGCTCGTGTCAGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGcc                                                                                                   <  1:153410/238‑1 (MQ=255)
                                                                                                                                                                                                                   tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAcc                                                                                                                                                                                                                             >  4:804853/1‑114 (MQ=255)
                                                                                                                                                                                                                     cTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                                                                                                                                                         <  6:353266/148‑1 (MQ=255)
                                                                                                                                                                                                                     cTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGAGATCACTCCCAGATGGAGTGATAACAGCGCACAGCCCGCCCAGGCAAGGCTGATGAACATAACCAGCGCAATTTTCATccgc                                                                                                                                                                                                           >  4:507856/1‑128 (MQ=255)
                                                                                                                                                                                                                      tCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                                                                                                                                                       >  8:578091/1‑149 (MQ=255)
                                                                                                                                                                                                                       cGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAgcg                                                                                                                                                                                                                         <  2:804855/114‑1 (MQ=255)
                                                                                                                                                                                                                       cGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                                                                                                                                                       >  5:606925/1‑148 (MQ=255)
                                                                                                                                                                                                                        gtgtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAgg                                                                                                                                                                                                                                                  >  1:48110/1‑88 (MQ=255)
                                                                                                                                                                                                                        gtgtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACaa                                                                                                                                                                                      <  5:310908/148‑1 (MQ=255)
                                                                                                                                                                                                                        gtgtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATg                                                                     >  2:146693/1‑261 (MQ=255)
                                                                                                                                                                                                                        gtgtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACAGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGCAGGCCAGCTCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCCGCAGCATGATCACTGGCAGTAc                                                       >  1:173964/1‑275 (MQ=255)
                                                                                                                                                                                                                          gtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGAATGGACCCGAGAAGCAGGCTACCTGGCAGTACg                                                      >  2:399796/1‑274 (MQ=255)
                                                                                                                                                                                                                            ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                                                                                                                                              <  3:48109/88‑1 (MQ=255)
                                                                                                                                                                                                                               gCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCACTCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTgcg                                                 >  2:456879/1‑274 (MQ=255)
                                                                                                                                                                                                                                 cTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGc                                                                                                                                                                             <  5:807803/148‑1 (MQ=255)
                                                                                                                                                                                                                                    ccAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGcc                                                                                                                                                                         >  8:88634/1‑149 (MQ=255)
                                                                                                                                                                                                                                         cGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCa                                                                                                                                                                    <  7:699604/149‑1 (MQ=255)
                                                                                                                                                                                                                                            tCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATa                                                                                                                                                                                                                                                                      <  5:685015/48‑1 (MQ=255)
                                                                                                                                                                                                                                            tCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATa                                                                                                                                                                                                                                                                      >  6:685015/1‑48 (MQ=255)
                                                                                                                                                                                                                                                   cgcgcgcCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAg                                                                                                                                                                   >  1:140402/1‑140 (MQ=255)
                                                                                                                                                                                                                                                    gcgcgcCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACt                                                                                                                                                         >  8:580748/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     cccgcCACGGCACGGGGATCACTCCCGTATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAg                                                                                        >  2:713303/3‑213 (MQ=255)
                                                                                                                                                                                                                                                              gcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCg                                                                                                                                                <  6:606925/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ggggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGaagct              >  2:481455/1‑270 (MQ=255)
                                                                                                                                                                                                                                                                    ttaTCACTCCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAt                                                                                                                                         >  7:719605/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATc                                                                                          >  2:389611/1‑196 (MQ=255)
                                                                                                                                                                                                                                                                    ggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGc            <  2:664086/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gagcACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCAAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTg                                                                                                                                                  >  5:219903/4‑139 (MQ=255)
                                                                                                                                                                                                                                                                     gATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTg                                                                                                                                                  <  6:219903/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tCACTCCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtctc                                                                                                                                      >  7:607796/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cactccATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtctc                                                                                                                                      >  5:312800/5‑148 (MQ=255)
                                                                                                                                                                                                                                                                         aCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGTTGAAAATAACCATCGCAATTTTCATCCTCGCCTCCTCACAAGCCATACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCg                                                                                                                                    >  8:664485/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGtt                                                    >  2:139346/1‑229 (MQ=255)
                                                                                                                                                                                                                                                                          cTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCGTCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                                                                                                                                                >  6:570133/1‑104 (MQ=255)
                                                                                                                                                                                                                                                                          actccATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTACCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCg                                                                                                                                    >  8:363618/4‑149 (MQ=255)
                                                                                                                                                                                                                                                                           tCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg                                                                                                                                  >  7:292678/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           tCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg                                                                                                                                  >  7:91580/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           ctccATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCg                                                                                                                                    >  6:790275/3‑148 (MQ=255)
                                                                                                                                                                                                                                                                           ctccATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCg                                                                                                                                    <  5:790275/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ccATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTgcgc                                                <  4:139345/229‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ccATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGa   <  2:477941/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ccATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCGCACTGGCAGTACGTTGCCCTGACCGAAGCCTGCCCAGAAGCCAGCCCGATGTGACACccgccacc  >  4:276338/1‑271 (MQ=255)
                                                                                                                                                                                                                                                                              cATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa                                                                                                                               <  8:767547/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               atatGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa                                                                                                                               <  5:1432/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               atatGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa                                                                                                                               <  5:395962/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                tatGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGc                                                                                                                                             >  1:46534/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                                tatGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGAt                 >  2:188594/1‑257 (MQ=255)
                                                                                                                                                                                                                                                                                 atGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAgg                                                                                                                             <  6:741639/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  tGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATc                                                                                                                                                                     <  4:778228/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                             
CAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGAC  >  NC_000913/4512076‑4512619

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