Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,898,541 A→T 6.5% I293N (ATT→AAT)  ygcU ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,898,5410AT 6.5% 533.3 / 21.7 155I293N (ATT→AAT) ygcUputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (75/69);  new base T (3/7);  total (79/76)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01

Only 100 of 157 total aligned reads displayed.

GAGGATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGTTGAAGC  >  NC_000913/2898396‑2898688
                                                                                                                                                 |                                                                                                                                                    
  ggATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGa                                                                                                                                                 <  3:526464/149‑1 (MQ=255)
   gATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAg                                                                                                                                                <  3:186499/149‑1 (MQ=255)
   gATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAg                                                                                                                                                <  3:927406/149‑1 (MQ=255)
    aTCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAgg                                                                                                                                               <  3:819983/149‑1 (MQ=255)
    aTCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAgg                                                                                                                                               <  1:705542/149‑1 (MQ=255)
     tCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTAACCCTCGCATTGCGAagg                                                                                                                                              <  3:635353/149‑1 (MQ=255)
      cTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCCTCGCATTGCGAAGGt                                                                                                                                             <  1:1041980/149‑1 (MQ=255)
       tGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGtt                                                                                                                                            <  3:1109486/149‑1 (MQ=255)
       tGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGtt                                                                                                                                            <  3:975183/149‑1 (MQ=255)
       tGCACACGTTCGGCAGCCACTTTAACCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGtt                                                                                                                                            >  2:1428548/1‑149 (MQ=255)
        gCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGTGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTa                                                                                                                                           >  2:112187/1‑149 (MQ=255)
         cacaCGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCCTCGCATTGCGAAGGt                                                                                                                                             >  3:1373675/1‑146 (MQ=255)
            acGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCtt                                                                                                                                       >  4:1180778/1‑149 (MQ=255)
              gTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCa                                                                                                                                     >  4:1141212/1‑149 (MQ=255)
              gTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCa                                                                                                                                     >  2:1575491/1‑149 (MQ=255)
                      gCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCAACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCACGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATaaa                                                                                                                             <  4:836378/149‑1 (MQ=255)
                       ccACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAg                                                                                                                            >  3:549632/1‑149 (MQ=255)
                       ccACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGATGTCCACGCGCTCGCATTGCGGTTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAg                                                                                                                            >  2:157943/1‑149 (MQ=255)
                           tttATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                        <  4:821364/149‑1 (MQ=255)
                            ttATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCCATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                        <  4:1489595/148‑1 (MQ=255)
                            ttATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                        >  3:1489595/1‑148 (MQ=255)
                               tCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATaaa                                                                                                                             >  4:649145/1‑140 (MQ=255)
                               tCCGGTCCCAAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGc                                                                                                                                   <  3:1147957/134‑1 (MQ=255)
                               tCCGGTCCCAAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGc                                                                                                                                   >  4:1147957/1‑134 (MQ=255)
                                 cGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGc                                                                                                                  >  1:420699/1‑149 (MQ=255)
                                      ccaaGTTCAGGTTGTTAATCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACTTTCGCAATGTGAGGGTTACCTTCAGCCATAAAGATCAGCACGCAtttt                                                                                                             >  4:121145/4‑149 (MQ=255)
                                       ccAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑c                                                                                                           >  3:1162106/1‑149 (MQ=255)
                                       ccAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑c                                                                                                           >  3:863228/1‑149 (MQ=255)
                                       ccAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑c                                                                                                           >  2:1424046/1‑149 (MQ=255)
                                         aGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCg                                                                                                         <  1:1424046/149‑1 (MQ=255)
                                          gTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAgg                                                                                                                                               <  4:415004/111‑1 (MQ=255)
                                          gTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGt                                                                                                        >  3:1119046/1‑149 (MQ=255)
                                          gTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGt                                                                                                        >  3:275773/1‑149 (MQ=255)
                                               ggTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCa                                                                                                   >  2:257818/1‑149 (MQ=255)
                                                  tgttAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCaaaa                                                                                                <  4:1445357/149‑1 (MQ=255)
                                                   gttAAACCAGGATTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAAt                                                                                               >  3:1013666/1‑149 (MQ=255)
                                                    ttAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCGAAATg                                                                                              >  3:90051/1‑149 (MQ=255)
                                                      aaaCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATggg                                                                                            <  1:1575491/149‑1 (MQ=255)
                                                       aaCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCGCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGt                                                                                           <  2:420699/149‑1 (MQ=255)
                                                       aaCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCTAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGt                                                                                           <  2:708715/149‑1 (MQ=255)
                                                       aaCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACCCATTTTGCCGACGG‑AAAATGCGt                                                                                           <  2:1524306/149‑1 (MQ=255)
                                                        aCCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTg                                                                                          >  1:306452/1‑149 (MQ=255)
                                                         ccAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGa                                                                                         >  1:1004202/1‑149 (MQ=255)
                                                          cAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGaa                                                                                        <  2:132454/149‑1 (MQ=255)
                                                             gTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAgtg                                                                                     >  1:158710/1‑149 (MQ=255)
                                                              tttCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAgtgt                                                                                    >  1:587055/1‑149 (MQ=255)
                                                               ttCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGtt                                                                                   <  1:897739/149‑1 (MQ=255)
                                                                 cGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTgg                                                                                 <  4:139647/149‑1 (MQ=255)
                                                                 cGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTgg                                                                                 <  2:443616/149‑1 (MQ=255)
                                                                  gATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTggg                                                                                >  1:1081507/1‑149 (MQ=255)
                                                                   aTCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGt                                                                               <  2:629134/149‑1 (MQ=255)
                                                                   aTCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGt                                                                               <  2:193035/149‑1 (MQ=255)
                                                                   aTCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGt                                                                               >  1:527490/1‑149 (MQ=255)
                                                                    tCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                        <  2:67820/108‑1 (MQ=255)
                                                                    tCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCa                                                                                                                        >  1:67820/1‑108 (MQ=255)
                                                                    tCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTg                                                                              >  4:546382/1‑149 (MQ=255)
                                                                      aGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGcc                                                                            <  1:1428548/149‑1 (MQ=255)
                                                                             tGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGa                                                                                                                                                 >  4:1287714/1‑74 (MQ=255)
                                                                             tGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGa                                                                                                                                                 <  3:1287714/74‑1 (MQ=255)
                                                                                 cACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAg                                                                                       >  1:1514337/1‑127 (MQ=255)
                                                                                  aCGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCACTGCGAGGGTTACCTTCAGCCATAAAGATCAGCCCGCATTTT‑CCGTCGGCAAAATGGGAGAAGTGTTGGGTGCCGTCTTCAGCGTc                                                                >  2:1247020/1‑149 (MQ=255)
                                                                                  aCGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTc                                                                >  3:349742/1‑149 (MQ=255)
                                                                                         ggCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTgg                                                                                 >  3:1546374/1‑125 (MQ=255)
                                                                                         ggCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTgg                                                                                 <  4:1546374/125‑1 (MQ=255)
                                                                                          gCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGaa                                                                                        >  2:1527253/1‑117 (MQ=255)
                                                                                          gCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGaa                                                                                        <  1:1527253/117‑1 (MQ=255)
                                                                                          gCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAAc                                                        >  3:55020/1‑149 (MQ=255)
                                                                                            aTTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGa                                                      >  4:890757/1‑149 (MQ=255)
                                                                                            aTTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGa                                                      <  2:1158126/149‑1 (MQ=255)
                                                                                             ttGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCa                                                               <  3:1416414/139‑1 (MQ=255)
                                                                                             ttGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCa                                                               >  4:1416414/1‑139 (MQ=255)
                                                                                                cGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGa                                                  <  4:1537425/149‑1 (MQ=255)
                                                                                                 gggTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGAt                                                 >  2:376982/1‑149 (MQ=255)
                                                                                                  ggTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGAtt                                                <  2:537193/149‑1 (MQ=255)
                                                                                                     aGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgac                                             <  3:121145/149‑1 (MQ=255)
                                                                                                      gCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgacg                                            <  2:327731/149‑1 (MQ=255)
                                                                                                         ggCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgacggac                                         >  4:137205/1‑149 (MQ=255)
                                                                                                         ggCAACGATTTCCGCAATCCCCTCGCCCGCCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTgaccgac                                         <  2:1048007/149‑1 (MQ=255)
                                                                                                           cAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGt                                                                                                        <  1:199697/84‑1 (MQ=255)
                                                                                                            aaCGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGAt                                      >  4:1474094/1‑149 (MQ=255)
                                                                                                             aCGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATa                                     >  2:797084/1‑149 (MQ=255)
                                                                                                                 tttCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCt                                 >  1:1463294/1‑149 (MQ=255)
                                                                                                                      gCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTccac                            >  1:301613/1‑149 (MQ=255)
                                                                                                                      gCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTccac                            >  4:565627/1‑149 (MQ=255)
                                                                                                                       cAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTccacc                           <  3:1521607/149‑1 (MQ=255)
                                                                                                                           cccTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGa                       <  2:121552/149‑1 (MQ=255)
                                                                                                                             cTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGAtt                     >  3:755488/1‑149 (MQ=255)
                                                                                                                              ttCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGAttt                    >  1:824664/1‑149 (MQ=255)
                                                                                                                                     gTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAg             >  1:799107/1‑149 (MQ=255)
                                                                                                                                     gTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAg             <  4:466902/149‑1 (MQ=255)
                                                                                                                                     gTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAg             <  1:112187/149‑1 (MQ=255)
                                                                                                                                        aCCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGAt          >  3:41949/1‑149 (MQ=255)
                                                                                                                                           ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGaa                                                                                        >  1:949652/1‑68 (MQ=255)
                                                                                                                                           ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGaa                                                                                        <  2:949652/68‑1 (MQ=255)
                                                                                                                                           ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGtt       >  2:613883/1‑149 (MQ=255)
                                                                                                                                           ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGtt       >  2:838657/1‑149 (MQ=255)
                                                                                                                                           ttCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGtt       >  3:770712/1‑149 (MQ=255)
                                                                                                                                              gCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGTTGaa    <  2:800955/149‑1 (MQ=255)
                                                                                                                                                aaTGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGTTGAAGc  >  2:239366/1‑149 (MQ=255)
                                                                                                                                                aaTGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGTTGAAGc  >  1:1446016/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GAGGATCTGCACACGTTCGGCAGCCACTTTATCCGGTCCCCAGTTCAGGTTGTTAAACCAGGTTTCGATCAGCTTGCTGTCCACGCGCTGGCATTGCGGGTAGCGGGCAACGATTTCCGCAATCCCTTCGCCCGTCACCTTCGCAATGCGAGGGTTACCTTCAGCCATAAAGATCAGCACGCATTTT‑CCGTCGGCAAAATGGGTGAAGTGTTGGGTGCCATCTTCAGCGTCATACAAACGAGCGATTGACGGACGATACCCTTCCACCATGATTTCACGCAGGATGTTGAAGC  >  NC_000913/2898396‑2898688

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