Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,050,408 A→G S264G (AGT→GGT)  hemN → O2‑independent coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,050,4080AG99.4% 1489.1 / ‑2.9 470S264G (AGT→GGT) hemNO2‑independent coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (252/214);  total (254/216)
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.57e-01

Only 100 of 486 total aligned reads displayed.

CTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGT  >  NC_000913/4050318‑4050500
                                                                                          |                                                                                            
  gAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTcc                                                                                          >  2:3568041/1‑93 (MQ=255)
  gAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTcc                                                                                          <  2:4250071/93‑1 (MQ=255)
      cccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCgcag                                                                                      >  2:1255046/1‑93 (MQ=255)
       ccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCgcagc                                                                                     >  1:8039232/1‑93 (MQ=255)
       ccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCgcagc                                                                                     <  1:9540851/93‑1 (MQ=255)
       ccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCgcagc                                                                                     >  2:8748353/1‑93 (MQ=255)
           tCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCaaaa                                                                                 >  1:7782666/1‑93 (MQ=255)
           tCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCaaaa                                                                                 <  1:2213192/93‑1 (MQ=255)
            cGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCaaaaa                                                                                <  1:4295058/93‑1 (MQ=255)
             gTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAAc                                                                               >  2:9025710/1‑93 (MQ=255)
              tCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACt                                                                              >  2:6801491/1‑93 (MQ=255)
              tCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACt                                                                              <  2:3199718/93‑1 (MQ=255)
              tCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACt                                                                              >  2:529999/1‑93 (MQ=255)
              tCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGCCCGCAGCAAAAACt                                                                              >  2:7118152/1‑93 (MQ=255)
                tGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCg                                                                            >  2:9232211/1‑93 (MQ=255)
                 gattgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGa                                                                           <  2:1996128/90‑1 (MQ=255)
                 gAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGa                                                                           >  1:6778129/1‑93 (MQ=255)
                 gAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGa                                                                           >  2:5995361/1‑93 (MQ=255)
                 gAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGa                                                                           <  1:1751250/93‑1 (MQ=255)
                  aGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGat                                                                          >  1:9786959/1‑93 (MQ=255)
                  aGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGat                                                                          <  2:9447653/93‑1 (MQ=255)
                  aGTATCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGat                                                                          >  1:1147382/1‑93 (MQ=255)
                   gtgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGata                                                                         <  2:266741/93‑1 (MQ=255)
                   gtgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGata                                                                         <  1:2787244/93‑1 (MQ=255)
                   gtgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGata                                                                         <  2:8205085/93‑1 (MQ=255)
                    tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGatat                                                                        >  2:1864487/1‑93 (MQ=255)
                    tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGatat                                                                        <  1:8502367/93‑1 (MQ=255)
                    tgtCTTTAACTACGCGCATCTGCAGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTTACCTGCCGGGTCCGCAGCAAAAACCCGatat                                                                        >  2:3497580/1‑93 (MQ=255)
                     gtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATc                                                                       <  2:9948300/93‑1 (MQ=255)
                      tCTTTAACTCCGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATcc                                                                      >  2:6731623/1‑93 (MQ=255)
                      tCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATcc                                                                      <  2:4783115/93‑1 (MQ=255)
                      gcgttAAATACGCGCATCTGCCGACCATGTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATcc                                                                      <  1:5741328/90‑1 (MQ=255)
                       cTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCt                                                                     <  1:6660352/93‑1 (MQ=255)
                       cTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCt                                                                     >  2:3894360/1‑93 (MQ=255)
                          tAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCa                                                                  >  1:3583725/1‑93 (MQ=255)
                             cTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGa                                                               >  2:3989006/1‑93 (MQ=255)
                              tACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGaa                                                              >  2:6205224/1‑93 (MQ=255)
                                cgcgCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAc                                                            >  1:443301/1‑93 (MQ=255)
                                 gcgcATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAcc                                                           >  1:7076525/1‑93 (MQ=255)
                                    cATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATc                                                        <  1:10024447/93‑1 (MQ=255)
                                    cATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATc                                                        >  2:5798722/1‑93 (MQ=255)
                                    cATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATCCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATc                                                        <  2:1212145/93‑1 (MQ=255)
                                      tCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGc                                                      >  1:4879252/1‑93 (MQ=255)
                                      tCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGc                                                      <  1:9249449/93‑1 (MQ=255)
                                      tCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGc                                                      >  1:8631587/1‑93 (MQ=255)
                                        tGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCt                                                    >  1:4744246/1‑93 (MQ=255)
                                          ccGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTc                                                  <  2:1564286/93‑1 (MQ=255)
                                            gACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCt                                                <  2:1418069/93‑1 (MQ=255)
                                             aCCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTg                                               >  1:3091633/1‑93 (MQ=255)
                                              ccATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGa                                              <  2:6502730/93‑1 (MQ=255)
                                              ccATTTTTGCTGCTCAGCGCAAAATCAAAGATCCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGa                                              <  1:5167122/93‑1 (MQ=255)
                                                aTTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACg                                            >  1:7233426/1‑93 (MQ=255)
                                                aTTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACg                                            >  1:6014956/1‑93 (MQ=255)
                                                 tttttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGc                                           >  1:9895588/1‑93 (MQ=255)
                                                 tttttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGc                                           >  2:2577462/1‑93 (MQ=255)
                                                    ttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAAt                                        <  2:2208514/93‑1 (MQ=255)
                                                    ttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAAt                                        <  1:8733777/93‑1 (MQ=255)
                                                    ttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAAt                                        <  1:2140595/93‑1 (MQ=255)
                                                     tgctgcTCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAACATCGCCTTCCTGACGCAATc                                       >  2:5679381/1‑93 (MQ=255)
                                                      gctgctCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCg                                      <  1:4973467/93‑1 (MQ=255)
                                                      gctgctCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCg                                      <  2:4408233/93‑1 (MQ=255)
                                                      gctgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCg                                      <  1:1887237/93‑1 (MQ=255)
                                                       ctgctCAGCGCAAAATCAAAGATGCTGACCTGCCAGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCgg                                     <  2:6923016/93‑1 (MQ=255)
                                                         gctCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGc                                   >  2:3371436/1‑93 (MQ=255)
                                                         gctCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGc                                   >  1:62422/1‑93 (MQ=255)
                                                         gctCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGc                                   >  1:2374830/1‑93 (MQ=255)
                                                           tCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTa                                 >  2:3685960/1‑93 (MQ=255)
                                                           tCAGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTa                                 >  2:496385/1‑93 (MQ=255)
                                                             aGCGCAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATc                               <  2:4579533/93‑1 (MQ=255)
                                                              gcgcAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCa                              <  1:3290890/93‑1 (MQ=255)
                                                              gcgcAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCa                              >  1:1522081/1‑93 (MQ=255)
                                                                gcAAAATCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGt                            >  2:7972268/1‑93 (MQ=255)
                                                                  aaaaTCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGttt                          <  1:644568/93‑1 (MQ=255)
                                                                  aaaaTCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGttt                          >  1:9017346/1‑93 (MQ=255)
                                                                   aaaTCAGAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTACCTGACGCAATCGGGCTATCAGTTTa                         >  2:7318394/1‑93 (MQ=255)
                                                                   aaaTCAAAGATGCTGACCTGCCGGGTCCGCAGCTAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                         <  1:9682343/93‑1 (MQ=255)
                                                                   aaaTCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                         >  1:2483157/1‑93 (MQ=255)
                                                                   aaaTCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACCCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                         <  1:7008057/93‑1 (MQ=255)
                                                                     aTCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATc                       <  2:6659027/93‑1 (MQ=255)
                                                                      tCAAAGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCg                      >  1:1635409/1‑93 (MQ=255)
                                                                         aaGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTa                   >  1:3735092/1‑93 (MQ=255)
                                                                          aGATGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTAt                  >  2:8327624/1‑93 (MQ=255)
                                                                           gATGCTAACCTGCCGGGTCCGCAGCAAAAAATCGATATCCTGCAGGAAACCATCGCCTTCCTGAAGCAATCGGGCTATCAGTTTATCGGTATg                 >  2:5149056/1‑93 (MQ=255)
                                                                             tGCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGa               <  2:1017111/93‑1 (MQ=255)
                                                                              gCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGAt              <  1:7527512/93‑1 (MQ=255)
                                                                              gCTGACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGAt              <  1:5861428/93‑1 (MQ=255)
                                                                                 gACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCAc           >  2:2747475/1‑93 (MQ=255)
                                                                                 gACCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCAc           >  2:968631/1‑93 (MQ=255)
                                                                                  aCCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGACACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACt          <  2:6324399/93‑1 (MQ=255)
                                                                                  aCCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACt          >  1:5476181/1‑93 (MQ=255)
                                                                                  aCCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACt          <  2:8471586/93‑1 (MQ=255)
                                                                                  aCCTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACt          >  1:8354448/1‑93 (MQ=255)
                                                                                    cTGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACttt        <  1:9121047/93‑1 (MQ=255)
                                                                                     tGCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTg       >  2:3002557/1‑93 (MQ=255)
                                                                                      gCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc      <  1:3874616/93‑1 (MQ=255)
                                                                                      gCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc      >  1:4138608/1‑93 (MQ=255)
                                                                                      gCCGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc      <  1:7352637/93‑1 (MQ=255)
                                                                                       ccGGGTCCGCAGCAAAACCTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGcc     <  2:7622849/93‑1 (MQ=255)
                                                                                       ccGGGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGcc     >  2:7781813/1‑93 (MQ=255)
                                                                                          aGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGt  >  2:8522168/1‑93 (MQ=255)
                                                                                          |                                                                                            
CTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGT  >  NC_000913/4050318‑4050500

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