Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 268,395 A→G V189A (GTC→GCC)  ykfA ← CP4‑6 prophage; putative GTP‑binding protein YkfA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913268,3950AG98.7% 258.5 / ‑5.3 76V189A (GTC→GCC) ykfACP4‑6 prophage; putative GTP‑binding protein YkfA
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (40/35);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGT  >  NC_000913/268260‑268530
                                                                                                                                       |                                                                                                                                       
aCCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCggg                                                                                                                                     >  2:1135489/1‑140 (MQ=255)
     cgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGa                                                                                                                                <  1:472485/140‑1 (MQ=255)
      gtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGaa                                                                                                                               >  1:398703/1‑140 (MQ=255)
       tGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAAc                                                                                                                              >  2:5059/1‑140 (MQ=255)
       tGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAAc                                                                                                                              >  1:678650/1‑140 (MQ=255)
               gaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATg                                                                                                                      <  2:547021/140‑1 (MQ=255)
                 aggaaggaGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTgg                                                                                                                                 <  1:1052136/127‑1 (MQ=255)
                 aggaaggaGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTgg                                                                                                                                 >  2:1052136/1‑127 (MQ=255)
                        aGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAgcaggca                                                                                                             >  1:170761/1‑140 (MQ=255)
                               acacCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTg                                                                                                      <  2:171121/140‑1 (MQ=255)
                                  ccaccGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCgg                                                                                                   >  2:1119123/1‑140 (MQ=255)
                                   caccGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGc                                                                                                  >  2:1139966/1‑140 (MQ=255)
                                   caccGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGc                                                                                                  <  2:551727/140‑1 (MQ=255)
                                     ccGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGAAACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTg                                                                                                >  2:500410/1‑140 (MQ=255)
                                          cTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATAtt                                                                                           <  1:547168/140‑1 (MQ=255)
                                             gtggCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATAttctg                                                                                        >  2:944824/1‑140 (MQ=255)
                                             gtggCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATAttctg                                                                                        >  1:477430/1‑140 (MQ=255)
                                                             ggcaCTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGAc                                                                        >  1:283566/1‑140 (MQ=255)
                                                               caCTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGTGGCTGATATTCTGCTTCTGCGCCGTGGACaa                                                                      >  1:822787/1‑140 (MQ=255)
                                                               caCTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACaa                                                                      <  2:398703/140‑1 (MQ=255)
                                                               caCTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACaa                                                                      >  1:996187/1‑140 (MQ=255)
                                                                    aaTCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTc                                                                 >  2:83729/1‑140 (MQ=255)
                                                                    aaTCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTc                                                                 >  2:616037/1‑140 (MQ=255)
                                                                       catcCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcccc                                                              <  2:170761/140‑1 (MQ=255)
                                                                         tcCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCa                                                            >  2:1078933/1‑140 (MQ=255)
                                                                            gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTg                                                                   >  2:148103/1‑130 (MQ=255)
                                                                            gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTg                                                                   <  1:148103/130‑1 (MQ=255)
                                                                            gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccc                                                               <  2:187149/134‑1 (MQ=255)
                                                                            gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccc                                                               >  1:187149/1‑134 (MQ=255)
                                                                              tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCt                                                       >  1:577301/1‑140 (MQ=255)
                                                                                ggCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTgg                                                     >  2:272557/1‑140 (MQ=255)
                                                                                ggCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTgg                                                     >  1:306003/1‑140 (MQ=255)
                                                                                 gCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGt                                                    <  1:1136402/140‑1 (MQ=255)
                                                                                         cTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCgg                                            <  1:1135489/140‑1 (MQ=255)
                                                                                         cTTAAGCCCCCATTGTAGACGGACCGAAACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGAGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCgg                                            <  1:500410/140‑1 (MQ=255)
                                                                                          ttAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGc                                           >  2:252552/1‑140 (MQ=255)
                                                                                            aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCt                                         >  1:14416/1‑140 (MQ=255)
                                                                                             aGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTc                                              >  1:869395/1‑134 (MQ=255)
                                                                                             aGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTc                                              <  2:869395/134‑1 (MQ=255)
                                                                                               cccccATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTAt                                      >  2:441978/1‑140 (MQ=255)
                                                                                                ccccATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATc                                     >  2:243538/1‑140 (MQ=255)
                                                                                                 cccATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATcc                                    >  1:303720/1‑140 (MQ=255)
                                                                                                   cATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGa                                  <  1:272557/140‑1 (MQ=255)
                                                                                                   cATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGa                                  >  2:51398/1‑140 (MQ=255)
                                                                                                   cATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGa                                  <  1:638807/140‑1 (MQ=255)
                                                                                                   cATTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGa                                  >  1:339820/1‑140 (MQ=255)
                                                                                                    aTTGTAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc                                 <  2:283566/140‑1 (MQ=255)
                                                                                                        tAGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGc                             <  2:577301/140‑1 (MQ=255)
                                                                                                         aGACGGACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCt                            <  2:477430/140‑1 (MQ=255)
                                                                                                             ggACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATa                        <  1:51398/140‑1 (MQ=255)
                                                                                                              gACCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaa                       <  2:959124/140‑1 (MQ=255)
                                                                                                               aCCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaac                      >  2:661081/1‑140 (MQ=255)
                                                                                                               aCCGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaac                      >  2:666091/1‑140 (MQ=255)
                                                                                                                ccGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGTTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaaca                     >  1:239094/1‑140 (MQ=255)
                                                                                                                 cGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCg                                   <  2:64590/125‑1 (MQ=255)
                                                                                                                 cGACACAGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCg                                   >  1:64590/1‑125 (MQ=255)
                                                                                                                    cacaGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaacaaaca                 <  2:303720/140‑1 (MQ=255)
                                                                                                                    cacaGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaacaaaca                 <  2:306003/140‑1 (MQ=255)
                                                                                                                     acaGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAg                <  1:661081/140‑1 (MQ=255)
                                                                                                                       aGCACACACCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCa              <  1:1078933/140‑1 (MQ=255)
                                                                                                                           cacacCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCtt          >  2:346235/1‑140 (MQ=255)
                                                                                                                            acacCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTa         <  1:5059/140‑1 (MQ=255)
                                                                                                                             cacCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt        <  1:944824/140‑1 (MQ=255)
                                                                                                                             cacCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt        <  1:83729/140‑1 (MQ=255)
                                                                                                                             cacCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt        <  1:616037/140‑1 (MQ=255)
                                                                                                                             cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt        <  2:996187/140‑1 (MQ=255)
                                                                                                                              acCGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATg       >  2:1132095/1‑140 (MQ=255)
                                                                                                                               ccGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGTGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGc      <  2:822787/140‑1 (MQ=255)
                                                                                                                               ccGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGc      >  2:152972/1‑140 (MQ=255)
                                                                                                                               ccGGATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGc      <  1:1132095/140‑1 (MQ=255)
                                                                                                                                 ggATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCg    >  1:838984/1‑140 (MQ=255)
                                                                                                                                  gATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGTTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCgg   <  2:239094/140‑1 (MQ=255)
                                                                                                                                  gATGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCgg   <  2:14416/140‑1 (MQ=255)
                                                                                                                                   aTGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGt  >  1:359494/1‑140 (MQ=255)
                                                                                                                                    tGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGcc     <  2:19672/136‑1 (MQ=255)
                                                                                                                                    tGGGCGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGcc     >  1:19672/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                       
ACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGT  >  NC_000913/268260‑268530

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