Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,621 A→C 7.9% A426A (GCT→GCG cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,890,6210AC 7.9% 203.3 / 7.6 73A426A (GCT→GCGcysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base C (2/4);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGAT  >  NC_000913/2890475‑2890766
                                                                                                                                                  |                                                                                                                                                 
gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGc                                                                                                                                                 <  3:171403/149‑1 (MQ=255)
  gcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGa                                                                                                                                               <  4:306120/149‑1 (MQ=255)
  gcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGa                                                                                                                                               <  3:131131/149‑1 (MQ=255)
    gTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAgg                                                                                                                                             <  4:296902/149‑1 (MQ=255)
    gTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAg                                                                                                                                              >  5:191015/1‑148 (MQ=255)
      tgctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGaa                                                                                                                                           <  7:155149/149‑1 (MQ=255)
       gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGaa                                                                                                                                           <  8:192905/148‑1 (MQ=255)
       gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGaa                                                                                                                                           <  2:168884/148‑1 (MQ=255)
        ctgcATATAGGCGCGGAACGGAGCAATACCGGTGCCTGGGGCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCTCCCTCTTCTTCCTCGCTGACCGCGTGGAAg                                                                                                                                          <  6:103252/148‑1 (MQ=255)
           cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTg                                                                                                                                       <  6:197247/148‑1 (MQ=255)
           cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTAGCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTg                                                                                                                                       <  8:307778/148‑1 (MQ=255)
            aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGaa                                                                                                                                           <  2:73414/143‑1 (MQ=255)
            aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGaa                                                                                                                                           >  1:73414/1‑143 (MQ=255)
             tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                    <  7:253547/149‑1 (MQ=255)
               aaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAgg                                                                                                                                   <  4:128428/148‑1 (MQ=255)
               aaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGGTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAgg                                                                                                                                   <  7:210048/148‑1 (MQ=255)
                 ggCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGcac                                                                                                                                <  8:10547/149‑1 (MQ=255)
                    gcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTACACGCGGTCAGCGAGGAAGCTGGAGGcaccac                                                                                                                             >  5:276881/1‑149 (MQ=255)
                         aCGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCGTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACGACCGGCa                                                                                                                        <  8:86169/149‑1 (MQ=255)
                         aCGGCGCAATACAGGTTCCTGTGCAAATCATAATCACCGGGGTTTCTGGATTGGCTTGCAGGCGGAAGTTATCGTTAAGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCa                                                                                                                        <  8:81318/149‑1 (MQ=255)
                                       gTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTc                                                                                                          >  3:107889/1‑149 (MQ=255)
                                                    tCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACg                                                                                              >  3:202965/1‑148 (MQ=255)
                                                     cATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGca                                                                                            >  4:239534/1‑149 (MQ=255)
                                                      aTAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGca                                                                                            <  1:75538/148‑1 (MQ=255)
                                                       tAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcac                                                                                           >  6:179730/1‑148 (MQ=255)
                                                                  gTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGt                                                                               <  6:191015/149‑1 (MQ=255)
                                                                      cTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGAc                                                                            <  1:75961/148‑1 (MQ=255)
                                                                       tGGATTGGCTGGCAGGCGGAAGTTAGCGGTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGATGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAGGCACCCCCCCAACGGTGACgt                                                                          >  6:265000/1‑149 (MQ=255)
                                                                         gATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACgtg                                                                         >  5:338183/1‑148 (MQ=255)
                                                                         gATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTACTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACgtg                                                                         >  8:218831/1‑148 (MQ=255)
                                                                              gctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTc                                                                   >  2:245876/1‑149 (MQ=255)
                                                                              gctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTc                                                                   >  3:232889/1‑149 (MQ=255)
                                                                              gctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGGCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACCCCAACGGGGACGTGTACTTc                                                                   <  7:243299/149‑1 (MQ=255)
                                                                                tggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGt                                                                 >  3:327164/1‑149 (MQ=255)
                                                                                  gcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtc                                                               >  8:248652/1‑149 (MQ=255)
                                                                                    aGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACATCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGTCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtctc                                                             <  4:4859/149‑1 (MQ=255)
                                                                                     ggCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCg                                                            <  8:132009/149‑1 (MQ=255)
                                                                                                tCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacac                                                                                      >  1:344916/1‑112 (MQ=255)
                                                                                                tCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacac                                                                                      <  2:344916/112‑1 (MQ=255)
                                                                                                   ttATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGACACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGc                                                >  4:290801/1‑147 (MQ=255)
                                                                                                     aTGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACAACCGGCACGGGCGCGGCCTTACACGTCGTAACGCACCACACAAACGTTGACGTGTAATTCGTTGTCGACTTCCGCCTGCgag                                             >  6:237636/1‑148 (MQ=255)
                                                                                                       gTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGGCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgagga                                           <  6:213270/148‑1 (MQ=255)
                                                                                                        ttCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACTGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                                         <  2:132189/149‑1 (MQ=255)
                                                                                                          cGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACCCCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGc                                        >  5:222300/1‑148 (MQ=255)
                                                                                                           gATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                      >  1:202599/1‑149 (MQ=255)
                                                                                                           gATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGAGAGGAAGCTGGATGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAAAGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                      >  1:207654/1‑149 (MQ=255)
                                                                                                            aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAt                                     <  2:200273/149‑1 (MQ=255)
                                                                                                             tAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATg                                    <  4:202965/149‑1 (MQ=255)
                                                                                                               aaTACGCGGACTTCGCCCTCTTCTTCCACGAGGACCGCGCGCAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCAACGTCGTAACGCACCACACCAACGGTGACGTGTACTGCGTTCTCGACTTCCGCCTGCGAGGAGGCGATgg                                   >  8:112115/1‑148 (MQ=255)
                                                                                                                  aCGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATa                               >  7:286832/1‑149 (MQ=255)
                                                                                                                    gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCCCGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATa                               >  3:206333/1‑147 (MQ=255)
                                                                                                                    gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGACACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                             <  3:290801/149‑1 (MQ=255)
                                                                                                                     cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAg                            <  7:331323/149‑1 (MQ=255)
                                                                                                                       gACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAc                          >  7:275145/1‑149 (MQ=255)
                                                                                                                       gACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAc                          <  4:232889/149‑1 (MQ=255)
                                                                                                                          ttCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgcg                       >  6:113841/1‑149 (MQ=255)
                                                                                                                          ttCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgcg                       >  5:75154/1‑149 (MQ=255)
                                                                                                                            cGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgcg                       <  8:275145/147‑1 (MQ=255)
                                                                                                                             gCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggc                     <  5:179730/148‑1 (MQ=255)
                                                                                                                              cccTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGt                   <  6:222300/149‑1 (MQ=255)
                                                                                                                                ctctTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGt                   >  7:70471/1‑147 (MQ=255)
                                                                                                                                 tcttctTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAAGAGGCGATGGAATCAAGACGCGGCGTCAg                >  8:238764/1‑149 (MQ=255)
                                                                                                                                  cttcttCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGCCGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGc               >  1:301897/1‑149 (MQ=255)
                                                                                                                                   ttcttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAg                <  6:338183/147‑1 (MQ=255)
                                                                                                                                      ttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCAACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCAGCCTGCGAGGAGGCGATGGAATACAGACGCGTCGTCAGCGGgc           >  5:107300/1‑149 (MQ=255)
                                                                                                                                      ttcCACGCGGACAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTACTCGACTTACGCCAGCGAGCAGGCGATGGAAGACAGACGCGGCGTCAGCGGgc           >  4:214543/1‑149 (MQ=255)
                                                                                                                                       tcCCCGCGGTCAGCGCGGAAGCTGGAGGCCCGCCCGGCGCGGGCGCGGCCTGCCCCGTCGTACCGGGCCACCCCAACGGTGACGTGTACTTCGGGCGCGCCTTCCGCCTGCGAGGCGGCGAGGGAAGACAGACGCGGCGTCAGCGGccg          >  5:79521/1‑149 (MQ=255)
                                                                                                                                       tcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGCCGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGgcg          >  5:79519/1‑149 (MQ=255)
                                                                                                                                       tcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGgcg          >  4:176406/1‑149 (MQ=255)
                                                                                                                                       tcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGACTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGgcg          >  5:132782/1‑149 (MQ=255)
                                                                                                                                         cACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgca        <  7:248652/149‑1 (MQ=255)
                                                                                                                                         cACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGCTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGgcg          >  7:276407/1‑147 (MQ=255)
                                                                                                                                               gTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGAACAACACCAACGGTGACGTTTACTGCGTTCTCGACTTCCGCCTGCTCAGATGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGat  <  7:238764/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGAT  >  NC_000913/2890475‑2890766

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