Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,086,966 A→C 6.7% L294R (CTG→CGG) ‡ pgaC ← biofilm PGA synthase PgaCD, catalytic subunit; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,086,9660AC 6.7% 313.5 / 12.6 89L294R (CTG→CGG) ‡pgaCbiofilm PGA synthase PgaCD, catalytic subunit; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine synthase
Reads supporting (aligned to +/- strand):  ref base A (43/40);  new base C (3/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

GACCAGGCAGGTGAAGGCCCATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGTGATCATATCGTCACTCCAGTAACCCAC  >  NC_000913/1086833‑1087114
                                                                                                                                     |                                                                                                                                                    
gACCAGGCAGGTGAAGGCCCATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCcaca                                                                                                                                       >  4:354187/1‑149 (MQ=255)
    aggcaggTGAAGGCCCATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGccc                                                                                                                                   <  3:354187/149‑1 (MQ=255)
     ggcaggTGAAGGCCCATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCt                                                                                                                                  >  4:1421179/1‑149 (MQ=255)
         ggTGAAGGCCCATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTa                                                                                                                              <  3:1112759/149‑1 (MQ=255)
                 cccATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCa                                                                                                                      >  3:1191849/1‑149 (MQ=255)
                  ccATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGccc                                                                                                                                   <  4:979904/135‑1 (MQ=255)
                  ccATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGccc                                                                                                                                   >  3:979904/1‑135 (MQ=255)
                  ccATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGc                                                                                                                                     <  1:119752/133‑1 (MQ=255)
                  ccATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGc                                                                                                                                     >  4:673978/1‑133 (MQ=255)
                  ccATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGc                                                                                                                                     >  2:119752/1‑133 (MQ=255)
                  ccATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGc                                                                                                                                     <  3:673978/133‑1 (MQ=255)
                  ccATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAg                                                                                                                     >  1:735309/1‑149 (MQ=255)
                    atatTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGc                                                                                                                   >  1:708184/1‑149 (MQ=255)
                     tatTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGttt                                                                                                                        <  1:917232/143‑1 (MQ=255)
                     tatTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGttt                                                                                                                        >  2:917232/1‑143 (MQ=255)
                       ttGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCAtt                                                                                                                <  2:735309/149‑1 (MQ=255)
                             tcAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATc                                                                                                          <  4:1191849/149‑1 (MQ=255)
                               aGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCa                                                                                                        <  3:1421179/149‑1 (MQ=255)
                                ggCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTc                                                                                                                       >  2:1272354/1‑133 (MQ=255)
                                ggCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTc                                                                                                                       <  1:1272354/133‑1 (MQ=255)
                                 gCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTaa                                                                                                                             <  1:454547/126‑1 (MQ=255)
                                 gCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTaa                                                                                                                             >  2:454547/1‑126 (MQ=255)
                                 gCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAAc                                                                                                                            >  1:1089063/1‑127 (MQ=255)
                                 gCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAAc                                                                                                                            <  2:1089063/127‑1 (MQ=255)
                                 gCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTc                                                                                                                       >  2:807529/1‑132 (MQ=255)
                                 gCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTc                                                                                                                       <  1:807529/132‑1 (MQ=255)
                                 gCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGc                                                                                                      >  4:1217715/1‑149 (MQ=255)
                                         cAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCa                                                                                                                      <  3:27874/125‑1 (MQ=255)
                                         cAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCa                                                                                                                      >  4:27874/1‑125 (MQ=255)
                                           aaaaaCAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCg                                                                                            >  3:239702/1‑149 (MQ=255)
                                             aaaCAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCg                                                                                            <  3:431259/147‑1 (MQ=255)
                                             aaaCAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCg                                                                                            >  4:431259/1‑147 (MQ=255)
                                              aaCAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTgg                                                                                         <  2:708184/149‑1 (MQ=255)
                                                    ggCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCa                                                                                                        >  3:187329/1‑128 (MQ=255)
                                                    ggCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCa                                                                                                        <  4:187329/128‑1 (MQ=255)
                                                          aTTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATc                                                                                                          <  3:21014/120‑1 (MQ=255)
                                                          aTTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATc                                                                                                          >  4:21014/1‑120 (MQ=255)
                                                           ttCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATc                                                                            <  4:239702/149‑1 (MQ=255)
                                                              gAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTGATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTc                                                                         >  4:125139/1‑149 (MQ=255)
                                                                 aGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACGTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTGATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCAc                                                                      <  3:125139/149‑1 (MQ=255)
                                                                   ttttCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTg                                                                    >  1:1334599/1‑149 (MQ=255)
                                                                        tttGCGCCACAACCTTGTCGTATTTTTGAGGAATACTTCTGCGCCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCa                                                               >  3:866147/1‑149 (MQ=255)
                                                                        tttGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCGGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTg                                                                    <  3:1251660/144‑1 (MQ=255)
                                                                        tttGCGCCACAACCTTGTCATATTTTCGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTg                                                                    >  4:1251660/1‑144 (MQ=255)
                                                                          tGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCGCCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAAc                                                             <  4:866147/149‑1 (MQ=255)
                                                                              ccACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCa                                                         <  3:873023/149‑1 (MQ=255)
                                                                              ccACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCa                                                         <  3:1217715/149‑1 (MQ=255)
                                                                              ccACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCa                                                         >  4:873023/1‑149 (MQ=255)
                                                                                  aaCCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTc                                                                                                                       >  4:55552/1‑83 (MQ=255)
                                                                                  aaCCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTc                                                                                                                       <  3:55552/83‑1 (MQ=255)
                                                                                     cTTGTCATATTTTTGTGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTaaaaaa                                                                              <  4:419841/121‑1 (MQ=255)
                                                                                     cTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTaaaaaa                                                                              >  3:419841/1‑121 (MQ=255)
                                                                                            tatTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAat                                                                                                             >  3:1303750/1‑83 (MQ=255)
                                                                                            tatTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAat                                                                                                             <  4:1303750/83‑1 (MQ=255)
                                                                                            tatTTTTGAGGAATACTTCTGCACAGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATCCAACTGCAGCTTCCAGCTaata                                           >  3:809121/1‑149 (MQ=255)
                                                                                             atTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTaatat                                          <  2:1334599/149‑1 (MQ=255)
                                                                                                tttGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTaatatcaa                                       >  4:186845/1‑149 (MQ=255)
                                                                                                      gAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGccc                                                                                             <  4:868126/89‑1 (MQ=255)
                                                                                                      gAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGccc                                                                                             >  3:868126/1‑89 (MQ=255)
                                                                                                          aCTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAata                                                                                                            >  1:1014801/1‑70 (MQ=255)
                                                                                                          aCTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAata                                                                                                            <  2:1014801/70‑1 (MQ=255)
                                                                                                          aCTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACt                                                                     <  1:1278826/109‑1 (MQ=255)
                                                                                                          aCTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACt                                                                     >  2:1278826/1‑109 (MQ=255)
                                                                                                            ttCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGttt                                                                                                                        <  4:1006107/56‑1 (MQ=255)
                                                                                                            ttCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGttt                                                                                                                        >  3:1006107/1‑56 (MQ=255)
                                                                                                            ttCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCa                                                                                                        >  3:1158811/1‑72 (MQ=255)
                                                                                                            ttCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCa                                                                                                        <  4:1158811/72‑1 (MQ=255)
                                                                                                             tCTGCACCGCCCTGAGCCAAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGTGAt                          <  4:809121/149‑1 (MQ=255)
                                                                                                                 cACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGTGATCata                      >  3:1414721/1‑149 (MQ=255)
                                                                                                                  aCCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTc                                >  1:990414/1‑138 (MQ=255)
                                                                                                                  aCCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTc                                <  2:990414/138‑1 (MQ=255)
                                                                                                                   ccGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGt                             >  4:283806/1‑140 (MQ=255)
                                                                                                                   ccGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACAGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGt                             <  3:283806/140‑1 (MQ=255)
                                                                                                                    cGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCtt                                                     <  2:542692/115‑1 (MQ=255)
                                                                                                                    cGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCtt                                                     >  1:542692/1‑115 (MQ=255)
                                                                                                                    cGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTcc                                                   <  1:128886/117‑1 (MQ=255)
                                                                                                                    cGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTcc                                                   >  2:128886/1‑117 (MQ=255)
                                                                                                                        cTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTc                                                    <  1:129426/112‑1 (MQ=255)
                                                                                                                        cTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTc                                                    >  2:129426/1‑112 (MQ=255)
                                                                                                                         tGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGTGATCATATCGTCACt              >  2:1280439/1‑149 (MQ=255)
                                                                                                                          gAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCaa                                                              <  3:695410/100‑1 (MQ=255)
                                                                                                                          gAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCaa                                                              >  4:695410/1‑100 (MQ=255)
                                                                                                                             aacAGCGCCTGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTaaaaaa                                                                              >  4:1229468/3‑81 (MQ=255)
                                                                                                                             aacAGCGCCTGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTaaaaaa                                                                              <  3:1229468/79‑1 (MQ=255)
                                                                                                                             aacAGCGCCTGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTaatatcaa                                       >  4:565886/3‑120 (MQ=255)
                                                                                                                             aacAGCGCCTGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTaatatcaa                                       <  3:565886/118‑1 (MQ=255)
                                                                                                                             aacAGCGCCTGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTaa                                             <  3:1290503/112‑1 (MQ=255)
                                                                                                                             aacAGCGCCTGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGTGATCatat                     >  4:1334598/3‑138 (MQ=255)
                                                                                                                              ccAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGTGATCATATCGTCACTCCAGt         >  4:1140118/1‑149 (MQ=255)
                                                                                                                                     aGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGTGATCATATCGTCACTCCAGTAACCCAc  >  3:1072396/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
GACCAGGCAGGTGAAGGCCCATATTGTCGTCAGGCAGTATTCAAAAAACAGCGGCCACATTCGAAAGTTTTCTTTGCGCCACAACCTTGTCATATTTTTGAGGAATACTTCTGCACCGCCCTGAGCCCAGCGCAGGCGCTGTTTCCACAGCCCTTTTAACGTTTCAGGCATTAATATCCAGCACAGTGCCCGTGGCTCGTAAAAAATCGTCCACTGATTCAACTGCAGCTTCCAGCTAATATCAATATCTTCGGTGATCATATCGTCACTCCAGTAACCCAC  >  NC_000913/1086833‑1087114

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