Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,888,575 G→T 8.5% R509S (CGC→AGC)  cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,888,5750GT 8.5% 236.1 / 5.8 70R509S (CGC→AGC) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base T (3/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTT  >  NC_000913/2888444‑2888710
                                                                                                                                   |                                                                                                                                       
cttacaGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCCCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGCGGGATACGTGTCCCAATGCGGTTgccgc                                                                                                                                 <  2:834496/135‑1 (MQ=255)
          ccGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAg                                                                                                                       >  2:1173730/1‑140 (MQ=255)
          ccGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTGCATCAAGCGACGCCAGGATTTCCGGCTCGGTGAGGTTGTCTTTATACAACCGTGGGATACGTGTCCCAATGCGGTTGCCGCCACGATGCAg                                                                                                                       >  2:964393/1‑140 (MQ=255)
          ccGAAGCCTTCACCCGCATCGCGCTCGTTCGCCCAGCTCCATATCATTTCTTCTAGAGACGCCAGGATTTCCGGCTCTGTGTTGTTTTGTTTATACTTCCGTGGGATACTTGTCCCAATGCGGTTGCCGCCAAGATGCAg                                                                                                                       >  2:803510/1‑140 (MQ=255)
                   tCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCTGTTGCCGCCAAGATGCAGGTTGTAGCg                                                                                                              >  1:894014/1‑140 (MQ=255)
                   tCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCTGTTGCCGCCAAGATGCAGGTTGTAGCg                                                                                                              >  1:894236/1‑140 (MQ=255)
                   tCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCTGTTGCCGCCAAGATGCAGGTTGTAGCg                                                                                                              >  1:292178/1‑140 (MQ=255)
                      cccGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGAcc                                                                                                           <  1:179/140‑1 (MQ=255)
                      cccGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGAcc                                                                                                           <  1:261/140‑1 (MQ=255)
                      cccGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGAcc                                                                                                           <  1:107/140‑1 (MQ=255)
                      cccGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGAcc                                                                                                           <  1:1173730/140‑1 (MQ=255)
                      cccGCTTCGCGCTATTTCGCCCAGCGCCCTATCAGTTAATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCACAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGAcc                                                                                                           <  1:14598/140‑1 (MQ=255)
                       ccGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGAccc                                                                                                          <  2:743650/140‑1 (MQ=255)
                       ccGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGAccc                                                                                                          <  2:743669/140‑1 (MQ=255)
                        cGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCg                                                                                                         <  1:334484/140‑1 (MQ=255)
                         gCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCgg                                                                                                        <  2:551684/140‑1 (MQ=255)
                         gCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCgg                                                                                                        <  2:292823/140‑1 (MQ=255)
                          cTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGgc                                                                                                       >  2:254999/1‑140 (MQ=255)
                          cTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGgc                                                                                                       >  2:256039/1‑140 (MQ=255)
                          cTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGgc                                                                                                       >  2:849013/1‑140 (MQ=255)
                          cTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGgc                                                                                                       >  2:544569/1‑140 (MQ=255)
                               cgcTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCttt                                                                                                  >  1:158853/1‑140 (MQ=255)
                               cgcTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCttt                                                                                                  >  2:59225/1‑140 (MQ=255)
                               cgcTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTGTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCttt                                                                                                  >  1:768575/1‑140 (MQ=255)
                                gcTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTa                                                                                                 <  1:308535/140‑1 (MQ=255)
                                gcTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTATTTATACATCCGTGGGATACGTGTCCCAATGCGGTTTCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTa                                                                                                 <  1:1176592/140‑1 (MQ=255)
                                                        gTTCATCAATCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGTGATACTTGTCCGAATTCGGTAGCCGCCAAGAAGCAGGTTGTAGAGACCCGGCGCGTTACACACCAGGCCCACTTccgcccgc                                                                         >  2:683605/1‑140 (MQ=255)
                                                        gTTCATCAAGCGACGCCAGGATATCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAACATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTACCGCCAGc                                                                         >  2:836854/1‑140 (MQ=255)
                                                           catcaAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATc                                                                      >  1:630011/1‑140 (MQ=255)
                                                             tcaAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATcgc                                                                    >  2:163001/1‑140 (MQ=255)
                                                                               tCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCgg                                                  >  1:765526/1‑140 (MQ=255)
                                                                               tCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCgg                                                  >  2:33383/1‑140 (MQ=255)
                                                                               tCCGGCTCGGTGATGTTTTCTTTATACAACCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCgg                                                  >  2:33452/1‑140 (MQ=255)
                                                                                  ggCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCa                                               <  1:330018/140‑1 (MQ=255)
                                                                                  ggCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCa                                               <  2:1008770/140‑1 (MQ=255)
                                                                                         tGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGtt                                        <  2:1187553/140‑1 (MQ=255)
                                                                                            tGTTTTCTTTCTACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTaca                                     <  1:726152/140‑1 (MQ=255)
                                                                                            tGTTTTCTTTATTCATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTaca                                     <  1:607208/140‑1 (MQ=255)
                                                                                            tGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTaca                                     <  1:1157074/140‑1 (MQ=255)
                                                                                            tGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTaca                                     <  1:1204519/140‑1 (MQ=255)
                                                                                            tGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTaca                                     <  1:203070/140‑1 (MQ=255)
                                                                                            tGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTaca                                     <  1:799002/140‑1 (MQ=255)
                                                                                              ttttCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACg                                   <  1:316790/140‑1 (MQ=255)
                                                                                                 tCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCAt                                <  1:364830/140‑1 (MQ=255)
                                                                                                  cTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATc                               <  1:809606/140‑1 (MQ=255)
                                                                                                  cTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATc                               <  1:809557/140‑1 (MQ=255)
                                                                                                             cGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTc                    <  1:465055/140‑1 (MQ=255)
                                                                                                               tGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCAt                  >  1:135129/1‑140 (MQ=255)
                                                                                                               tGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCAt                  >  1:135303/1‑140 (MQ=255)
                                                                                                                 ggATTCGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGCTGTAGCGACCCGGCGCTTTACCCACCAGGCCCTCTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCg                <  1:818291/140‑1 (MQ=255)
                                                                                                                 ggATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCg                <  2:534590/140‑1 (MQ=255)
                                                                                                                 ggATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCg                <  1:907425/140‑1 (MQ=255)
                                                                                                                 ggATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCg                <  2:551902/140‑1 (MQ=255)
                                                                                                                    tACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCAGTTACACGCATCACGATATGCTCATCGCTg             >  2:949687/1‑140 (MQ=255)
                                                                                                                      cGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACACCCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGac           >  2:918299/1‑140 (MQ=255)
                                                                                                                      cGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGac           >  2:239299/1‑140 (MQ=255)
                                                                                                                      cGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGac           >  1:533011/1‑140 (MQ=255)
                                                                                                                      cGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGac           >  1:62362/1‑140 (MQ=255)
                                                                                                                      cGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGac           >  2:239198/1‑140 (MQ=255)
                                                                                                                        tgtCCCAATGCTGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGacac         <  2:292178/140‑1 (MQ=255)
                                                                                                                        tgtCCCAATGCTGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCACTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGacac         <  2:894014/140‑1 (MQ=255)
                                                                                                                        tgtCCCAATGCTGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCACTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGacac         <  2:894236/140‑1 (MQ=255)
                                                                                                                            ccAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATg     >  1:73296/1‑140 (MQ=255)
                                                                                                                            ccAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATg     >  1:641027/1‑140 (MQ=255)
                                                                                                                            ccAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATg     >  1:73294/1‑140 (MQ=255)
                                                                                                                            ccAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATg     >  1:737226/1‑140 (MQ=255)
                                                                                                                            ccAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATg     >  1:73758/1‑140 (MQ=255)
                                                                                                                            ccAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATg     >  1:963125/1‑140 (MQ=255)
                                                                                                                               aTGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGttt  <  2:1158177/140‑1 (MQ=255)
                                                                                                                               aTGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGttt  >  2:414068/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                       
CGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTTT  >  NC_000913/2888444‑2888710

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