Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,539,929 C→T 5.4% V296I (GTA→ATA)  cysA ← sulfate/thiosulfate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,539,9290CT 5.4% 233.1 / 2.4 74V296I (GTA→ATA) cysAsulfate/thiosulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base T (2/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

AGCCGCGCATGTTGCAGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATG  >  NC_000913/2539781‑2540077
                                                                                                                                                    |                                                                                                                                                    
aGCCGCGCATGTTGCAGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTAc                                                                                                                                                      >  1:123372/1‑149 (MQ=255)
  ccGCGCATGTTGCAGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCt                                                                                                                                                    <  5:217216/149‑1 (MQ=255)
   cgcgCATGTTGCAGACCAACGAATAAGCGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTg                                                                                                                                                   <  5:76099/149‑1 (MQ=255)
     cgcATGTTGCAGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTa                                                                                                                                                 >  3:85933/1‑149 (MQ=255)
         tGTTGCAGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCgg                                                                                                                                             <  8:367555/149‑1 (MQ=255)
           ttGCAGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCgg                                                                                                                                             >  5:160201/1‑147 (MQ=255)
            tGCAGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAg                                                                                                                                          >  7:136969/1‑149 (MQ=255)
               aGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAgcgg                                                                                                                                       >  8:242585/1‑149 (MQ=255)
                gACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAgcggc                                                                                                                                      >  1:141863/1‑149 (MQ=255)
                 aCCAACGAATAAACGCTCGCTACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAgcggc                                                                                                                                      >  3:17539/1‑148 (MQ=255)
                    aaCGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAAt                                                                                                                                  >  6:110159/1‑149 (MQ=255)
                       gAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATc                                                                                                                                 >  6:98816/1‑147 (MQ=255)
                         aTAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAgg                                                                                                                             <  3:157210/149‑1 (MQ=255)
                         aTAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAgg                                                                                                                             >  8:182315/1‑149 (MQ=255)
                          tAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTAcc                                                                                                                                               >  2:169145/1‑130 (MQ=255)
                            aaCGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTg                                                                                                                          >  8:10352/1‑149 (MQ=255)
                                     cACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcggcgg                                                                                                                 <  7:242585/149‑1 (MQ=255)
                                             ggggCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATc                                                                                                         <  7:64754/149‑1 (MQ=255)
                                             ggggCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATc                                                                                                         <  6:160201/149‑1 (MQ=255)
                                             ggggCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATc                                                                                                         <  3:261110/149‑1 (MQ=255)
                                                   tcgtcgCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAg                                                                                                                              >  4:339565/1‑122 (MQ=255)
                                                   tcgtcgCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAg                                                                                                                              <  3:339565/122‑1 (MQ=255)
                                                     gtcgCCATGCATCACGACCGTCAGCGGTTCGTTGTACCCCCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTcc                                                                                                  >  5:135414/1‑148 (MQ=255)
                                                         ccATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAgg                                                                                             <  8:136969/149‑1 (MQ=255)
                                                          cATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCt                                                                                                               <  2:273346/130‑1 (MQ=255)
                                                                   cGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAag                                                                                   <  7:10352/149‑1 (MQ=255)
                                                                   cGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAag                                                                                   <  2:123372/149‑1 (MQ=255)
                                                                     aCCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAagag                                                                                 <  7:182315/149‑1 (MQ=255)
                                                                      ccGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAgaga                                                                                <  5:98816/149‑1 (MQ=255)
                                                                        gTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAag                                                                                   >  3:50918/1‑144 (MQ=255)
                                                                        gTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAag                                                                                   <  4:50918/144‑1 (MQ=255)
                                                                           aGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCAc                                                                           <  2:141863/149‑1 (MQ=255)
                                                                               gTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTAc                                                                                                                                                      <  6:59443/70‑1 (MQ=255)
                                                                               gTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTAc                                                                                                                                                      >  5:59443/1‑70 (MQ=255)
                                                                               gTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGaa                                                                                    >  3:177766/1‑136 (MQ=255)
                                                                               gTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGGGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGc                                                                       <  7:366628/149‑1 (MQ=255)
                                                                               gTTCGTGGTAGCACCCCAGCGGCTGGACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGATGCTGGTACGGCGGCTGATATCCACTTGCCAAGGGCGCAGGaa                                                                                    <  4:177766/136‑1 (MQ=255)
                                                                                       tACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATa                                                               >  1:333139/1‑149 (MQ=255)
                                                                                         ccaccCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACg                                                             >  6:446878/1‑149 (MQ=255)
                                                                                          caccCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAgcggc                                                                                                                                      <  1:14024/75‑1 (MQ=255)
                                                                                          caccCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAgcggc                                                                                                                                      >  2:14024/1‑75 (MQ=255)
                                                                                          caccCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACg                                                             >  3:402393/1‑148 (MQ=255)
                                                                                          caccCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACg                                                             >  4:191764/1‑148 (MQ=255)
                                                                                          caccCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACg                                                             >  4:161094/1‑148 (MQ=255)
                                                                                            ccccAGCGGCTGCACCAATAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCTGCGAATCGAGTCTGGTACGTCGGCTGATATACACTTCCCAAGGGCTCAGGAAGAGATCCACCGGCCCCTGATACTCAg                                                          <  4:165173/149‑1 (MQ=255)
                                                                                              gcAGCGGCTGCACCACTTATTGGGTGTAGTGACCTTTCGGGTTGGCTTCCTGTACCTGTACCGGCAGCGGCGAATCGAGTCTGGTACGGCTGCTGATATACACTTCCCAAGGGCGCAGGAAGAGATCCACCGGACCCTGATACGCAGgt                                                        <  4:308350/148‑1 (MQ=255)
                                                                                              ccAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGgt                                                        <  5:110159/149‑1 (MQ=255)
                                                                                              ccAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGgt                                                        >  8:282120/1‑149 (MQ=255)
                                                                                               cAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGATGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGgtg                                                       >  2:368336/1‑149 (MQ=255)
                                                                                                  cGGCTGCACCACTAATTGGGTGTAGTGACCTTAAGGGATGCCTTCCAGTACCTGTACCGCCAGCGGCGAATCGAGCCTGGAACGCCGGCTCATATCCACTTCCCAAGGCAGCAGAAAGAGAACCCCCGCCCCCTGATAAGCAGgtgtg                                                     >  6:211903/1‑148 (MQ=255)
                                                                                                        caccacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGcc                                                <  2:333139/147‑1 (MQ=255)
                                                                                                        caccacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCa                                              <  4:17539/149‑1 (MQ=255)
                                                                                                         accacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCg                                                                                                                                <  1:83628/66‑1 (MQ=255)
                                                                                                         accacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCg                                                                                                                                >  2:83628/1‑66 (MQ=255)
                                                                                                         accacTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTAACTGTACCGGAAGCGGCGAATAGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACTCAGGTGTGTAGCCCAg                                             >  1:160163/1‑149 (MQ=255)
                                                                                                            acTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCgg                                          >  3:429656/1‑149 (MQ=255)
                                                                                                              tAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGcc                                                <  7:77333/141‑1 (MQ=255)
                                                                                                              tAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGcc                                                >  8:77333/1‑141 (MQ=255)
                                                                                                               aaTTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTATCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAgaga                                                                                >  3:169368/1‑108 (MQ=255)
                                                                                                               aaTTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTATCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAgaga                                                                                >  3:175662/1‑108 (MQ=255)
                                                                                                               aaTTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTATCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAgaga                                                                                <  4:169368/108‑1 (MQ=255)
                                                                                                               aaTTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTATCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAgaga                                                                                <  4:175662/108‑1 (MQ=255)
                                                                                                                aTTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCa                                       <  4:85933/148‑1 (MQ=255)
                                                                                                                      tgtAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATgc                                >  4:217645/1‑149 (MQ=255)
                                                                                                                       gtAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATgcg                               >  8:21127/1‑149 (MQ=255)
                                                                                                                           tGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCt                                                                   <  7:196042/109‑1 (MQ=255)
                                                                                                                           tGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCt                                                                   >  8:196042/1‑109 (MQ=255)
                                                                                                                           tGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATgcgcgc                            >  2:407528/1‑148 (MQ=255)
                                                                                                                             aCCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCaa                         >  4:385492/1‑149 (MQ=255)
                                                                                                                                          gctTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGccc            >  4:207622/1‑149 (MQ=255)
                                                                                                                                            ttCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGACAGCGATGCGCGCCAACATGGAACTGCccgc          >  2:153257/1‑149 (MQ=255)
                                                                                                                                              ccAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACAt                      <  6:160154/135‑1 (MQ=255)
                                                                                                                                              ccAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACAt                      >  5:160154/1‑135 (MQ=255)
                                                                                                                                              ccAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCccgccg        >  2:351479/1‑149 (MQ=255)
                                                                                                                                                aGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCg                                                                         >  8:266219/1‑82 (MQ=255)
                                                                                                                                                aGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCg                                                                         <  7:266219/82‑1 (MQ=255)
                                                                                                                                                 gTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGa     <  7:282120/149‑1 (MQ=255)
                                                                                                                                                    ccTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATg  <  3:161094/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCCGCGCATGTTGCAGACCAACGAATAAACGCTCGCCACGCTGCGGGGCATCGTCGCCATGCATCACGACCGTCAGCGGTTCGTTGTACCACCCCAGCGGCTGCACCACTAATTGGGTGTAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGGCGGCTGATATCCACTTCCCAAGGGCGCAGGAAGAGATCCACCGGCCCCTGATACGCAGGTGTGTAGCCCAGCGGCCAGCGATGCGCGCCAACATGGAACTGCCCGCCGCGAATG  >  NC_000913/2539781‑2540077

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: