Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,131,917 A→T 5.2% Y52N (TAT→AAT)  wcaC ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,131,9170AT 5.2% 255.1 / 3.3 77Y52N (TAT→AAT) wcaCputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base T (2/2);  total (38/39)
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.42e-01

AATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAAT  >  NC_000913/2131773‑2132064
                                                                                                                                                |                                                                                                                                                   
aaTGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTc                                                                                                                                                 <  4:79129/149‑1 (MQ=255)
 aTGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCt                                                                                                                                                <  6:139429/149‑1 (MQ=255)
    cAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAAATCATTGAAATTGCCAATCAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGa                                                                                                                                              <  7:377063/148‑1 (MQ=255)
     aGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATg                                                                                                                                            <  1:65671/149‑1 (MQ=255)
       gACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGc                                                                                                                                          >  8:301448/1‑149 (MQ=255)
         ccaccGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCt                                                                                                                                         >  8:326238/1‑148 (MQ=255)
          caccGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGa                                                                                                                                       >  5:306579/1‑149 (MQ=255)
                          cGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTgccgcc                                                                                                                       >  6:265550/1‑149 (MQ=255)
                            aGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCtt                                                                                                                     >  7:213364/1‑149 (MQ=255)
                             gTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCttt                                                                                                                    >  8:125236/1‑149 (MQ=255)
                              tAATGGTGCGATATAACTCATTGAAATTGCCAAAGAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCttt                                                                                                                    <  7:270836/148‑1 (MQ=255)
                              tAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGTATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTg                                                                                                                   <  7:372495/149‑1 (MQ=255)
                              tAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCCTTg                                                                                                                   <  7:123714/149‑1 (MQ=255)
                                aTGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGc                                                                                                                  <  8:213364/148‑1 (MQ=255)
                                aTGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGc                                                                                                                  <  7:186649/148‑1 (MQ=255)
                                   gTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGACACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATACTTCTGATGGCGGCCGATCTCTTTGCCGCCTTTGCCga                                                                                                               >  2:83127/1‑147 (MQ=255)
                                   gTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGt                                                                                                               >  1:134912/1‑148 (MQ=255)
                                      cGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAAcc                                                                                                           <  3:197766/149‑1 (MQ=255)
                                       gATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCg                                                                                                          >  4:290994/1‑149 (MQ=255)
                                        atatAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGt                                                                                                         <  5:143616/149‑1 (MQ=255)
                                         tataACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTa                                                                                                        <  7:209726/149‑1 (MQ=255)
                                          ataACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCg                                                                                                          >  4:149795/1‑146 (MQ=255)
                                           taACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGACTTTGCCGTAACCGTAg                                                                                                       >  7:144609/1‑148 (MQ=255)
                                              cTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACa                                                                                                    >  3:82288/1‑148 (MQ=255)
                                               tCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACaa                                                                                                   >  3:328257/1‑148 (MQ=255)
                                                 aTTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAAt                                                                                                <  2:111976/149‑1 (MQ=255)
                                                    gAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAAtgt                                                                                              <  2:134912/148‑1 (MQ=255)
                                                    gAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCAGTAACCGTAGACAAAAtgt                                                                                              >  1:262578/1‑148 (MQ=255)
                                                      aaTTGCCATACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTGTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGAc                                                                                           <  7:342342/149‑1 (MQ=255)
                                                             aaaCAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGc                                                                                     >  4:39221/1‑148 (MQ=255)
                                                             aaaCAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGc                                                                                     <  3:39221/148‑1 (MQ=255)
                                                              aaCAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGccc                                                                                   <  7:301448/149‑1 (MQ=255)
                                                               actgATCGCGATTTAACTGAGTACCCAGCGCAATATTCGCCTTCGCGGTCTTCGGAGGGGTATGTTTGATGCCCTGCGGATATTTCGGCTGGCTGAGGGTCTCTTTGCCGCGTTTTCAGTAACCGTAGACAAAATGTGACGCCAGcc                                                                                    <  4:283887/144‑1 (MQ=255)
                                                               aCAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGcccc                                                                                  <  2:45434/149‑1 (MQ=255)
                                                                     cgcgATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGc                                                                                                                  <  4:318432/111‑1 (MQ=255)
                                                                     cgcgATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGc                                                                                                                  >  3:318432/1‑111 (MQ=255)
                                                                            aaaCAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCa                                                                      <  7:326238/148‑1 (MQ=255)
                                                                             aaCAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATTGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACg                                                                    >  8:276197/1‑149 (MQ=255)
                                                                             aaCAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATTGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACg                                                                    >  8:276195/1‑149 (MQ=255)
                                                                             aaCAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGtt                                                                              <  5:361744/139‑1 (MQ=255)
                                                                             aaCAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGtt                                                                              >  6:361744/1‑139 (MQ=255)
                                                                             aaCAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACg                                                                    >  1:34045/1‑149 (MQ=255)
                                                                              aCAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGAAGCGCACGc                                                                   >  2:361044/1‑149 (MQ=255)
                                                                              aCAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTCATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGc                                                                   >  6:176656/1‑149 (MQ=255)
                                                                                  aCGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCtggt                                                               <  8:84030/149‑1 (MQ=255)
                                                                                      aaCAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCtggtg                                                              >  2:39382/1‑146 (MQ=255)
                                                                                       aCAGAGCAATATTCGCCATCGCGGGCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGgaga                                                          <  2:79837/149‑1 (MQ=255)
                                                                                          gagCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCt                                                       <  5:265550/149‑1 (MQ=255)
                                                                                           agCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAAc                                                                                                            >  6:225963/1‑95 (MQ=255)
                                                                                           agCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAAc                                                                                                            <  5:225963/95‑1 (MQ=255)
                                                                                             cAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCtggtg                                                              >  7:80232/1‑139 (MQ=255)
                                                                                             cAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCtggtg                                                              <  8:80232/139‑1 (MQ=255)
                                                                                                tatTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGcc                                                 >  5:240621/1‑149 (MQ=255)
                                                                                                     gCCATCGCGATCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACAcc                                            >  7:231001/1‑149 (MQ=255)
                                                                                                      ccATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCGCGCTGGTGTAGATCTAACGCCACACct                                           <  1:270029/149‑1 (MQ=255)
                                                                                                      ccATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACct                                           >  4:76057/1‑149 (MQ=255)
                                                                                                           gcgGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACctgctg                                       <  7:260835/148‑1 (MQ=255)
                                                                                                                  tCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGCTCTAACGCCACACCTGCTGCCCCGCCt                               >  7:127411/1‑149 (MQ=255)
                                                                                                                  tCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCt                               >  8:365198/1‑149 (MQ=255)
                                                                                                                  tCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCt                               >  3:157694/1‑149 (MQ=255)
                                                                                                                    cgcgGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTc                             >  5:318197/1‑149 (MQ=255)
                                                                                                                     gcggcgTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAAc                                                    <  7:145984/125‑1 (MQ=255)
                                                                                                                     gcggcgTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAAc                                                    >  8:145984/1‑125 (MQ=255)
                                                                                                                     gcggcgTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTcc                            >  2:215342/1‑149 (MQ=255)
                                                                                                                      cggcgTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCg                           >  6:342878/1‑149 (MQ=255)
                                                                                                                       ggcgTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGc                          <  7:35140/149‑1 (MQ=255)
                                                                                                                       ggcgTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGACAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGc                          <  5:141898/149‑1 (MQ=255)
                                                                                                                       ggcgTATGTTTGATGACCTGCGGATAGTGCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGc                          <  3:318949/149‑1 (MQ=255)
                                                                                                                         cgTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGc                          >  6:22358/1‑147 (MQ=255)
                                                                                                                             tGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTc                     <  6:306579/148‑1 (MQ=255)
                                                                                                                                      agcttcGGATAGTTCCGATGTCTGACGCTCTCTTTGCCGCCTTTGCCTTAACCGTATACAAAATGTGACGCCAGCCCCTGTGGCATCGCACGCTGGTGGAGATCTAAGGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTaaa           <  5:352041/144‑1 (MQ=255)
                                                                                                                                      aCCTGCGGATTGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTaaa           <  7:276195/149‑1 (MQ=255)
                                                                                                                                      aCCTGCGGATTGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTaaa           <  7:276197/149‑1 (MQ=255)
                                                                                                                                       ccTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAAt          <  3:368349/149‑1 (MQ=255)
                                                                                                                                       ccTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAAt          >  8:83073/1‑149 (MQ=255)
                                                                                                                                             gATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCaa     >  8:256178/1‑148 (MQ=255)
                                                                                                                                               tAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAat  <  3:290994/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAAT  >  NC_000913/2131773‑2132064

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: