Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,625,306 C→A 6.8% intergenic (+15/+29) 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,3060CA 6.8% 241.5 / 5.9 74intergenic (+15/+29)ydeJ/dcpinactive PncC family protein/dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base A (3/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

GTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGA  >  NC_000913/1625168‑1625452
                                                                                                                                          |                                                                                                                                                   
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                                                                                                                             <  1:317051/149‑1 (MQ=255)
             atTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggt                                                                                                                             >  3:316674/1‑149 (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                                                                                                                                  <  6:163284/137‑2 (MQ=255)
                    gCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGCCAGgccagcct                                                                                                                                  >  5:163284/1‑136 (MQ=255)
                      cAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTtat                                                                                                                    <  3:92610/149‑1 (MQ=255)
                        gAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTatat                                                                                                                  <  7:327112/149‑1 (MQ=255)
                         aaCTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTAtgtg                                                                                                                 >  1:89974/1‑149 (MQ=255)
                           cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGgtgt                                                                                                                        <  5:362703/140‑1 (MQ=255)
                           cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGgtgt                                                                                                                        >  6:362703/1‑140 (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                                                                                                >  5:277167/1‑147 (MQ=255)
                                            aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTg                                                                                                >  7:133034/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                                                                                       <  6:303148/78‑1 (MQ=255)
                                                                                                                          aaCCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCgg                                                                                       >  5:303148/1‑78 (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)
                                                                                                                                          |                                                                                                                                                   
GTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGA  >  NC_000913/1625168‑1625452

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: