Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,133,826 G→C 5.9% 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 5.9% 265.1 / 4.8 79R609R (CGC→CGG) wzcprotein‑tyrosine kinase Wzc
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base C (3/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAG  >  NC_000913/2133680‑2133967
                                                                                                                                                  |                                                                                                                                             
tGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACcgcgc                                                                                                                                            >  1:937377/1‑150 (MQ=255)
    gacCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGa                                                                                                                                        <  2:1845627/150‑1 (MQ=255)
    gacCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGa                                                                                                                                        >  1:24338/1‑150 (MQ=255)
     acCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGAt                                                                                                                                       >  1:900239/1‑150 (MQ=255)
      ccAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATc                                                                                                                                      >  1:352736/1‑150 (MQ=255)
       caacaaTTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCa                                                                                                                                     <  1:1299712/150‑1 (MQ=255)
       caacaaTTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTAGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCa                                                                                                                                     >  2:1101233/1‑150 (MQ=255)
          caaTTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCATTTCGGTAAAGCGTTGGCTCATAAACAGTTCTTAAGGATTTGTCTGCACCTGACCGCGCGGGTTCAGGt                                                                                                                                  <  2:782877/150‑1 (MQ=255)
               gCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGATCCCGCGCGGTCAGGtgccgc                                                                                                                             <  1:1417246/150‑6 (MQ=255)
                ccGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACGTGATCCCGCGCGGTCAGGtgccgcc                                                                                                                            <  2:473721/150‑7 (MQ=255)
                ccGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAgg                                                                                                                                   >  2:1640298/1‑143 (MQ=255)
                ccGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAgg                                                                                                                                   <  1:1640298/143‑1 (MQ=255)
                ccGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAAtt                                                                                                                            >  1:1664778/1‑150 (MQ=255)
                ccGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAAtt                                                                                                                            >  2:1862876/1‑150 (MQ=255)
                 cGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTTACCGCGCGGGATCAGTTCAAAttt                                                                                                                           <  2:1450878/150‑1 (MQ=255)
                    aTCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGc                                                                                                                        <  2:21826/150‑1 (MQ=255)
                    aTCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGc                                                                                                                        <  1:1776739/150‑1 (MQ=255)
                     tCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCa                                                                                                                       <  1:1510622/150‑1 (MQ=255)
                     tCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCa                                                                                                                       <  1:1503518/150‑1 (MQ=255)
                      cGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCaa                                                                                                                      <  1:1355912/150‑1 (MQ=255)
                           aCTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATagag                                                                                                                 <  2:934697/150‑1 (MQ=255)
                              gCCAGAATCGGCGTCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTc                                                                                                              >  1:1786173/1‑150 (MQ=255)
                              gCCAGAATCGGCGGCGTATCAATCAACTCCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTc                                                                                                              >  2:626673/1‑150 (MQ=255)
                              gCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCACAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGtg                                                                                                              >  2:187597/1‑149 (MQ=255)
                               ccAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCg                                                                                                             >  1:41958/1‑150 (MQ=255)
                                 aGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGt                                                                                                           <  2:1070926/150‑1 (MQ=255)
                                   aaTCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGttt                                                                                                         <  2:765572/150‑1 (MQ=255)
                                      cggcggCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTagc                                                                                                      <  1:804573/150‑1 (MQ=255)
                                       ggcggcGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTagca                                                                                                     >  1:932112/1‑150 (MQ=255)
                                           gcgTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCgg                                                                                                 <  2:443822/150‑1 (MQ=255)
                                            cgTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGgt                                                                                                >  2:526644/1‑150 (MQ=255)
                                            cgTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGgt                                                                                                <  1:526644/150‑1 (MQ=255)
                                            cgTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGgt                                                                                                <  1:1767363/150‑1 (MQ=255)
                                                  aatcaaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAat                                                                                          >  2:1810068/1‑150 (MQ=255)
                                                  aatcaaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAat                                                                                          >  2:285217/1‑150 (MQ=255)
                                                     caaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATc                                                                                       >  1:1493448/1‑150 (MQ=255)
                                                      aaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATct                                                                                      <  2:932112/150‑2 (MQ=255)
                                                           cAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTg                                                                                 <  2:1386207/150‑1 (MQ=255)
                                                             ggTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGAc                                                                               <  1:285217/150‑1 (MQ=255)
                                                              gTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGAcc                                                                              >  2:263632/1‑150 (MQ=255)
                                                                cATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCaa                                                                            >  2:1679635/1‑150 (MQ=255)
                                                                cATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCaa                                                                            >  2:1664958/1‑150 (MQ=255)
                                                                   agcttttACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCa                                                                         <  2:1596608/147‑1 (MQ=255)
                                                                    gTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAg                                                                        <  2:352736/150‑1 (MQ=255)
                                                                    gTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAg                                                                        <  1:1467476/150‑1 (MQ=255)
                                                                      ttttACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGaa                                                                      >  1:330212/1‑150 (MQ=255)
                                                                           cTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATc                                                                                                                                      <  1:461972/81‑1 (MQ=255)
                                                                           cTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATc                                                                                                                                      >  2:461972/1‑81 (MQ=255)
                                                                            tCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGATGTCGGGTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCg                                                                <  1:754558/150‑1 (MQ=255)
                                                                            tCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCg                                                                <  1:434496/150‑1 (MQ=255)
                                                                             cGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGa                                                               >  1:1155170/1‑150 (MQ=255)
                                                                                  aGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGc                                                          <  2:261367/150‑1 (MQ=255)
                                                                                    ttCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCa                                                        <  1:1101233/150‑1 (MQ=255)
                                                                                      cACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGcc                                                         >  2:8517/1‑147 (MQ=255)
                                                                                      cACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGcc                                                         <  1:8517/147‑1 (MQ=255)
                                                                                       aCCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTa                                                     >  2:1214579/1‑150 (MQ=255)
                                                                                       aCCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTa                                                     >  1:682411/1‑150 (MQ=255)
                                                                                           gTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGt                                                 >  2:351437/1‑150 (MQ=255)
                                                                                              cGGCAAAGCGTGCGCGCATCAACCGTGCAGAACGAGTTGGCGGCACCTGACCGCGCGGGCGCAGGTCAAATTTGGCACTAGAGCGCGGCTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGttat                                              <  2:1523369/150‑1 (MQ=255)
                                                                                                    aGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGc                                        <  2:900239/150‑1 (MQ=255)
                                                                                                       gTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACCTTATTAGTGCCCa                                     <  2:24338/150‑1 (MQ=255)
                                                                                                          cGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACa                                  <  1:1664958/150‑1 (MQ=255)
                                                                                                              catcaACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTc                              >  2:1132918/1‑150 (MQ=255)
                                                                                                              catcaACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTc                              >  2:342516/1‑150 (MQ=255)
                                                                                                                tcaACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGt                            >  1:285607/1‑150 (MQ=255)
                                                                                                                tcaACAGTTCAGAAGGATTTGGCGACACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGt                            >  2:912822/1‑150 (MQ=255)
                                                                                                                       ttCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTa                     >  1:1503772/1‑150 (MQ=255)
                                                                                                                        taaGAAGTATTTGGAGGCACCTGACACCGCTGGATAAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAg                    >  2:1921200/3‑150 (MQ=255)
                                                                                                                         cAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGc                   >  2:1915814/1‑150 (MQ=255)
                                                                                                                         cAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGCTGTAGc                   >  2:559382/1‑150 (MQ=255)
                                                                                                                          aGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGcc                  <  1:1679635/150‑1 (MQ=255)
                                                                                                                            aaGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCtt                >  2:742241/1‑150 (MQ=255)
                                                                                                                                aTTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTgcg            >  2:1743531/1‑150 (MQ=255)
                                                                                                                                 tttGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTCgcgc           <  2:41958/150‑1 (MQ=255)
                                                                                                                                    ggcggcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCata        <  1:1743531/150‑1 (MQ=255)
                                                                                                                                    aatttgaCCTGATCCCGCGCGGTCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCata        >  1:493893/7‑150 (MQ=255)
                                                                                                                                     gcggcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCatat       >  1:735925/1‑150 (MQ=255)
                                                                                                                                         gaCCTGATCCCGCGCGGTCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCa   >  1:1283193/2‑150 (MQ=255)
                                                                                                                                         cACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCa   >  2:1981799/1‑150 (MQ=255)
                                                                                                                                          aCCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAg  <  1:1707366/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAG  >  NC_000913/2133680‑2133967

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