Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,580,004 A→T 5.2% C263S (TGC→AGC)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,0040AT 5.2% 249.3 / 3.0 77C263S (TGC→AGC) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (39/34);  new base T (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCA  >  NC_000913/1579863‑1580137
                                                                                                                                             |                                                                                                                                     
gTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTcc                                                                                                                                >  8:264001/1‑149 (MQ=255)
 ttttttttGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCt                                                                                                                               <  2:10301/149‑1 (MQ=255)
    ttgttGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTgg                                                                                                                             <  8:403653/145‑1 (MQ=255)
          cgcTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATa                                                                                                                       <  6:416727/148‑1 (MQ=255)
          cgcTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATAc                                                                                                                      <  3:10645/149‑1 (MQ=255)
              gTCAGTACTACACTGATCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGAGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACt                                                                                                                  >  2:345718/1‑149 (MQ=255)
               tCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTTCAGGACTCCTGGAAAATACaaatt                                                                                                                 >  8:88265/1‑146 (MQ=255)
               tCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTg                                                                                                                 <  1:302278/149‑1 (MQ=255)
               tCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTg                                                                                                                 <  7:137530/149‑1 (MQ=255)
               tCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTg                                                                                                                 <  7:68732/149‑1 (MQ=255)
                  gTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATa                                                                                                                       <  7:181645/140‑1 (MQ=255)
                  gTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATa                                                                                                                       >  8:181645/1‑140 (MQ=255)
                   ttGTACACTTTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACTACTGGTAt                                                                                                             <  6:315612/149‑1 (MQ=255)
                      tACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGcc                                                                                                          >  2:65024/1‑149 (MQ=255)
                              tCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTccc                                                                                                   >  3:267205/1‑148 (MQ=255)
                              tCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGACGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTccc                                                                                                   >  7:72274/1‑148 (MQ=255)
                                 tCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCaa                                                                                                >  7:246037/1‑148 (MQ=255)
                                  cGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCAGTCCCCaaa                                                                                               >  4:218797/1‑148 (MQ=255)
                                    ttttGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAAc                                                                                           >  4:344326/1‑149 (MQ=255)
                                        gAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACgg                                                                                         <  1:198106/148‑1 (MQ=255)
                                         aGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCt                                                                                       <  6:236278/149‑1 (MQ=255)
                                         aGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCt                                                                                       <  6:243028/149‑1 (MQ=255)
                                           ttCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGa                                                                                     >  5:213840/1‑149 (MQ=255)
                                             cAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGAc                                                                                    >  7:309424/1‑148 (MQ=255)
                                                 tttAATAATGTTTCCAATTATGTACTCTGCATATACAAAATGGTCGCATTCGAAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGAATCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAg                                                                                <  7:349755/148‑1 (MQ=255)
                                                     ttaATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCttt                                                                            <  3:106126/148‑1 (MQ=255)
                                                             tCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGa                                                                    <  7:88265/148‑1 (MQ=255)
                                                              ccAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCt                                                                                                                               >  1:76344/1‑88 (MQ=255)
                                                              ccAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCt                                                                                                                               <  2:76344/88‑1 (MQ=255)
                                                                aaTTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGa                                                                >  2:70030/1‑149 (MQ=255)
                                                                       tACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGaaaa                                                          <  1:345718/148‑1 (MQ=255)
                                                                            tgGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCa                                                     <  7:264001/148‑1 (MQ=255)
                                                                            tgGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCAc                                                    <  1:65024/149‑1 (MQ=255)
                                                                               aTAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCACCAACCAGATTATCCTTGCAGTACTCCTGTAAAATACAACTGGTATGCCCATTCACCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCca                                                 >  3:209278/1‑149 (MQ=255)
                                                                                tAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCaacaa                                                 >  3:64083/1‑144 (MQ=255)
                                                                                   aCAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATc                                             >  4:218949/1‑149 (MQ=255)
                                                                                   aCAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATc                                             >  2:326017/1‑149 (MQ=255)
                                                                                       aaTGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGttt                                         >  1:118133/1‑149 (MQ=255)
                                                                                       aaTGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGttt                                         >  6:20001/1‑149 (MQ=255)
                                                                                        aTGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGtttt                                        <  7:183471/149‑1 (MQ=255)
                                                                                        aTGGACGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGtttt                                        <  8:72274/149‑1 (MQ=255)
                                                                                           gTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTa                                      <  8:246037/148‑1 (MQ=255)
                                                                                            tCGCATTCGTAAATGTCTCCATTACTTTCAAAAACCAGATTATCCTTGCAGTACTCCTGGAAAATACAACTGTTATGCGCATTCCCCAAAAAACGGCTGTCTAAGCTTTGAAAATGACTGATGAAAATTTCACCCACATCGTTTTTAAc                                    >  6:128553/1‑149 (MQ=255)
                                                                                              gCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATATCACCCACATCGTTTTTaaac                                   >  6:52808/1‑146 (MQ=255)
                                                                                               cATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCa                                 >  1:392128/1‑149 (MQ=255)
                                                                                                aTTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCAt                                <  2:200153/149‑1 (MQ=255)
                                                                                                aTTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTACCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCAt                                >  1:200153/1‑149 (MQ=255)
                                                                                                    gTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAACCGGCTGCCAAAGCTTACAAACTGACGGAGGAAAATTTCACCCACATCGTTTTTAACCCATTGCa                            >  7:367319/1‑149 (MQ=255)
                                                                                                    gTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCCCATCTTTTTTATCCCATTGCa                            >  2:252028/1‑149 (MQ=255)
                                                                                                      aaaTGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCAt                           >  4:232980/1‑148 (MQ=255)
                                                                                                        aTGTCTCCATTACTTTCAACAACCATATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGATTTCAAACTGACGTATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATaaa                        <  5:52808/149‑1 (MQ=255)
                                                                                                        aTGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATaa                         >  1:345541/1‑148 (MQ=255)
                                                                                                         tGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATaaa                        >  7:434782/1‑148 (MQ=255)
                                                                                                               cATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATg                    <  1:326017/146‑1 (MQ=255)
                                                                                                               cATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGtt                 <  6:213840/149‑1 (MQ=255)
                                                                                                               cATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGtt                 <  3:232980/149‑1 (MQ=255)
                                                                                                               cATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGt                  <  4:267205/148‑1 (MQ=255)
                                                                                                                aTTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGtt                                          <  6:214802/123‑1 (MQ=255)
                                                                                                                aTTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGtt                                          >  5:214802/1‑123 (MQ=255)
                                                                                                                aTTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTCACAAAGATATCAAACTGACGGAGGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATAGTTg                >  6:83132/1‑149 (MQ=255)
                                                                                                                    cTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCAtt                               <  6:35780/130‑1 (MQ=255)
                                                                                                                    cTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCAtt                               >  5:35780/1‑130 (MQ=255)
                                                                                                                     tttCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACa             >  4:356856/1‑147 (MQ=255)
                                                                                                                     tttCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACa             <  3:356856/147‑1 (MQ=255)
                                                                                                                     tttCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACAt            >  3:379259/1‑148 (MQ=255)
                                                                                                                     tttCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCAATTGCATAAAAATGGTTGACAt            >  3:162774/1‑148 (MQ=255)
                                                                                                                       tCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATaaa         >  3:140042/1‑149 (MQ=255)
                                                                                                                         aacaacCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc        <  3:382650/148‑1 (MQ=255)
                                                                                                                         aacaacCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt       >  3:280252/1‑149 (MQ=255)
                                                                                                                            aacCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGc    >  3:349197/1‑149 (MQ=255)
                                                                                                                              ccAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa  <  5:20001/149‑1 (MQ=255)
                                                                                                                               cAGATTATCCTTGCTGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg     >  2:18525/1‑145 (MQ=255)
                                                                                                                               cAGATTATCCTTGCTGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg     >  2:25159/1‑145 (MQ=255)
                                                                                                                               cAGATTATCCTTGCTGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg     <  1:25159/145‑1 (MQ=255)
                                                                                                                               cAGATTATCCTTGCTGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg     <  1:18525/145‑1 (MQ=255)
                                                                                                                                aGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACtt      >  1:365497/1‑143 (MQ=255)
                                                                                                                                aGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACtt      <  2:365497/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                     
GTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCA  >  NC_000913/1579863‑1580137

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