Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,625,313 T→C 6.7% intergenic (+22/+22) ydeJ → / ← dcp inactive PncC family protein/dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,625,3130TC 6.7% 231.4 / 5.4 75intergenic (+22/+22)ydeJ/dcpinactive PncC family protein/dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base T (38/32);  new base C (3/2);  total (41/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 = 8.27e-01

GTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAG  >  NC_000913/1625168‑1625460
                                                                                                                                                  |                                                                                                                                                   
gTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATAtctc                                                                                                                                                    <  1:265195/147‑1 (MQ=255)
gTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCa                                                                                                                                                   <  8:107663/148‑1 (MQ=255)
gTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCa                                                                                                                                                   <  7:128359/148‑1 (MQ=255)
gTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAg                                                                                                                                                  >  5:34910/1‑149 (MQ=255)
 tttGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAgg                                                                                                                                                 >  6:130043/1‑149 (MQ=255)
       tGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATTTCTCAGgccagc                                                                                                                                            >  8:327733/1‑148 (MQ=255)
       tGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGgccagc                                                                                                                                            >  5:279902/1‑148 (MQ=255)
         gCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCAt                                                                                                                                         <  6:34910/149‑1 (MQ=255)
          cATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGTTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCAtt                                                                                                                                        <  1:81767/149‑1 (MQ=255)
           atatTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtg                                                                                                                                       >  4:92610/1‑149 (MQ=255)
            tatTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtgg                                                                                                                                      <  5:130043/149‑1 (MQ=255)
             atTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggt                                                                                                                                     >  3:316674/1‑149 (MQ=255)
             atTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggt                                                                                                                                     <  1:317051/149‑1 (MQ=255)
              ttAATGGCCAGAAGTTAACTGCGGTTATTCATTTTGCTGGGGAGTGCGATACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGGTGCTGCAATTTCTGGTATATCCAGGCTGGCCTGGCGATATCTGAGGCCAGCCATtggtg                                                                                                                                    <  4:316674/149‑1 (MQ=255)
              ttAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggtg                                                                                                                                    >  4:303629/1‑149 (MQ=255)
               tAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGCTATCTCAGGCCAGCCATtggtgg                                                                                                                                   >  6:239488/1‑149 (MQ=255)
               tAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGCCAGgccagcctggt                                                                                                                                       >  3:149574/1‑141 (MQ=255)
               tAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGCCAGgccagcctggt                                                                                                                                       <  4:149574/145‑5 (MQ=255)
                aaaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggtgg                                                                                                                                   <  8:323772/148‑1 (MQ=255)
                  agcccAGAACTACACTGCGGCTATGCATTTCGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTCAGGGCAGCCCCCGCCCAGCCGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGgtgt                                                                                                                                <  6:279902/146‑1 (MQ=255)
                    gCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGtt                                                                                                                               >  8:327112/1‑148 (MQ=255)
                    gCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGCCAGgccagcct                                                                                                                                          >  5:163284/1‑136 (MQ=255)
                    gCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGCCAGgccagcct                                                                                                                                          <  6:163284/137‑2 (MQ=255)
                      cAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTtat                                                                                                                            <  3:92610/149‑1 (MQ=255)
                        gAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTatat                                                                                                                          <  7:327112/149‑1 (MQ=255)
                         aaCTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTAtgtg                                                                                                                         >  1:89974/1‑149 (MQ=255)
                           cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGgtgt                                                                                                                                >  6:362703/1‑140 (MQ=255)
                           cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGgtgt                                                                                                                                <  5:362703/140‑1 (MQ=255)
                              cacTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                                      <  8:33000/147‑1 (MQ=255)
                                 tGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTGTGCCCTCGACCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCGGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGc                                                                                                                  >  5:253999/1‑148 (MQ=255)
                                      ttATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGgcc                                                                                                                                               <  8:262521/114‑1 (MQ=255)
                                      ttATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGgcc                                                                                                                                               >  7:262521/1‑114 (MQ=255)
                                      ttATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggtg                                                                                                                                    >  5:98203/1‑125 (MQ=255)
                                      ttATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggtg                                                                                                                                    <  6:98203/125‑1 (MQ=255)
                                            aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTg                                                                                                        >  5:195080/1‑147 (MQ=255)
                                            aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTg                                                                                                        >  7:133034/1‑147 (MQ=255)
                                            aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTg                                                                                                        >  5:277167/1‑147 (MQ=255)
                                                  cTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCg                                                                                                >  3:226455/1‑149 (MQ=255)
                                                   tGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCgg                                                                                               >  2:137960/1‑149 (MQ=255)
                                                        aCTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATaa                                                                                          >  6:275683/1‑149 (MQ=255)
                                                            cGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCtt                                                                                      <  4:192980/149‑1 (MQ=255)
                                                            cGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCAAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCTGCATAATCt                                                                                       >  6:333738/1‑148 (MQ=255)
                                                             gAAACGGTATTAGCTTTAGCGGTGAGGGTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGGTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTa                                                                                     <  2:60956/149‑1 (MQ=255)
                                                               aaCGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAgg                                                                                   <  4:226455/149‑1 (MQ=255)
                                                                   gTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGcc                                                                               >  3:289625/1‑149 (MQ=255)
                                                                      ttAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTa                                                                            >  6:151230/1‑149 (MQ=255)
                                                                       tAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTAc                                                                           >  2:388662/1‑149 (MQ=255)
                                                                        aGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTAcc                                                                          >  4:322687/1‑149 (MQ=255)
                                                                           tttAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACcgcg                                                                       >  6:368993/1‑149 (MQ=255)
                                                                             tAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCgcgc                                                                      <  5:151230/148‑1 (MQ=255)
                                                                                cGGGGAGGGTTGCCCTCGCCCAGCTGCTGCAATTACTGCGATAACCAGGCGGGCCGGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATTTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCAtt                                                                  <  6:253999/149‑1 (MQ=255)
                                                                                  gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGt                                                                >  1:14969/1‑149 (MQ=255)
                                                                                  gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAAACTTAGGTGCCTTACCGCGCCATTGt                                                                >  2:150706/1‑149 (MQ=255)
                                                                                      ggTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGAt                                                            <  5:333738/149‑1 (MQ=255)
                                                                                       gTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATa                                                           <  6:277167/149‑1 (MQ=255)
                                                                                          tGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATGTGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACa                                                         <  2:89974/148‑1 (MQ=255)
                                                                                            cccTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGc                                                      <  4:201977/149‑1 (MQ=255)
                                                                                             ccTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCg                                                     >  5:261018/1‑149 (MQ=255)
                                                                                              cTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGt                                                    <  8:133034/149‑1 (MQ=255)
                                                                                               tctcccAGCTTCTGCAATTAATTCTTTAACAAGGCTGGCCTGGCGATATCTCTTGACAGCCATTGGTTGTGTTTATATGTTCAAGGCACGATTTTGCAGCATCGGCATAATCTTATGTGCCTTTCCGCGCCATTGTCTATACAGGCGt                                                    <  6:321523/145‑1 (MQ=255)
                                                                                               tCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGtt                                                   >  5:326341/1‑149 (MQ=255)
                                                                                                  cccATCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCTCGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCa                                                <  7:154660/149‑1 (MQ=255)
                                                                                                   ccAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATTTCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAg                                               <  7:327733/149‑1 (MQ=255)
                                                                                                   ccAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTc                                                  >  3:142966/1‑146 (MQ=255)
                                                                                                   ccAGCTGCTGCAATTACTGCTATAACCAGGCTGGCATGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATCTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCTCCATAGTCGATACAGGCGt                                                    >  5:321523/1‑144 (MQ=255)
                                                                                                    cAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGa                                              <  6:195080/149‑1 (MQ=255)
                                                                                                      gctgctGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATc                                            <  3:322687/149‑1 (MQ=255)
                                                                                                                ttACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTa                                  >  7:94211/1‑149 (MQ=255)
                                                                                                                      cTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTg                            <  3:303629/149‑1 (MQ=255)
                                                                                                                        ataACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTTGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTgg                           >  6:320647/1‑148 (MQ=255)
                                                                                                                          aaCCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCgg                                                                                               >  5:303148/1‑78 (MQ=255)
                                                                                                                          aaCCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCgg                                                                                               <  6:303148/78‑1 (MQ=255)
                                                                                                                           caatggCTGGCCTGAGATATCGCCAGGCCAGCCTTTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAg                       >  8:173374/5‑149 (MQ=255)
                                                                                                                               ggCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATc                   >  4:79278/1‑149 (MQ=255)
                                                                                                                                       tggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTcgc            >  4:233807/1‑148 (MQ=255)
                                                                                                                                         gcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGa          >  3:317701/1‑148 (MQ=255)
                                                                                                                                                ctcAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAg  >  6:61164/1‑149 (MQ=255)
                                                                                                                                                  tcAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTcg             <  3:1276/137‑1 (MQ=255)
                                                                                                                                                  tcAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTcg             >  4:1276/1‑137 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAG  >  NC_000913/1625168‑1625460

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: