Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 266,465 G→A 5.0% A30V (GCC→GTC)  yafY ← lipoprotein, inner membrane; degP regulator; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,4650GA 5.0% 619.8 / 14.5 179A30V (GCC→GTC) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (81/89);  new base A (4/5);  total (85/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 186 total aligned reads displayed.

CAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAAT  >  NC_000913/266318‑266612
                                                                                                                                                   |                                                                                                                                                   
cAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGc                                                                                                                                                    <  4:901855/149‑1 (MQ=255)
 aGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGcc                                                                                                                                                   <  3:795704/149‑1 (MQ=255)
  gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGGCCTTAAATTTACTGATGGCCt                                                                                                                                                  >  3:1410001/1‑149 (MQ=255)
     aTCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCa                                                                                                                                               <  4:1298497/149‑1 (MQ=255)
     aTCGGGCGACTGCCAGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCa                                                                                                                                               >  3:105010/1‑149 (MQ=255)
        gggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGt                                                                                                                                            >  2:811216/1‑149 (MQ=255)
            gACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatc                                                                                                                                        >  2:614184/1‑149 (MQ=255)
                  cGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt                                                                                                                                  <  1:811216/149‑1 (MQ=255)
                         gTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatca                                                                                                                                       >  2:733994/1‑137 (MQ=255)
                         gTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCg                                                                                                                           >  3:635472/1‑149 (MQ=255)
                         gTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCg                                                                                                                           >  1:552285/1‑149 (MQ=255)
                          tAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGc                                                                                                                          >  3:578159/1‑149 (MQ=255)
                           aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                                                                                                                         <  1:541944/149‑1 (MQ=255)
                                 tGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt                                                                                                                   <  4:578159/149‑1 (MQ=255)
                                 tGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt                                                                                                                   <  4:1410001/149‑1 (MQ=255)
                                   cAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATga                                                                                                                 <  2:552285/149‑1 (MQ=255)
                                     aGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgaga                                                                                                               >  1:39541/1‑149 (MQ=255)
                                        tGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgacagaa                                                                                                            <  3:1152154/149‑6 (MQ=255)
                                           gACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATga                                                                                                                 <  3:323597/141‑1 (MQ=255)
                                           gACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACAGCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATga                                                                                                                 >  2:1333496/1‑141 (MQ=255)
                                            aCCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc                                                                                                            >  2:668077/1‑145 (MQ=255)
                                            aCCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccatcagt                                                                                                                                     >  2:188400/1‑118 (MQ=18)
                                            aCCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccatcagt                                                                                                                                     <  1:188400/120‑3 (MQ=18)
                                               ggATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTc                                                                                                                               >  2:77940/1‑123 (MQ=255)
                                                gATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgac                                                                                                                >  1:1144556/1‑136 (MQ=255)
                                                gATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATga                                                                                                                 >  3:1341783/1‑136 (MQ=255)
                                                gATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggt                                                                                                    <  4:467888/149‑1 (MQ=255)
                                                 aTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggtg                                                                                                   <  4:105010/149‑1 (MQ=255)
                                                   aaGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg                                                                                                                                 <  3:1360651/117‑1 (MQ=255)
                                                     gtagtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATgg                                                                                                                                                     >  1:1282590/1‑95 (MQ=255)
                                                     gtagtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCa                                                                                                                                               <  1:1427137/101‑1 (MQ=255)
                                                     gtagtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACtt                                                                                                                                    <  1:192727/112‑1 (MQ=255)
                                                     gtagtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCaa                                                                                               <  4:1439931/149‑1 (MQ=255)
                                                       agtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccatca                                                                                                                                       >  1:703471/1‑107 (MQ=18)
                                                        gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat                                                                                                                                         <  3:1316853/104‑1 (MQ=255)
                                                        gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat                                                                                                                                         >  4:1316853/1‑104 (MQ=255)
                                                        gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccatcagta                                                                                                                                    >  1:872831/1‑106 (MQ=17)
                                                         tCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatca                                                                                                                                       >  4:845994/1‑105 (MQ=255)
                                                         tCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggt                                                                                                    <  3:953888/140‑1 (MQ=255)
                                                         tCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggt                                                                                                    >  4:953888/1‑140 (MQ=255)
                                                           aaCCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAgcg                                                                                         >  4:936456/1‑149 (MQ=255)
                                                           aaCCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccatcagtaaa                                                                                                                                  >  2:229784/1‑103 (MQ=17)
                                                              cTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTGAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgac                                                                                                                <  1:626665/123‑2 (MQ=255)
                                                               tgtgACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATga                                                                                                                 >  1:72205/1‑121 (MQ=255)
                                                               tgtgACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcag                                                                                     <  3:936456/149‑1 (MQ=255)
                                                               tgtgACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACACAGTGGTGGCAACCAGCGCcag                                                                                     >  4:1381400/1‑149 (MQ=255)
                                                                   aCCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAg                                                                                                                                              >  2:1466344/1‑88 (MQ=255)
                                                                     cATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATga                                                                                                                 <  2:1353636/115‑1 (MQ=255)
                                                                      aTTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTc                                                                                                                                                >  3:385406/1‑83 (MQ=255)
                                                                      aTTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGc                                                                                                 <  1:836149/130‑1 (MQ=255)
                                                                          ttcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCAc                                                                                                                                      >  4:1456954/1‑89 (MQ=255)
                                                                            cGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgacag                                                                                                              >  1:1314066/1‑108 (MQ=255)
                                                                               tGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAgcgc                                                                                        <  3:1169270/130‑1 (MQ=255)
                                                                                gACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACtt                                                                                                                                    >  2:583580/1‑85 (MQ=255)
                                                                                 aCATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTa                                                                     >  2:1269854/1‑147 (MQ=255)
                                                                                   atcatcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGtt                                                                 <  2:39541/149‑1 (MQ=255)
                                                                                    tcatcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt                                                                                                                   >  4:454019/1‑98 (MQ=255)
                                                                                    tcatcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaaa                                                                                                             >  3:505732/1‑104 (MQ=255)
                                                                                      atcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCa                                                              <  4:635472/149‑1 (MQ=255)
                                                                                       tcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgtgt                                                                         >  3:1014348/1‑137 (MQ=255)
                                                                                       tcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgtgt                                                                         <  4:1014348/137‑1 (MQ=255)
                                                                                       tcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCat                                                             >  1:459134/1‑149 (MQ=255)
                                                                                       tcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCat                                                             <  2:63537/149‑1 (MQ=255)
                                                                                         cTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACagag                                                                                                        <  1:662170/104‑1 (MQ=255)
                                                                                         cTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACagag                                                                                                        <  1:464702/104‑1 (MQ=255)
                                                                                         cTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACagag                                                                                                        >  2:464702/1‑104 (MQ=255)
                                                                                         cTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCatta                                                           >  1:391452/1‑149 (MQ=255)
                                                                                           ggtggCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCattatt                                                         <  3:1241282/149‑1 (MQ=255)
                                                                                             tggCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgt                                                                           <  2:113792/129‑1 (MQ=255)
                                                                                              ggCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTAttttt                                                      <  1:99130/149‑1 (MQ=255)
                                                                                              ggCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTAttttt                                                      <  4:157926/149‑1 (MQ=255)
                                                                                                cTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTAttttttt                                                    <  2:459134/149‑1 (MQ=255)
                                                                                                 tGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCatta                                                           <  2:989196/141‑1 (MQ=255)
                                                                                                 tGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCatta                                                           >  1:989196/1‑141 (MQ=255)
                                                                                                 tGACCACATCACTGAAACGCGGTGGGGTGAGGCTCTTAAATTTACTCATGGCCTTCAGGTCATCACTTCTGTCATCGCACGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTGTACGTTTCATTATTTTTTTc                                                   >  1:80122/1‑149 (MQ=18)
                                                                                                  gACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTcc                                                  <  3:146840/149‑1 (MQ=255)
                                                                                                   aCCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggtg                                                                                                   <  2:1295919/99‑1 (MQ=255)
                                                                                                     cacaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCattat                                                          >  3:227723/1‑138 (MQ=255)
                                                                                                     cacaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCa                                                              >  2:851128/1‑134 (MQ=255)
                                                                                                      acaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatc                                                                                                                                        >  3:1438672/1‑59 (MQ=255)
                                                                                                       caTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGaaa                                             <  1:597998/149‑1 (MQ=255)
                                                                                                        aTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAAt                                            <  1:1341635/149‑1 (MQ=255)
                                                                                                          cACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCa                                          >  1:1021779/1‑149 (MQ=255)
                                                                                                            cTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGa                                        <  4:1352010/149‑1 (MQ=255)
                                                                                                                 aCGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTcc                                                  <  4:360397/134‑1 (MQ=255)
                                                                                                                 aCGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAAc                                   <  3:1455385/149‑1 (MQ=255)
                                                                                                                             gAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCaaa                       <  2:907981/149‑1 (MQ=255)
                                                                                                                                gTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggtg                                                                                                   <  3:886634/70‑1 (MQ=255)
                                                                                                                                 tCCTTAAATTTACTGATGGCCTTCAGGTCGTCACTTCTGTCATCGCa                                                                                                                         >  1:109802/1‑47 (MQ=255)
                                                                                                                                    ttAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGtt                                                                        >  4:301194/1‑93 (MQ=255)
                                                                                                                                    ttAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCaaa                       >  1:641815/1‑142 (MQ=255)
                                                                                                                                    ttAAATTTACTGATGGCCTTCAGGTCATCACTTCAGTCATCGCATGCGATGAGAAACAGAGTTGTCGCTACCAGCGCCAGCAGAGGCAGTGCTTTACGTTTCATTATTTTTTTCATGAAATCAGACGAACCGCTTTGGCAAAAACATaa                <  2:80122/149‑1 (MQ=17)
                                                                                                                                      aaaTTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAAt              <  4:1153480/149‑1 (MQ=255)
                                                                                                                                        aTTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACaga                                                                                                         >  2:644863/1‑56 (MQ=255)
                                                                                                                                        aTTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTg                            >  1:743439/1‑133 (MQ=255)
                                                                                                                                         tttACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgt                                                                           <  1:277038/85‑1 (MQ=255)
                                                                                                                                         tttACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgt                                                                           >  2:277038/1‑85 (MQ=255)
                                                                                                                                              tGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGTAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTAttttt                                                      <  4:1235563/101‑1 (MQ=255)
                                                                                                                                                  ggCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTAc                                                                    >  2:714778/1‑83 (MQ=255)
                                                                                                                                                  ggCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTAc                                                                    <  1:714778/83‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAAT  >  NC_000913/266318‑266612

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: