Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,133,824 C→T 8.8% G610D (GGT→GAT)  wzc ← protein‑tyrosine kinase Wzc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,133,8240CT 8.8% 197.2 / 7.6 68G610D (GGT→GAT) wzcprotein‑tyrosine kinase Wzc
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base T (2/4);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

CGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATC  >  NC_000913/2133675‑2133973
                                                                                                                                                     |                                                                                                                                                     
cGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGAc                                                                                                                                                       >  2:1461625/1‑150 (MQ=255)
cGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGAc                                                                                                                                                       <  1:357759/150‑1 (MQ=255)
       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:151487/150‑5 (MQ=255)
                   tGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGATCCCGCGCGGTCAGGtgccg                                                                                                                                    <  1:578968/150‑5 (MQ=255)
                       gCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTg                                                                                                                                <  2:337593/150‑1 (MQ=255)
                       gCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTg                                                                                                                                >  2:1253731/1‑150 (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:873905/1‑150 (MQ=255)
                                                                                cTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTcc                                                                       <  2:176214/150‑1 (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                 >  2:1332342/4‑144 (MQ=255)
                                                                                                                                            ttgaCCTGATCCCGCGCGGTCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTgcgc                 <  1:1332342/141‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
CGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATC  >  NC_000913/2133675‑2133973

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: ‑