Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,334,493 T→G 6.1% Y156S (TAC→TCC)  yfaT ← DUF1175 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,334,4930TG 6.1% 230.5 / 5.6 66Y156S (TAC→TCC) yfaTDUF1175 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base G (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCC  >  NC_000913/2334345‑2334640
                                                                                                                                                    |                                                                                                                                                   
cAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGt                                                                                                                                                     >  1:444804/1‑149 (MQ=255)
   aaaaTTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAAc                                                                                                                                                  >  1:410117/1‑149 (MQ=255)
            aCGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATcc                                                                                                                                         <  3:580965/149‑1 (MQ=255)
             cGATAAAGGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGATGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCa                                                                                                                                        >  4:434788/1‑149 (MQ=255)
              gATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAg                                                                                                                                       >  2:107168/1‑149 (MQ=255)
              gATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAg                                                                                                                                       >  2:547252/1‑149 (MQ=255)
               aTAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGa                                                                                                                                      <  3:203551/149‑1 (MQ=255)
                 aaaTGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGAcc                                                                                                                                    >  2:301142/1‑149 (MQ=255)
                         aTGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTg                                                                                                                            >  4:20368/1‑149 (MQ=255)
                           gAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCt                                                                                                                          <  3:6291/149‑1 (MQ=255)
                                         aaTCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGa                                                                                                            <  2:240648/149‑1 (MQ=255)
                                         aaTCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGa                                                                                                            >  3:135193/1‑149 (MQ=255)
                                         aaTCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGa                                                                                                            >  3:530235/1‑149 (MQ=255)
                                         aaTCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGa                                                                                                            >  4:156006/1‑149 (MQ=255)
                                           tCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATc                                                                                                          >  1:666712/1‑149 (MQ=255)
                                           tCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATc                                                                                                          >  2:371698/1‑149 (MQ=255)
                                                    aTCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAAt                                                                                                 <  1:335635/149‑1 (MQ=255)
                                                       cATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCat                                                                                              >  1:38561/1‑149 (MQ=255)
                                                         tCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCatat                                                                                            <  1:354210/149‑1 (MQ=255)
                                                          cGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATc                                                                                           <  3:585931/149‑1 (MQ=255)
                                                               gtCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAg                                                                                      >  2:394242/1‑149 (MQ=255)
                                                                tCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAgg                                                                                     >  4:647360/1‑149 (MQ=255)
                                                                  cTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCa                                                                                   >  3:156666/1‑149 (MQ=255)
                                                                  cTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCa                                                                                   <  4:562761/149‑1 (MQ=255)
                                                                   ttCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCAt                                                                                                                    >  1:152122/1‑115 (MQ=255)
                                                                   ttCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCAt                                                                                                                    <  2:152122/115‑1 (MQ=255)
                                                                   ttCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGcgc                                                                                  >  4:573587/1‑147 (MQ=255)
                                                                    tCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAgc                                                                                 <  1:646089/149‑1 (MQ=255)
                                                                         gTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTg                                                                            <  1:325680/149‑1 (MQ=255)
                                                                            aTAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGtt                                                                         <  1:107168/149‑1 (MQ=255)
                                                                            aTAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGtt                                                                         >  4:110230/1‑149 (MQ=255)
                                                                              aaGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCGCTGCCTggcgat                                                                       <  3:573587/149‑5 (MQ=255)
                                                                                gttgttGCAGAGTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATg                                                                     >  1:513273/1‑149 (MQ=255)
                                                                                    ttgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCa                                                                                               <  1:82255/119‑1 (MQ=255)
                                                                                    ttgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCa                                                                                               >  2:82255/1‑119 (MQ=255)
                                                                                    ttgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCt                                                                 >  3:274440/1‑149 (MQ=255)
                                                                                     tgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTg                                                                <  2:10905/149‑1 (MQ=255)
                                                                                         gactgactGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGttt                                                                        <  1:299966/137‑1 (MQ=255)
                                                                                         gactgactGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGttt                                                                        >  2:299966/1‑137 (MQ=255)
                                                                                                tgCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACt                                                     <  1:497302/149‑1 (MQ=255)
                                                                                                  cgcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTgg                                                   <  1:160369/149‑1 (MQ=255)
                                                                                                   gtgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGc                                                  <  4:407550/147‑1 (MQ=255)
                                                                                                     gcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCtg                                                >  2:213895/1‑149 (MQ=255)
                                                                                                        ttCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCtgttg                                             >  1:569666/1‑149 (MQ=255)
                                                                                                           cGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTAc                                          >  1:297892/1‑149 (MQ=255)
                                                                                                               gTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGaa                                      >  4:499508/1‑149 (MQ=255)
                                                                                                                tCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAAt                                     >  4:580195/1‑149 (MQ=255)
                                                                                                                  aGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCa                                   >  3:572973/1‑149 (MQ=255)
                                                                                                                    ttttCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCaaa                                 <  2:410117/149‑1 (MQ=255)
                                                                                                                    ttttCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCaaa                                 <  3:647360/149‑1 (MQ=255)
                                                                                                                    ttttCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCaaa                                 <  1:547252/149‑1 (MQ=255)
                                                                                                                              gcTTCCGGAGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCgg                       >  2:564700/1‑149 (MQ=255)
                                                                                                                                 tCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTg                                                    <  3:642058/117‑1 (MQ=255)
                                                                                                                                 tCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTg                                                    >  4:642058/1‑117 (MQ=255)
                                                                                                                                     gtgtGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATa                >  1:360508/1‑149 (MQ=255)
                                                                                                                                       gtgGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGAGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTCTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAgg              <  1:564700/149‑1 (MQ=255)
                                                                                                                                        tgGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAggg             >  3:611414/1‑149 (MQ=255)
                                                                                                                                          gTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGc           <  2:297892/149‑1 (MQ=255)
                                                                                                                                             gATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCgg        <  1:301142/149‑1 (MQ=255)
                                                                                                                                               tGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGtt      <  4:135193/149‑1 (MQ=255)
                                                                                                                                               tGACGGAACGCCCCATCCAGACCATTAAGTGCTGGTCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCaaa                                 >  3:306383/1‑122 (MQ=255)
                                                                                                                                               tGACGGAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCaaa                                 <  2:180781/122‑1 (MQ=255)
                                                                                                                                               tGACGGAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCaaa                                 <  4:306383/122‑1 (MQ=255)
                                                                                                                                               tGACGGAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCaaa                                 >  1:180781/1‑122 (MQ=255)
                                                                                                                                                  cGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTc   <  3:434788/149‑1 (MQ=255)
                                                                                                                                                  cGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTc   <  3:20368/149‑1 (MQ=255)
                                                                                                                                                   gTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTcc  <  4:530235/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCC  >  NC_000913/2334345‑2334640

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