Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,891,288 G→A 5.9% A204V (GCG→GTG)  cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,891,2880GA 5.9% 220.8 / 2.3 68A204V (GCG→GTG) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base A (2/2);  total (32/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 = 5.51e-01

GAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTA  >  NC_000913/2891144‑2891434
                                                                                                                                                |                                                                                                                                                  
gagTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCAt                                                                                                                                                 <  1:185533/148‑1 (MQ=255)
   ttACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGCCGGAGTGCGAAGGTGCCGCGACAGGCGCACGCGATTTCAGCgcagcac                                                                                                                                              >  4:86336/1‑147 (MQ=255)
    tACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAAc                                                                                                                                             >  6:311491/1‑148 (MQ=255)
           ggTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCgg                                                                                                                                      >  3:156196/1‑148 (MQ=255)
            gTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCggg                                                                                                                                     >  3:89828/1‑148 (MQ=255)
              aaTTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGAGTCTTTGCTGTACGGGCTGGGGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGCCGCACGCGATTTAAGCGGATCAACCACGCGGgcg                                                                                                                                   >  4:40825/1‑148 (MQ=255)
                  ttCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGgcgcgc                                                                                                                                <  2:163262/147‑1 (MQ=255)
                       gTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCg                                                                                                                          <  4:89828/148‑1 (MQ=255)
                        ttAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGc                                                                                                                         <  8:115152/148‑1 (MQ=255)
                             agagagGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGcc                                                                                                                               <  1:12737/137‑1 (MQ=255)
                             agagagGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGcc                                                                                                                               >  2:12737/1‑137 (MQ=255)
                               agagGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGgcagc                                                                                                                  <  1:33265/148‑1 (MQ=255)
                               agagGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCACATCAACCACGCGGGCGCGCCACTCGc                                                                                                                         <  8:7570/141‑1 (MQ=255)
                               agagGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCACATCAACCACGCGGGCGCGCCACTCGc                                                                                                                         >  7:7570/1‑141 (MQ=255)
                                    cTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCt                                                                                                             <  3:277707/148‑1 (MQ=255)
                                       gCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTgg                                                                                                           >  4:291327/1‑147 (MQ=255)
                                       gCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGt                                                                                                          <  1:67410/148‑1 (MQ=255)
                                       gCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGt                                                                                                          <  7:232172/148‑1 (MQ=255)
                                       gCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGt                                                                                                          <  3:23929/148‑1 (MQ=255)
                                         caccaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCCCTCGCTGGCAGCAGCCTGGTAt                                                                                                        >  2:354456/1‑148 (MQ=255)
                                         caccaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTAt                                                                                                        >  1:262383/1‑148 (MQ=255)
                                          accaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTg                                                                                                                       >  4:304372/1‑132 (MQ=255)
                                          accaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTg                                                                                                                       <  3:304372/132‑1 (MQ=255)
                                          accaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTAtt                                                                                                       <  1:340830/148‑1 (MQ=255)
                                          accaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTAtt                                                                                                       <  7:366315/148‑1 (MQ=255)
                                           ccaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTc                                                                                                      <  3:315728/148‑1 (MQ=255)
                                           ccaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCAACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTAtt                                                                                                       >  5:332112/1‑147 (MQ=255)
                                            caGCGGCGCGTCTTTGCTGTACTGGCTTGTGTGGATTTCATTTACCGCGCCAGTAGCGCCGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAATCGCATCAACCAGGCGGGCGCGCCTCTCGCTGGCAGCAGCATGGTATTCa                                                                                                     <  3:40825/148‑1 (MQ=255)
                                            caGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTc                                                                                                      >  3:292529/1‑147 (MQ=255)
                                                         ttGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGa                                                                                        >  2:187394/1‑148 (MQ=255)
                                                          tGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATc                                                                                          <  1:26170/145‑1 (MQ=255)
                                                               tACGGGCTGGTGTTGATTTCATTTACCGCGCCATTAGCTACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCAAGCGGGCGCGCCACTAGCTGGCAGCAGCCTTGTATTCAACATCGGCAACGACACgg                                                                                   <  1:312286/147‑1 (MQ=255)
                                                               tACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGt                                                                                  <  5:311491/148‑1 (MQ=255)
                                                                  ggggTGGTGGGGTTTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCTCGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTTTTCAACATCGGCATCGACACGGTCGa                                                                               <  4:292529/148‑1 (MQ=255)
                                                                       ggTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAgg                                                                          >  3:118149/1‑148 (MQ=255)
                                                                         tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTc                                                                                                                           <  3:131104/97‑1 (MQ=255)
                                                                         tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTc                                                                                                                           >  4:131104/1‑97 (MQ=255)
                                                                              aTTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTcac                                                                   >  8:283188/1‑148 (MQ=255)
                                                                               tttCATTTACCGCGACAGTAGCGACGGATTGCGAAGCTGCCGCGACAGGCGCACGCGATATAAGCGAATCAAACAAGCGGGCCCGCCACTCGCTGGCAGCAGCCTGCTATTCAACAACTCCATCGACACGGACGAGCAGGCGTTcac                                                                   >  2:239167/1‑147 (MQ=255)
                                                                                ttCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTcacc                                                                  <  2:262383/147‑1 (MQ=255)
                                                                                   aTTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCa                                                                                                     <  6:279786/109‑1 (MQ=255)
                                                                                   aTTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCa                                                                                                     >  5:279786/1‑109 (MQ=255)
                                                                                   aTTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCAccc                                                              <  1:187394/148‑1 (MQ=255)
                                                                                    tttACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCa                                                             <  6:165151/148‑1 (MQ=255)
                                                                                           cgcCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGc                                                      >  3:234804/1‑148 (MQ=255)
                                                                                              cAGTAGCGACGGATTGCGAAGGTGCCGCAACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAg                                                   <  6:332112/148‑1 (MQ=255)
                                                                                                  aGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCa                                                                                                     >  5:123955/1‑94 (MQ=255)
                                                                                                  aGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCa                                                                                                     <  6:123955/94‑1 (MQ=255)
                                                                                                      aCGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGt                                           <  8:387568/148‑1 (MQ=255)
                                                                                                          aTTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGa                                        >  5:58921/1‑147 (MQ=255)
                                                                                                          aTTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGa                                        >  5:345440/1‑147 (MQ=255)
                                                                                                               gAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCtt                                  >  8:268949/1‑148 (MQ=255)
                                                                                                                aaGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAActctgt                                 >  2:239143/1‑146 (MQ=255)
                                                                                                                 aGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTc                                <  6:345440/148‑1 (MQ=255)
                                                                                                                 aGGTGCCGCGACAGGCGCACGCGATTTAAGCGCAGCACCCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTc                                <  3:86336/148‑1 (MQ=255)
                                                                                                                    tGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCCCTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCg                             >  6:299600/1‑148 (MQ=255)
                                                                                                                    tGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCa                                                             >  5:257817/1‑116 (MQ=255)
                                                                                                                    tGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCa                                                             <  6:257817/116‑1 (MQ=255)
                                                                                                                    tGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCg                             >  7:308118/1‑148 (MQ=255)
                                                                                                                    tGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCAGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCg                             >  4:65487/1‑148 (MQ=255)
                                                                                                                            cAGGCGCACGCGATTTAAGCACATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTc                                                                     >  7:318902/1‑100 (MQ=255)
                                                                                                                            cAGGCGCACGCGATTTAAGCACATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTc                                                                     <  8:318902/100‑1 (MQ=255)
                                                                                                                               gcgcACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGa                   >  2:101008/1‑147 (MQ=255)
                                                                                                                                    cgcgATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAAtt             >  6:131861/1‑148 (MQ=255)
                                                                                                                                        aTTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATa         <  6:58921/148‑1 (MQ=255)
                                                                                                                                           tAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAag       <  4:118149/147‑1 (MQ=255)
                                                                                                                                            aaGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGc                                              <  1:203202/107‑1 (MQ=255)
                                                                                                                                            aaGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGc                                              >  2:203202/1‑107 (MQ=255)
                                                                                                                                            aaGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAaga      <  8:308118/147‑1 (MQ=255)
                                                                                                                                               cgcATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTa  >  8:301686/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTA  >  NC_000913/2891144‑2891434

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