Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,133,826 G→C 8.7% R609R (CGC→CGG)  wzc ← protein‑tyrosine kinase Wzc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,133,8260GC 8.7% 215.1 / 10.3 69R609R (CGC→CGG) wzcprotein‑tyrosine kinase Wzc
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base C (2/4);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

ACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCA  >  NC_000913/2133682‑2133974
                                                                                                                                                |                                                                                                                                                    
acgacCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCgg                                                                                                                                                 <  2:326239/150‑1 (MQ=255)
 cgacCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCggg                                                                                                                                                >  2:1055730/1‑150 (MQ=255)
    ccAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATc                                                                                                                                             >  1:398432/1‑150 (MQ=255)
    ccAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATc                                                                                                                                             >  1:1237105/1‑150 (MQ=255)
       acaaTTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCTGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAgg                                                                                                                                          <  2:473351/150‑1 (MQ=255)
          aTTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCa                                                                                                                                       >  1:837100/1‑150 (MQ=255)
            tGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGATCCCGCGCGGTCATGtgccg                                                                                                                                     <  1:1292060/150‑5 (MQ=255)
            tGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGATCCCGCGCGGTCAGGtgccg                                                                                                                                     <  1:578968/150‑5 (MQ=255)
            tGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGATCCCGCGCGGTCAGGtgccg                                                                                                                                     <  1:151487/150‑5 (MQ=255)
                gCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTg                                                                                                                                 >  2:1253731/1‑150 (MQ=255)
                gCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTg                                                                                                                                 <  2:337593/150‑1 (MQ=255)
                gCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACAGCGCGGGATCAGGTCAAATTTg                                                                                                                                 >  1:298263/1‑150 (MQ=255)
                gCATCGGTCACTGCCAGAATCGGCGGCGTATCAATAAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTg                                                                                                                                 >  2:1206900/1‑150 (MQ=255)
                 cATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTgg                                                                                                                                <  1:1338034/150‑1 (MQ=255)
                 cATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTgg                                                                                                                                <  2:550953/150‑1 (MQ=255)
                     ggTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAAt                                                                                                                            >  2:120655/1‑150 (MQ=255)
                       tCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATag                                                                                                                          <  2:837100/150‑1 (MQ=255)
                        cACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATaga                                                                                                                         <  2:1424902/150‑1 (MQ=255)
                               aGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGt                                                                                                                  <  1:1427091/150‑1 (MQ=255)
                                gAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGtt                                                                                                                 <  2:117654/150‑1 (MQ=255)
                                 aaTCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGttt                                                                                                                <  2:1279925/150‑1 (MQ=255)
                                    cggcggCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTagc                                                                                                             <  1:471970/150‑1 (MQ=255)
                                               caatcaaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTaa                                                                                                  <  1:79668/150‑1 (MQ=255)
                                                 atcaaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAATGCGTTCGCTGATCTACAGTTCAGAAGGTTTTGGCGGGACCTGACCGCGCGGGATCAGGTGAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAata                                                                                                <  1:1495878/150‑1 (MQ=255)
                                                    aaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCg                                                                                             <  1:1461625/150‑1 (MQ=255)
                                                     acacCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGc                                                                                            >  2:604501/1‑150 (MQ=255)
                                                       acCATGTCATAGTTTTTACTCGGCCAGTTCACCATTTCGGCAAAGCGTTCGCTCATCAACACTTCAGAATGATTTGGCGGCACCTGACCGCGGGGGATCAGGTCAAATTTGGCAATAGAGGGCGGTTTAGCAGCGGTAGTAATATCGccc                                                                                          <  2:834072/150‑1 (MQ=255)
                                                          aGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTgc                                                                                       >  2:237623/1‑149 (MQ=255)
                                                          aGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGa                                                                                       >  2:945305/1‑150 (MQ=255)
                                                           ggTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGAc                                                                                      <  1:1253731/150‑1 (MQ=255)
                                                           ggTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGACCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGAc                                                                                      >  1:800492/1‑150 (MQ=255)
                                                            gTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGAcc                                                                                     <  2:863928/150‑1 (MQ=255)
                                                                    ttttACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGaa                                                                             <  1:1321074/150‑1 (MQ=255)
                                                                         cTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTcc                                                                        <  2:176214/150‑1 (MQ=255)
                                                                         cTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTcc                                                                        >  2:873905/1‑150 (MQ=255)
                                                                         cTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCAAGAACAGTTCAGAAGGATTTGGCGGCACATGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCGGCGGTAGAAATATCGCCCTGACCAATCAGAATTTcc                                                                        >  2:1318746/1‑150 (MQ=255)
                                                                           cGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGgtagt                                                                                                    <  2:799338/120‑1 (MQ=255)
                                                                           cGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGgtagt                                                                                                    >  1:799338/1‑120 (MQ=255)
                                                                           cGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGa                                                                      >  1:801649/1‑150 (MQ=255)
                                                                            gCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGAc                                                                     >  1:772906/1‑150 (MQ=255)
                                                                                 gTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGcc                                                                <  2:801649/150‑1 (MQ=255)
                                                                                    cACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCAtt                                                             >  1:907673/1‑150 (MQ=255)
                                                                                      ccAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTaa                                                           >  1:1214068/1‑150 (MQ=255)
                                                                                       cAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAAc                                                          >  1:622360/1‑150 (MQ=255)
                                                                                         gtacGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGt                                                        >  1:510017/4‑150 (MQ=255)
                                                                                          ttCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGtt                                                       <  2:273502/150‑1 (MQ=255)
                                                                                             ggCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGttatt                                                    <  2:1440975/150‑1 (MQ=255)
                                                                                                   gCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGGCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGcc                                              <  2:257315/150‑1 (MQ=255)
                                                                                                    cGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGccc                                             <  1:958529/150‑1 (MQ=255)
                                                                                                       tCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAAc                                          <  1:497141/150‑1 (MQ=255)
                                                                                                             atgaaCAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGAAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCg                                    >  2:412696/4‑150 (MQ=255)
                                                                                                               caACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTg                                  >  1:174533/1‑150 (MQ=255)
                                                                                                                 aCAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTggg                                <  1:1351679/150‑1 (MQ=255)
                                                                                                                 aCAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTggg                                <  2:398432/150‑1 (MQ=255)
                                                                                                                    gTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGgtgt                             >  2:1467958/1‑150 (MQ=255)
                                                                                                                          aaGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCtt                       >  2:1512912/1‑150 (MQ=255)
                                                                                                                           aGGATTTGGAGCCACCTGACCGCGCGGGATCAGGTCAAATTTGGAAATAGAGGTCGGTTTAGCAGCGGTACTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCttt                      >  2:362595/1‑150 (MQ=255)
                                                                                                                            ggATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTg                     >  1:1000288/1‑150 (MQ=255)
                                                                                                                                 tGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCat                >  2:647903/1‑150 (MQ=255)
                                                                                                                                    cggcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCAAACAGCTCGTGGGTGTAGCCTTTGCGCATATc             >  1:113615/1‑150 (MQ=255)
                                                                                                                                    cggcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATAGACCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTCGGTTTACACTTTGCGCatata             >  2:886267/1‑149 (MQ=255)
                                                                                                                                     ttgaCCTGATCCCGCGCGGTCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTgcgc                  <  1:1332342/141‑1 (MQ=255)
                                                                                                                                     ttgaCCTGATCCCGCGCGGTCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTgcgc                  >  2:1332342/4‑144 (MQ=255)
                                                                                                                                     ggcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCg            >  1:179951/1‑150 (MQ=255)
                                                                                                                                      tgaCCTGATCCCGCGCGGTCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGc           >  1:262027/3‑150 (MQ=255)
                                                                                                                                      gcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGc           <  1:308454/150‑1 (MQ=255)
                                                                                                                                      gcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGc           <  1:647903/150‑1 (MQ=255)
                                                                                                                                      gcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGc           <  1:1467958/150‑1 (MQ=255)
                                                                                                                                              ccGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATc   <  2:772906/150‑1 (MQ=255)
                                                                                                                                               cgcgcgGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATca  >  1:903056/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCA  >  NC_000913/2133682‑2133974

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