Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,430,453 +C coding (745/747 nt) smf ← putative Rossmann‑fold nucleotide‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,430,4531.C99.5% 676.3 / ‑5.8 198P249A (CCT→GCT) smfputative Rossmann‑fold nucleotide‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (113/84);  minor base A (1/0);  total (114/84)
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.97e-01

Only 100 of 219 total aligned reads displayed.

GAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTT  >  NC_000913/3430361‑3430539
                                                                                             |                                                                                      
gAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAgg                                                                                         >  1:6066022/1‑93 (MQ=255)
 aaTTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAggc                                                                                        <  1:5831460/93‑2 (MQ=255)
 aaTTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAggc                                                                                        <  2:4479093/93‑2 (MQ=255)
    ttcttcCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGcaaa                                                                                    >  2:4455860/1‑93 (MQ=255)
     tcttcCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGcaaaa                                                                                   >  1:414793/1‑93 (MQ=255)
        tcCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACt                                                                                 <  2:2371723/93‑1 (MQ=255)
        tcCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACt                                                                                 <  2:5443200/93‑1 (MQ=255)
          cGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtc                                                                               >  1:4280607/1‑93 (MQ=255)
           ggTTCCGTCCCAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtct                                                                              >  2:5627205/1‑93 (MQ=255)
           ggTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtct                                                                              >  2:5766428/1‑93 (MQ=255)
           ggTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtct                                                                              <  2:4120222/93‑1 (MQ=255)
           ggTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtct                                                                              <  1:5016328/93‑1 (MQ=255)
            gTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                                             >  1:4729104/1‑93 (MQ=255)
            gTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                                             >  2:3037297/1‑93 (MQ=255)
            gTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTGCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                                             <  2:680169/93‑1 (MQ=255)
             ttCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCg                                                                            <  1:5516012/93‑1 (MQ=255)
             ttCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCg                                                                            >  2:2892114/1‑93 (MQ=255)
               ccGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGcc                                                                          <  1:1341626/93‑1 (MQ=255)
                 gccacaAGAATCGCACCTTTTTTTATTAACCAGGGAGGCCCTTCGCTCCCCGGATTCCCTATTGGCCCTGGCAAGGCAAAAACTTCTCGcccc                                                                        <  1:1202705/91‑1 (MQ=255)
                 gTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGcccc                                                                        >  1:320839/1‑93 (MQ=255)
                  tCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCt                                                                       <  2:6342882/93‑1 (MQ=255)
                  tCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCt                                                                       <  2:884936/93‑1 (MQ=255)
                  tCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGAAAAAACTTCTCGCCCCt                                                                       >  2:1621786/1‑93 (MQ=255)
                   cacaAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTg                                                                      >  1:2090438/1‑93 (MQ=255)
                    acaAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGc                                                                     >  2:358525/1‑93 (MQ=255)
                    acaAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGc                                                                     >  2:6384495/1‑93 (MQ=255)
                    acaAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGc                                                                     >  1:5636801/1‑93 (MQ=255)
                      aaGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTc                                                                   <  2:2112395/93‑1 (MQ=255)
                       aGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCa                                                                  >  1:2655270/1‑93 (MQ=255)
                        gAATCGCCCCTTGTTCTATTAAGCAGTGAGGCCCTTCGCTCCCCGGATTCCCTAGGGGTCCGGGCAAGGCAAAAACTCCTCGCCCCTGCTaaa                                                                 <  1:1065301/93‑1 (MQ=255)
                         aaTCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAg                                                                >  1:277219/1‑93 (MQ=255)
                          aTCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAgc                                                               >  1:4963308/1‑93 (MQ=255)
                            cGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAgcgc                                                             >  2:554246/1‑93 (MQ=255)
                             gCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGca                                                            >  2:3868578/1‑93 (MQ=255)
                             gCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGca                                                            >  1:2172378/1‑93 (MQ=255)
                              cACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGcac                                                           >  2:1291995/1‑93 (MQ=255)
                               aCCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGcaca                                                          >  2:1538961/1‑93 (MQ=255)
                                 cTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAAc                                                        <  2:5937915/93‑1 (MQ=255)
                                  ttGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACg                                                       >  2:3020838/1‑93 (MQ=255)
                                  ttGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACg                                                       <  1:4035155/93‑1 (MQ=255)
                                  ttGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACg                                                       <  2:4953696/93‑1 (MQ=255)
                                   tGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGt                                                      <  2:4238446/93‑1 (MQ=255)
                                   tGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGt                                                      <  1:744151/93‑1 (MQ=255)
                                    gTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTg                                                     >  1:5098118/1‑93 (MQ=255)
                                    gTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTg                                                     >  1:1471989/1‑93 (MQ=255)
                                     ttttATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGc                                                    <  1:5772657/93‑1 (MQ=255)
                                        tattaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGt                                                 <  2:11353/93‑1 (MQ=255)
                                         attaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTc                                                >  1:2930015/1‑93 (MQ=255)
                                         attaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTc                                                >  2:1339839/1‑93 (MQ=255)
                                          ttaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCa                                               >  2:1404581/1‑93 (MQ=255)
                                           taACCAGTGAGCCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCAc                                              <  1:4534914/93‑1 (MQ=255)
                                            aaCCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTTCTGTCAcc                                             <  1:873671/93‑1 (MQ=255)
                                             aCCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCa                                            <  2:3219510/93‑1 (MQ=255)
                                               cAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGc                                          <  1:3764549/93‑1 (MQ=255)
                                               cAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGc                                          <  1:4361275/93‑1 (MQ=255)
                                               cAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGc                                          <  1:1058877/93‑1 (MQ=255)
                                                aGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCg                                         >  2:812589/1‑93 (MQ=255)
                                                aGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCg                                         >  2:4327743/1‑93 (MQ=255)
                                                aGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCg                                         >  1:4803108/1‑93 (MQ=255)
                                                aGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCg                                         >  2:4805674/1‑93 (MQ=255)
                                                 gTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCGAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGa                                        >  2:6177174/1‑93 (MQ=255)
                                                   gAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAACCTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAAc                                      <  1:2285009/93‑1 (MQ=255)
                                                   gAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACCACGTGCTGTCACCAGCGAAc                                      >  1:3367122/1‑93 (MQ=255)
                                                      gCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGGGCTGTCACCAGCGAACCAc                                   >  1:2353970/1‑93 (MQ=255)
                                                       cccTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACt                                  >  2:1965907/1‑93 (MQ=255)
                                                          ttCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACg                               >  1:5204388/1‑93 (MQ=255)
                                                          ttCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACg                               >  2:6268734/1‑93 (MQ=255)
                                                           tCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGc                              <  2:3701165/93‑1 (MQ=255)
                                                           tCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGc                              >  1:825306/1‑93 (MQ=255)
                                                              cTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCaaa                           >  1:3077342/1‑93 (MQ=255)
                                                              cTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCaaa                           >  2:5240525/1‑93 (MQ=255)
                                                              cTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCaaa                           <  1:1425072/93‑1 (MQ=255)
                                                               ttCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAg                          >  2:2431606/1‑93 (MQ=255)
                                                                tCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGc                         <  1:4880049/93‑1 (MQ=255)
                                                                tCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGc                         >  2:5214452/1‑93 (MQ=255)
                                                                 cccGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGcc                        <  2:5125013/93‑1 (MQ=255)
                                                                   cGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGc                      >  2:750204/1‑93 (MQ=255)
                                                                    ggATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCt                     >  1:1911725/1‑93 (MQ=255)
                                                                     gATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCtt                    >  2:3320576/1‑93 (MQ=255)
                                                                      aTTCCCTATTGGACCTGGCAAGGCAAAACCTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTc                   <  2:4420535/93‑1 (MQ=255)
                                                                      aTTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTc                   >  2:935229/1‑93 (MQ=255)
                                                                       ttCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTcc                  <  2:2688920/93‑1 (MQ=255)
                                                                       ttCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTcc                  >  2:6118372/1‑93 (MQ=255)
                                                                       ttCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCACAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTcc                  >  2:3449063/1‑93 (MQ=255)
                                                                         cccTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTccac                <  1:4462282/93‑1 (MQ=255)
                                                                          ccTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTTTCACCAGCGAACCACTACGCAAAGCCGCTTccacc               <  1:1777204/93‑1 (MQ=255)
                                                                           cTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTccacca              >  2:2181176/1‑93 (MQ=255)
                                                                             aTTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTccaccacc            >  1:5574500/1‑93 (MQ=255)
                                                                             aTTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTccaccacc            >  1:5081245/1‑93 (MQ=255)
                                                                              ttGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTccaccacca           <  2:1008013/93‑1 (MQ=255)
                                                                              ttGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTccaccacca           <  2:528118/93‑1 (MQ=255)
                                                                                ggACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGt         >  2:1778906/1‑93 (MQ=255)
                                                                                ggACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGACGCTTCCACCACCAGt         >  1:5916427/1‑93 (MQ=255)
                                                                                 gACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCGACCACCAgca        >  1:356381/1‑91 (MQ=255)
                                                                                   ccTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTaca      >  1:5280661/1‑93 (MQ=255)
                                                                                   ccTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTaca      >  1:374747/1‑93 (MQ=255)
                                                                                     tGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACAcc    <  1:6412122/93‑1 (MQ=255)
                                                                                      ggcaaggcaaAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCt   >  2:2542471/1‑93 (MQ=255)
                                                                                       gCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCtt  >  1:188114/1‑93 (MQ=255)
                                                                                       gCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCtt  <  1:5700401/93‑1 (MQ=255)
                                                                                             |                                                                                      
GAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTT  >  NC_000913/3430361‑3430539

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