Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,057,859 T→C 6.0% K59R (AAA→AGA)  glnA ← glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,057,8590TC 6.0% 218.6 / 4.0 67K59R (AAA→AGA) glnAglutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base C (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

AATGGAGCGCGGGTCACGGTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGT  >  NC_000913/4057711‑4058003
                                                                                                                                                    |                                                                                                                                                
aaTGGAGCGCGGGTCACGGTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCtt                                                                                                                                                  >  4:470904/1‑149 (MQ=255)
 aTGGAGCGCGGGTCACGGTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCttt                                                                                                                                                 >  2:33644/1‑149 (MQ=255)
     aGCGCGGGTCACGGTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAg                                                                                                                                             <  3:470904/149‑1 (MQ=255)
                cGGTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCga                                                                                                                                  <  1:261067/149‑1 (MQ=255)
                 ggTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgag                                                                                                                                 >  3:278294/1‑149 (MQ=255)
                   tCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATTAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgaggc                                                                                                                               <  2:403458/149‑2 (MQ=255)
                      tAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGcc                                                                                                                            >  4:20557/1‑149 (MQ=255)
                          cTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCa                                                                                                                        <  3:102954/149‑1 (MQ=255)
                              cAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACa                                                                                                                    >  2:121973/1‑149 (MQ=255)
                               aGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACAt                                                                                                                   >  2:132471/1‑149 (MQ=255)
                                          ttCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGTGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttctt                                                                                                        >  3:16755/1‑149 (MQ=255)
                                              aGGATGTGGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaa                                                                                                    <  1:25687/149‑1 (MQ=255)
                                              aGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaa                                                                                                    <  4:459002/149‑1 (MQ=255)
                                              aGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaa                                                                                                    <  2:346283/149‑1 (MQ=255)
                                                  tGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAAt                                                                                                <  4:384029/149‑1 (MQ=255)
                                                   gTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAAtt                                                                                               <  1:121973/149‑1 (MQ=255)
                                                    tCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCCTCGAAGACTTc                                                                                              <  1:387742/149‑1 (MQ=255)
                                                     cGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCa                                                                                             <  2:382734/149‑1 (MQ=255)
                                                       cAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGc                                                                                           >  4:51901/1‑149 (MQ=255)
                                                        aaCGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGcc                                                                                                             >  1:170404/1‑130 (MQ=255)
                                                        aaCGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGcc                                                                                                             <  2:170404/130‑1 (MQ=255)
                                                           ggATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTc                                                                                       >  4:170146/1‑149 (MQ=255)
                                                               aaTCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCt                                                                                   <  1:63945/149‑1 (MQ=255)
                                                                aTCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTCTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCtg                                                                                  >  3:203047/1‑149 (MQ=255)
                                                                aTCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTCTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCtg                                                                                  >  1:346002/1‑149 (MQ=255)
                                                                 tCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCtga                                                                                 >  3:42169/1‑149 (MQ=255)
                                                                  cAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCtgat                                                                                <  3:20557/149‑1 (MQ=255)
                                                                  cAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCtgat                                                                                <  4:278294/149‑1 (MQ=255)
                                                                   aGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCtgatg                                                                               >  4:444309/1‑149 (MQ=255)
                                                                         gAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAgg                                                                         >  4:121915/1‑149 (MQ=255)
                                                                              ggCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAg                                                                    <  3:115885/149‑1 (MQ=255)
                                                                               gCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGt                                                                   >  3:387404/1‑149 (MQ=255)
                                                                                 gaagaaCGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgag                                                                                                                                 <  2:237711/85‑1 (MQ=255)
                                                                                 gaagaaCGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgag                                                                                                                                 >  1:237711/1‑85 (MQ=255)
                                                                                       cGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCt                                                           <  2:74852/149‑1 (MQ=255)
                                                                                        gggTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAggg                                                                        <  3:367583/135‑1 (MQ=255)
                                                                                        gggTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAggg                                                                        >  4:367583/1‑135 (MQ=255)
                                                                                           tCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTc                                                       <  2:301474/149‑1 (MQ=255)
                                                                                           tCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTc                                                       >  3:143130/1‑149 (MQ=255)
                                                                                             aaTCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCtt                                                     <  1:297098/149‑1 (MQ=255)
                                                                                                 aCTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTAcc                                                 >  3:211394/1‑149 (MQ=255)
                                                                                                     cGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTa                                             <  1:33644/149‑1 (MQ=255)
                                                                                                      ggTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAg                                            <  3:51901/149‑1 (MQ=255)
                                                                                                         ggATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTAt                                         >  4:218772/1‑149 (MQ=255)
                                                                                                                tCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGaa                                  <  2:72519/149‑1 (MQ=255)
                                                                                                                tCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTCTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGaa                                  <  2:346002/149‑1 (MQ=255)
                                                                                                                tCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTCTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGaa                                  <  4:203047/149‑1 (MQ=255)
                                                                                                                  tGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAgc                                <  3:121915/149‑1 (MQ=255)
                                                                                                                    gCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCa                                                                                             <  2:34583/86‑1 (MQ=255)
                                                                                                                    gCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCa                                                                                             >  1:34583/1‑86 (MQ=255)
                                                                                                                          gCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAgg                                                                         >  1:324425/1‑100 (MQ=255)
                                                                                                                          gCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAgg                                                                         <  2:324425/100‑1 (MQ=255)
                                                                                                                            aCCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATcaac                      >  2:33919/1‑149 (MQ=255)
                                                                                                                              cATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATcaacaa                    >  4:218738/1‑149 (MQ=255)
                                                                                                                                tGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGTCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACaaa                   <  3:1665/148‑1 (MQ=255)
                                                                                                                                tGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGTCTTCTTCGAAGAATCCAGCATTCACCTGATGAGCAGGGATAGTTCCGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACaaa                   >  4:1665/1‑148 (MQ=255)
                                                                                                                                    gggcTCGTTAATGCCGTTCCAGCCGCCAAGCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGCCGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAacttca              <  1:132471/146‑1 (MQ=255)
                                                                                                                                    ggACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTc                                                       >  2:117524/1‑108 (MQ=255)
                                                                                                                                    ggACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTc                                                       <  1:117524/108‑1 (MQ=255)
                                                                                                                                    ggACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAacttca              <  1:412808/149‑1 (MQ=255)
                                                                                                                                        tCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAgcgc                              >  2:153867/1‑129 (MQ=255)
                                                                                                                                        tCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAgcgc                              <  1:153867/129‑1 (MQ=255)
                                                                                                                                        tCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAacttcacttc          >  2:337145/1‑149 (MQ=255)
                                                                                                                                          gTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAAc                                                                                                                     <  2:394641/40‑1 (MQ=255)
                                                                                                                                          gTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAAc                                                                                                                     >  1:394641/1‑40 (MQ=255)
                                                                                                                                          gTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCt                                                           >  2:57791/1‑98 (MQ=255)
                                                                                                                                          gTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCt                                                           <  1:57791/98‑1 (MQ=255)
                                                                                                                                              aTGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTc    >  2:197285/1‑149 (MQ=255)
                                                                                                                                                gCCTTTCCAGCCGCCAATCGAGGAGCAGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGt  >  2:422152/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AATGGAGCGCGGGTCACGGTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGT  >  NC_000913/4057711‑4058003

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