Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,383,959 A→T 7.1% intergenic (+35/+36) elaD → / ← yfbK protease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,383,9590AT 7.1% 348.1 / 11.4 122intergenic (+35/+36)elaD/yfbKprotease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (54/59);  new base T (3/6);  total (57/65)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 127 total aligned reads displayed.

AGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTG  >  NC_000913/2383811‑2384108
                                                                                                                                                    |                                                                                                                                                     
aGAAATTGTTATTTAGATAAATATGGCGATTCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACAATCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAatat                                                                                                                                                      <  1:1095219/150‑1 (MQ=255)
   aaTTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCt                                                                                                                                                   >  2:128696/1‑150 (MQ=255)
       gTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACt                                                                                                                                               <  2:902712/150‑1 (MQ=255)
        ttATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTg                                                                                                                                              <  1:354207/150‑1 (MQ=255)
         tATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTgg                                                                                                                                             <  2:1613794/150‑1 (MQ=255)
         tATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTgg                                                                                                                                             <  2:1266282/150‑1 (MQ=255)
         tATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTgg                                                                                                                                             >  1:1033969/1‑150 (MQ=255)
         tATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTgg                                                                                                                                             <  1:915651/150‑1 (MQ=255)
          aTTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGt                                                                                                                                            <  2:303049/150‑1 (MQ=255)
           tttAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTGGTc                                                                                                                                           <  1:1583800/150‑1 (MQ=255)
           tttAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTc                                                                                                                                           <  1:347784/150‑1 (MQ=255)
            ttAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCg                                                                                                                                          <  1:578674/150‑1 (MQ=255)
                 tAAATATGGCGATGCAAATATCAATCATTATTAGAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATATGCTGACTGGTCTAACAt                                                                                                                                     >  1:1525440/1‑150 (MQ=255)
                      atGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGttt                                                                                                                                <  2:771402/150‑1 (MQ=255)
                        ggCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGc                                                                                                                              <  2:86635/150‑1 (MQ=255)
                        ggCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGc                                                                                                                              <  1:1502723/150‑1 (MQ=255)
                         gCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCt                                                                                                                             <  1:1386606/150‑1 (MQ=255)
                          cGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGACCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTa                                                                                                                            <  1:1200082/150‑1 (MQ=255)
                           gATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAATATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTaa                                                                                                                           >  1:1111969/1‑150 (MQ=255)
                             tGCAAATATCAAACATTATTATAGAAACTAAGAAATAAAACACTCACACCCCAAAAATAGAGCATCCGGCAAACGAGTGAGTTAACATCATTGTTAATAAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAAc                                                                                                                         <  1:1504121/150‑1 (MQ=255)
                              tcAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCAAAAATTACAGCATCCGGCAAAAGAGTGTGTTAACATCATTTTTAATTAATTCGACCAGTCAGGAAATTTGCTGACTGGTCGAACATAGTTTGCTAAACa                                                                                                                        <  1:305354/149‑1 (MQ=255)
                              acAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACa                                                                                                                        <  2:604736/149‑1 (MQ=255)
                                  atatCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtca                                                                                                                    >  2:771811/1‑150 (MQ=255)
                                   tatCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGCCCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGacat                                                                                                                   >  2:313758/1‑150 (MQ=255)
                                    atcaatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTAAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtcatc                                                                                                                  <  1:1622419/150‑1 (MQ=255)
                                      aaatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCac                                                                                                                <  2:1249609/149‑1 (MQ=255)
                                       aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAATGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGCCTGGTCGAACATAGTTTGCTAAACAGTCCTCact                                                                                                               >  2:1129813/1‑150 (MQ=255)
                                       aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCact                                                                                                               >  2:1673845/1‑150 (MQ=255)
                                        atcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactg                                                                                                              <  2:156133/150‑1 (MQ=255)
                                          caTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgac                                                                                                            <  2:28141/150‑1 (MQ=255)
                                          caTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgac                                                                                                            >  1:1197659/1‑150 (MQ=255)
                                          caTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgaa                                                                                                            >  2:306287/1‑149 (MQ=255)
                                           attattATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgact                                                                                                           >  1:226614/1‑150 (MQ=255)
                                           attattATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgact                                                                                                           >  2:544353/1‑150 (MQ=255)
                                            ttattaTAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCTAAAGAGTGAGTTAACATCATTGTTATTTAATTCGACCAGTCAGCAAATATGCTGTCTGGTGGAACATAGTTTTATTAACTGTCATAactgactg                                                                                                          >  1:325438/1‑150 (MQ=255)
                                               ttatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATg                                                                                                       <  1:20267/150‑1 (MQ=255)
                                               taataGAAACTTAGAATTAAAACACTCACAACACAAAATTATAGCATCCGGCAAAAGTGTGAGTTAACATCATTGTTAATTAATTCGACCCGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATg                                                                                                       <  2:1271710/148‑1 (MQ=255)
                                                 ataGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGTGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTTACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGAATGATgtc                                                                                                     <  2:624537/150‑1 (MQ=255)
                                                 ataGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtc                                                                                                     <  1:1360829/150‑1 (MQ=255)
                                                 ataGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtc                                                                                                     >  2:1095347/1‑150 (MQ=255)
                                                     aaaCTTATAAATAAAACCCTCACAACCCAAAAATAGAGCATCCTGCAAAAGATTGAGTTAACATCATTGTTAATTAATTCGACTATTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACATTCTTCACTGAATGATgtcagtc                                                                                                 <  2:79151/150‑1 (MQ=255)
                                                     aaaCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtcat                                                                                                                   >  2:1162211/1‑132 (MQ=255)
                                                     aaaCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtcat                                                                                                                   <  1:1162211/132‑1 (MQ=255)
                                                      aaCTTAGAAATAAAACAATAACAACCAAAAAATAGATCATCCGGCAAATGAGTGAGTTAACATCAATGTTAATTAATTCGACCAGTCATCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagtca                                                                                                <  1:312753/150‑1 (MQ=255)
                                                        cTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcaca                                                                                              >  1:519311/1‑150 (MQ=255)
                                                        cTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGGcaca                                                                                              >  2:1118757/1‑150 (MQ=255)
                                                          tAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACAcc                                                                                            >  2:152658/1‑150 (MQ=255)
                                                           aGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCa                                                                                           <  1:771811/150‑1 (MQ=255)
                                                           aGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATAAGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCa                                                                                           <  1:295530/150‑1 (MQ=255)
                                                              aaTAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATcc                                                                                        <  2:1469033/150‑1 (MQ=255)
                                                              aaTAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATcc                                                                                        >  2:118713/1‑150 (MQ=255)
                                                                 aaaaCACTCCCAACCCAAAAATAGAGCATCCGGCAAATGAGTGAGTTAACATCATTGTTAATTAATTCCACCAGTCAGCAAATATGCTGACTGGTAGCACATAGTTAGCTAACGAGTCAACACTGACTGATGTCAGGCACACCATCCGcc                                                                                     >  2:1676507/1‑150 (MQ=255)
                                                                 aaaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGcc                                                                                     <  1:895612/150‑1 (MQ=255)
                                                                 aaaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGATGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCATTCACACAATCCGcc                                                                                     >  1:1043788/1‑150 (MQ=255)
                                                                  aaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTTAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCa                                                                                    <  2:824534/150‑1 (MQ=255)
                                                                  aaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCa                                                                                    >  1:1208374/1‑150 (MQ=255)
                                                                       cTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCa                                                                               <  2:849294/150‑1 (MQ=255)
                                                                         cacaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCaaa                                                                             >  2:554218/1‑149 (MQ=255)
                                                                         cacaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACAATCCGCCAGTTCAAt                                                                             >  2:711835/1‑150 (MQ=255)
                                                                           caACCCAAAAATAGAGCATCCGGCAAAAGTGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCa                                                                           >  2:541172/1‑150 (MQ=255)
                                                                               ccAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCg                                                                       <  2:886153/150‑1 (MQ=255)
                                                                               ccAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCg                                                                       >  1:1014871/1‑150 (MQ=255)
                                                                                   aaaTAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATa                                                                   >  1:787794/1‑150 (MQ=255)
                                                                                    aaTAGATCATGAGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGCCCAGTCAGCAAATATGCAGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCCATCAGGCGAATaa                                                                  <  1:284349/150‑1 (MQ=255)
                                                                                            aTCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGcc                                                          <  1:1673845/150‑1 (MQ=255)
                                                                                            aTCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGcc                                                          >  1:1449024/1‑150 (MQ=255)
                                                                                            aTCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTACGcc                                                          >  1:1730988/1‑150 (MQ=255)
                                                                                               cGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTg                                                       >  2:180324/1‑150 (MQ=255)
                                                                                                ggCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGt                                                      <  1:708160/150‑1 (MQ=255)
                                                                                                  cAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTaa                                                    <  2:1693376/150‑1 (MQ=255)
                                                                                                  cAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGATAAACAGTCATCACTGACTGATTTCAGTCACACCAACCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTaa                                                    >  1:1040239/1‑150 (MQ=255)
                                                                                                         gtgagtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgtg                                             <  1:222047/150‑1 (MQ=255)
                                                                                                          tgagtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTgg                                            <  2:1294484/150‑1 (MQ=255)
                                                                                                            agtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTAAATCAGGCGAATAAATTCCGCCCTGTAACCATGTGGAt                                          >  1:1188707/1‑150 (MQ=255)
                                                                                                                aaCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                                      >  2:1340064/1‑150 (MQ=255)
                                                                                                                aaCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTCGACGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCCTGCGAATAAATTCCGCCCAGTAACCCTTTGGATCtgc                                      >  2:865025/1‑148 (MQ=255)
                                                                                                                 aCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTca                                     <  1:1095347/150‑2 (MQ=255)
                                                                                                                  catcatTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTccc                                    >  1:1369958/1‑150 (MQ=255)
                                                                                                                    tcatTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCt                                  <  2:1601563/150‑1 (MQ=255)
                                                                                                                         gTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGc                             <  1:609468/150‑1 (MQ=255)
                                                                                                                          ttaattaaTTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGcc                            <  2:808455/150‑1 (MQ=255)
                                                                                                                           taattaatTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCGTCCCCTTTTGCct                           >  2:1171143/1‑150 (MQ=255)
                                                                                                                            aattaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctg                          >  2:1403136/1‑150 (MQ=255)
                                                                                                                             attaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctgc                         <  1:1116400/150‑1 (MQ=255)
                                                                                                                                 agtCGACCAGTCAGCAAATATGCTGACTGGTCGAACATATTTTGCTAAACAGTCAGCACTGACTGATGTCAGTCACAACTTCTGACAGATCAATCAGGCGAATAATTTCCTCCCTGTAACCCTGTGGATCTACCCCTTTTGCctgctgtg                     <  2:211939/148‑3 (MQ=255)
                                                                                                                                    cGACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCACTTTTGCCTGCTGAGccc                  >  2:678267/1‑150 (MQ=255)
                                                                                                                                    cGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCCTCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGccc                  >  1:1281953/1‑150 (MQ=255)
                                                                                                                                    cGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCATCCCCTTTTGCCTGCTGAGccc                  >  2:1428588/1‑150 (MQ=255)
                                                                                                                                    cGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACGGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTACCCCTGTGGATCTTCCCCGTGTGCCTGCTGAGccc                  >  1:329939/1‑150 (MQ=255)
                                                                                                                                      aCCAGGCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCCATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAc                <  1:1405950/150‑1 (MQ=255)
                                                                                                                                       ccAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACt               >  2:1676273/1‑150 (MQ=255)
                                                                                                                                       ccAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAcg               >  1:914542/1‑149 (MQ=255)
                                                                                                                                         aGTCAGCAAATATGCTGACTGGTCGAAAATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCGACTGt             >  2:90108/1‑150 (MQ=255)
                                                                                                                                           tCAGCTAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCTCTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGttt           <  2:706640/150‑1 (MQ=255)
                                                                                                                                              gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGAt        >  1:914591/1‑150 (MQ=255)
                                                                                                                                              gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGAt        >  2:789712/1‑150 (MQ=255)
                                                                                                                                               cAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATc       <  1:783359/150‑1 (MQ=255)
                                                                                                                                                aaaTATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCAATGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATct      <  2:850072/150‑1 (MQ=255)
                                                                                                                                                   tatGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATctgct   <  2:787794/150‑1 (MQ=255)
                                                                                                                                                    atGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCATTTTGCCTGCTGAGCCCACTGTTTGATctgctg  <  2:274474/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCTG  >  NC_000913/2383811‑2384108

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: