Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,625,304 Δ1 bp 6.8% intergenic (+13/+31) 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,3040G. 6.8% 250.9 / 11.1 74intergenic (+13/+31)ydeJ/dcpinactive PncC family protein/dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base . (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 = 4.57e-01

GTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGC  >  NC_000913/1625168‑1625450
                                                                                                                                        |                                                                                                                                                   
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)
                      cAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGctggcctgg                                                                                                                                                    <  6:99499/116‑1 (MQ=255)
                      cAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGctggcctgg                                                                                                                                                    >  5:99499/1‑116 (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                                                                                                                                     >  7:262521/1‑114 (MQ=255)
                                      ttATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGgcc                                                                                                                                     <  8:262521/114‑1 (MQ=255)
                                      ttATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggtg                                                                                                                          >  5:98203/1‑125 (MQ=255)
                                      ttATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggtg                                                                                                                          <  6:98203/125‑1 (MQ=255)
                                            aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTg                                                                                              >  5:277167/1‑147 (MQ=255)
                                            aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTg                                                                                              >  7:133034/1‑147 (MQ=255)
                                            aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTg                                                                                              >  5:195080/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)
                                                                                                                                        |                                                                                                                                                   
GTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGC  >  NC_000913/1625168‑1625450

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