Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,540,831 T→A 5.7% A283A (GCA→GCT cysW ← sulfate/thiosulfate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,540,8310TA 5.7% 206.3 / 3.3 70A283A (GCA→GCTcysWsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base A (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

CAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGC  >  NC_000913/2540683‑2540978
                                                                                                                                                    |                                                                                                                                                   
cagcGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCtgt                                                                                                                                                     >  1:84038/1‑149 (MQ=255)
            cccGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTc                                                                                                                                          <  8:47542/148‑1 (MQ=255)
             ccGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGTGACCAAACGACTTCGTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCt                                                                                                                                        >  6:166496/1‑149 (MQ=255)
             ccGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTcc                                                                                                                                         <  2:379681/148‑1 (MQ=255)
             ccGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCt                                                                                                                                        >  7:247867/1‑149 (MQ=255)
              cGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCACAATGCTAATGATGTTCATCCTCCTGTGCGCGTTTTTCCt                                                                                                                                        >  6:116034/1‑148 (MQ=255)
                       aCGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAg                                                                                                                              <  4:30815/149‑1 (MQ=255)
                         ggCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGgc                                                                                                                            <  2:84038/149‑1 (MQ=255)
                           ccccAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGgcgc                                                                                                                          >  1:106881/1‑149 (MQ=255)
                             ccAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGAtt                                                                                                                                    >  1:333197/1‑137 (MQ=255)
                             ccAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGAtt                                                                                                                                    <  2:333197/137‑1 (MQ=255)
                             ccAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACAAAACGACTTCTTAATATTGGCAATCTCAATTCTCACGAGGTTACTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCa                                                                                                                        >  1:339093/1‑149 (MQ=255)
                                   aCGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTcc                                                                                                                                         <  7:172535/126‑1 (MQ=255)
                                   aCGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTcc                                                                                                                                         >  8:172535/1‑126 (MQ=255)
                                    cgcgACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAAc                                                                                                                 <  8:187226/149‑1 (MQ=255)
                                      cgACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACAt                                                                                                               <  6:47752/149‑1 (MQ=255)
                                           aTCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGTGACCAAACGACTTCTTAATTTTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACttt                                                                                                          <  5:166496/149‑1 (MQ=255)
                                                aCCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTaa                                                                                                      <  8:107357/148‑1 (MQ=255)
                                                aCCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTa                                                                                                       <  7:306738/147‑1 (MQ=255)
                                                 ccTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTaaa                                                                                                     >  7:225907/1‑148 (MQ=255)
                                                     aaGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAg                                                                                                >  1:309532/1‑149 (MQ=255)
                                                      aGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAgg                                                                                               <  4:100404/149‑1 (MQ=255)
                                                       ggAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAggg                                                                                              <  7:211955/149‑1 (MQ=255)
                                                             tacAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATa                                                                                         >  4:38339/3‑148 (MQ=255)
                                                             tCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAgg                                                                                               >  5:724/1‑142 (MQ=255)
                                                             tCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAgg                                                                                               <  6:724/142‑1 (MQ=255)
                                                             tCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATaa                                                                                        >  1:243161/1‑149 (MQ=255)
                                                               cAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCa                                                                                                                        <  6:346492/115‑1 (MQ=255)
                                                                    agaTATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCAt                                                                                 >  6:311827/1‑149 (MQ=255)
                                                                       tatCGTACAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACa                                                                                                                <  1:50890/115‑1 (MQ=255)
                                                                       tatCGTACAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACa                                                                                                                >  2:50890/1‑115 (MQ=255)
                                                                         tCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCg                                                                            <  2:309532/149‑1 (MQ=255)
                                                                                 cACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAgc                                                                    >  4:255854/1‑149 (MQ=255)
                                                                                 cACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAgc                                                                    >  4:32089/1‑149 (MQ=255)
                                                                                    cTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGgc                                                                                                                            >  5:133890/1‑90 (MQ=255)
                                                                                    cTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGgc                                                                                                                            <  6:133890/90‑1 (MQ=255)
                                                                                               cAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTaaa                                                      <  8:225907/149‑1 (MQ=255)
                                                                                                 aaCGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAg                                                     >  7:68642/1‑148 (MQ=255)
                                                                                                 aaCGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAg                                                     >  4:186335/1‑148 (MQ=255)
                                                                                                  aCGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGa                                                   <  3:32089/149‑1 (MQ=255)
                                                                                                          ttAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACgg                                           >  3:241870/1‑149 (MQ=255)
                                                                                                            aaTCTTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACg                                            <  6:219494/146‑1 (MQ=255)
                                                                                                             atatTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGgtg                                         <  2:28602/148‑1 (MQ=255)
                                                                                                                ttGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGGTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGtgttg                                      >  8:49501/1‑148 (MQ=255)
                                                                                                                    cAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTc                                 <  8:68642/149‑1 (MQ=255)
                                                                                                                     aaTCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTcc                                >  4:188136/1‑149 (MQ=255)
                                                                                                                        ctcAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGgtg                                         <  5:38532/137‑1 (MQ=255)
                                                                                                                        ctcAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGgtg                                         >  6:38532/1‑137 (MQ=255)
                                                                                                                           aaTGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAg                                   >  8:39533/1‑140 (MQ=255)
                                                                                                                           aaTGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAg                                   <  7:39533/140‑1 (MQ=255)
                                                                                                                           aaTGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTc                           >  3:228307/1‑148 (MQ=255)
                                                                                                                           aaTGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCCCCACTGCAACATACTTTTTAAAAACAGGGGGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGc                             <  8:247867/146‑1 (MQ=255)
                                                                                                                                tCATGATGTTCCTCCTGCTGAGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTcc                                >  5:340929/1‑138 (MQ=255)
                                                                                                                                tCATGATGTTCCTCCTGCTGAGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTcc                                <  6:335459/138‑1 (MQ=255)
                                                                                                                                tCATGATGTTCCTCCTGCTGAGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTcc                                <  6:340929/138‑1 (MQ=255)
                                                                                                                                tCATGATGTTCCTCCTGCTGAGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTcc                                >  5:335459/1‑138 (MQ=255)
                                                                                                                                  atgatgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAAtt                   <  7:49501/149‑1 (MQ=255)
                                                                                                                                  atgatgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAAtt                   <  3:255854/149‑1 (MQ=255)
                                                                                                                                  atgatgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAAt                    <  3:186335/148‑1 (MQ=255)
                                                                                                                                   tgatgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAAtt                   >  7:220760/1‑148 (MQ=255)
                                                                                                                                    gatgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCa                         <  5:229321/141‑1 (MQ=255)
                                                                                                                                    gatgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCa                         >  6:229321/1‑141 (MQ=255)
                                                                                                                                    gatgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCa                         >  6:224147/1‑141 (MQ=255)
                                                                                                                                    gatgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCa                         <  5:224147/141‑1 (MQ=255)
                                                                                                                                     atgTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCACCAGGCGCCACTGCAAGATACTTTTTAAAAACAGGGTGATAATCGCCTTCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCa                 <  2:106881/148‑1 (MQ=255)
                                                                                                                                      tggtcctccTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAAt               <  4:79354/146‑1 (MQ=255)
                                                                                                                                           ctcctGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGTAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAAt               >  4:77338/1‑144 (MQ=255)
                                                                                                                                           ctcctGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGTAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAAt               <  3:77338/144‑1 (MQ=255)
                                                                                                                                           ctcctGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAAt               >  2:117609/1‑144 (MQ=255)
                                                                                                                                           ctcctGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAAt               <  1:117609/144‑1 (MQ=255)
                                                                                                                                              ctgctgTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGACATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTaa         <  2:339093/147‑1 (MQ=255)
                                                                                                                                                   gtgCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGc  >  4:148523/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTGAAGGAATATCCAGTGAGATATCGTTCAGCACCTGGGTGCGACCAAACGACTTCTTAATATTGGCAATCTCAATGCTCATGATGTTCCTCCTGCTGTGCGCGTTTTTCCTGATTCTCCAGGCGCCACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGC  >  NC_000913/2540683‑2540978

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