Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,898,535 T→C 6.7% K295R (AAG→AGG)  ygcU ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,898,5350TC 6.7% 473.7 / 17.6 150K295R (AAG→AGG) ygcUputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (74/66);  new base C (3/7);  total (77/73)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

Only 100 of 155 total aligned reads displayed.

CGGTTTTGAGGATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGTT  >  NC_000913/2898389‑2898683
                                                                                                                                                  |                                                                                                                                                     
cGGTTTTGAGGATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTc                                                                                                                                                     >  2:617071/1‑149 (MQ=255)
   ttttGAGGATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCa                                                                                                                                                  >  2:309554/1‑149 (MQ=255)
       gAGGATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGc                                                                                                                                              >  1:405675/1‑149 (MQ=255)
         ggATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCTGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGa                                                                                                                                            <  2:760410/149‑1 (MQ=255)
         ggATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGa                                                                                                                                            >  3:545150/1‑149 (MQ=255)
         ggATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGa                                                                                                                                            <  3:526464/149‑1 (MQ=255)
          gATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAg                                                                                                                                           <  3:186499/149‑1 (MQ=255)
          gATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAg                                                                                                                                           <  3:927406/149‑1 (MQ=255)
          gATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAg                                                                                                                                           >  2:1537188/1‑149 (MQ=255)
            tCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCCTCGCATTGCGAagg                                                                                                                                         <  3:86544/149‑1 (MQ=255)
            tCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTAACCCTCGCATTGCGAagg                                                                                                                                         <  3:635353/149‑1 (MQ=255)
            tCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGACCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAggg                                                                                                                                         >  1:1488974/1‑149 (MQ=255)
             cTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCCTCGCATTGCGAAGGt                                                                                                                                        <  1:1041980/149‑1 (MQ=255)
              tGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGtt                                                                                                                                       <  3:1109486/149‑1 (MQ=255)
               gCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGTGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTa                                                                                                                                      >  2:112187/1‑149 (MQ=255)
                cacaCGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCCTCGCATTGCGAAGGt                                                                                                                                        <  4:1373675/146‑1 (MQ=255)
                cacaCGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCCTCGCATTGCGAAGGt                                                                                                                                        >  3:1373675/1‑146 (MQ=255)
                cacaCGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGACTTCCGCAATCCCTTCGCCCGTCACCCTCGCATTGCGAAGGTGAc                                                                                                                                     <  1:1577593/149‑1 (MQ=255)
                     gTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCa                                                                                                                                >  3:203070/1‑149 (MQ=255)
                     gTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCa                                                                                                                                >  2:1575491/1‑149 (MQ=255)
                     gTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCa                                                                                                                                <  2:1273308/149‑1 (MQ=255)
                        cGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGcc                                                                                                                             >  2:41836/1‑149 (MQ=255)
                             gCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCAACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCACGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATaaa                                                                                                                        <  4:836378/149‑1 (MQ=255)
                                aCTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGAt                                                                                                                     <  3:806750/149‑1 (MQ=255)
                                aCTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGAt                                                                                                                     <  2:1188174/149‑1 (MQ=255)
                                  tttATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                   <  4:821364/149‑1 (MQ=255)
                                   ttATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCCATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                   <  4:1489595/148‑1 (MQ=255)
                                   ttATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                   >  3:1489595/1‑148 (MQ=255)
                                      tCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATaaa                                                                                                                        <  3:649145/140‑1 (MQ=255)
                                      tCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATaaa                                                                                                                        >  4:649145/1‑140 (MQ=255)
                                      tCCGGTCCCAAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGc                                                                                                                              <  3:1147957/134‑1 (MQ=255)
                                      tCCGGTCCCAAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGc                                                                                                                              >  4:1147957/1‑134 (MQ=255)
                                             cccAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCCTCGCATTGCGAAGGt                                                                                                                                        >  2:1053754/1‑117 (MQ=255)
                                             ccaaGTTCAGGTTGTTAATCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACTTTCGCAATGTGAGGGTTACCTTCAGCCATAAAGATCAGCACGCAtttt                                                                                                        >  4:121145/4‑149 (MQ=255)
                                              ccAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTc                                                                                                      >  2:1424046/1‑149 (MQ=255)
                                              ccAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTc                                                                                                      >  3:1162106/1‑149 (MQ=255)
                                              ccAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTc                                                                                                      >  3:863228/1‑149 (MQ=255)
                                                aGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCg                                                                                                    <  1:1424046/149‑1 (MQ=255)
                                                 gTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAgg                                                                                                                                          <  4:415004/111‑1 (MQ=255)
                                                 gTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAgg                                                                                                                                          >  3:415004/1‑111 (MQ=255)
                                                 gTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGt                                                                                                   >  3:1119046/1‑149 (MQ=255)
                                                      ggTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCa                                                                                              >  2:257818/1‑149 (MQ=255)
                                                          gttAAACCAGGATTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAAt                                                                                          >  3:1013666/1‑149 (MQ=255)
                                                             aaaCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATggg                                                                                       <  1:257818/149‑1 (MQ=255)
                                                              aaCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCGCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGt                                                                                      <  2:420699/149‑1 (MQ=255)
                                                              aaCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACCCATTTTGCCGACGGAAAATGCGt                                                                                      <  2:1524306/149‑1 (MQ=255)
                                                               aCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTg                                                                                     >  1:306452/1‑149 (MQ=255)
                                                                ccAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGa                                                                                    >  1:1004202/1‑149 (MQ=255)
                                                                 cAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGaa                                                                                   <  4:275773/149‑1 (MQ=255)
                                                                 cAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGaa                                                                                   <  2:132454/149‑1 (MQ=255)
                                                                 cAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGaa                                                                                   <  4:1119046/149‑1 (MQ=255)
                                                                         gATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTggg                                                                           >  1:1081507/1‑149 (MQ=255)
                                                                          aTCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGt                                                                          <  2:193035/149‑1 (MQ=255)
                                                                          aTCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGt                                                                          >  1:527490/1‑149 (MQ=255)
                                                                           tCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                   <  2:67820/108‑1 (MQ=255)
                                                                           tCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                   >  1:67820/1‑108 (MQ=255)
                                                                           tCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTg                                                                         >  4:546382/1‑149 (MQ=255)
                                                                             aGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGcc                                                                       <  1:1428548/149‑1 (MQ=255)
                                                                               cTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCAt                                                                     >  1:765869/1‑149 (MQ=255)
                                                                                ttGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATc                                                                    <  3:1237173/149‑1 (MQ=255)
                                                                                     gTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAg                                                               >  1:905836/1‑149 (MQ=255)
                                                                                      tCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGc                                                              <  2:527490/149‑1 (MQ=255)
                                                                                        cACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAg                                                                                  <  2:1514337/127‑1 (MQ=255)
                                                                                         aCGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTc                                                           >  3:349742/1‑149 (MQ=255)
                                                                                          cgcgCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCa                                                          <  1:891362/149‑1 (MQ=255)
                                                                                             gcTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATAc                                                       >  1:60272/1‑149 (MQ=255)
                                                                                                ggCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTgg                                                                            >  3:1546374/1‑125 (MQ=255)
                                                                                                 gCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGaa                                                                                   <  1:1527253/117‑1 (MQ=255)
                                                                                                 gCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAAc                                                   >  3:55020/1‑149 (MQ=255)
                                                                                                   aTTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGa                                                 >  4:890757/1‑149 (MQ=255)
                                                                                                   aTTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGa                                                 <  2:158710/149‑1 (MQ=255)
                                                                                                    ttGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCa                                                          <  3:1416414/139‑1 (MQ=255)
                                                                                                    ttGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCa                                                          >  4:1416414/1‑139 (MQ=255)
                                                                                                    ttGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAg                                                >  1:448648/1‑149 (MQ=255)
                                                                                                       cGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGa                                             <  4:1537425/149‑1 (MQ=255)
                                                                                                         ggTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGAtt                                           <  2:537193/149‑1 (MQ=255)
                                                                                                            aGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgac                                        <  3:121145/149‑1 (MQ=255)
                                                                                                             gCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgacg                                       <  2:327731/149‑1 (MQ=255)
                                                                                                               gggCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgacgga                                     >  4:1281469/1‑149 (MQ=255)
                                                                                                                ggCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgacggac                                    <  4:726014/149‑1 (MQ=255)
                                                                                                                ggCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgacggac                                    >  4:137205/1‑149 (MQ=255)
                                                                                                                  cAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGt                                                                                                   <  1:199697/84‑1 (MQ=255)
                                                                                                                   aaCGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGAt                                 >  4:1474094/1‑149 (MQ=255)
                                                                                                                   aaCGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCACCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTTAAGTGTTGGGTGCCATCTTCAGCGACATACAAACGAGCGATTGACGGACGAt                                 >  2:649292/1‑149 (MQ=255)
                                                                                                                    aCGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATa                                >  2:797084/1‑149 (MQ=255)
                                                                                                                        tttCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCt                            >  1:1463294/1‑149 (MQ=255)
                                                                                                                              cAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTccacc                      <  3:1521607/149‑1 (MQ=255)
                                                                                                                                  cccTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGa                  <  2:121552/149‑1 (MQ=255)
                                                                                                                                    cTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTcca                        <  4:1400701/141‑1 (MQ=255)
                                                                                                                                     ttCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGAttt               >  1:824664/1‑149 (MQ=255)
                                                                                                                                           cGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCa         >  3:971376/1‑149 (MQ=255)
                                                                                                                                            gTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAg        <  1:112187/149‑1 (MQ=255)
                                                                                                                                            gTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAg        <  4:466902/149‑1 (MQ=255)
                                                                                                                                              cACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGa      <  4:1162106/149‑1 (MQ=255)
                                                                                                                                                 cTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGt   >  1:15271/1‑149 (MQ=255)
                                                                                                                                                  ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGaa                                                                                   >  1:949652/1‑68 (MQ=255)
                                                                                                                                                  ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGaa                                                                                   <  2:949652/68‑1 (MQ=255)
                                                                                                                                                  ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGtt  >  2:613883/1‑149 (MQ=255)
                                                                                                                                                  ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGtt  >  2:838657/1‑149 (MQ=255)
                                                                                                                                                  ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGtt  >  3:770712/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGTTTTGAGGATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTTCCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGTT  >  NC_000913/2898389‑2898683

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