Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,208,930 G→T 5.3% R64S (CGT→AGT)  tfaE ← e14 prophage; putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,9300GT 5.3% 334.9 / 4.8 114R64S (CGT→AGT) tfaEe14 prophage; putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base G (56/51);  new base T (2/5);  total (58/56)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 115 total aligned reads displayed.

ATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCTGCC  >  NC_000913/1208781‑1209078
                                                                                                                                                     |                                                                                                                                                     
aTCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACg                                                                                                                                                       >  2:536723/1‑150 (MQ=255)
   gTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCa                                                                                                                                                    >  1:1199528/1‑150 (MQ=255)
    tATGGAGATAAAGGGGCGTTAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTGTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAg                                                                                                                                                   <  1:190815/150‑1 (MQ=255)
    tATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAg                                                                                                                                                   >  1:882357/1‑150 (MQ=255)
      tGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAt                                                                                                                                                 >  2:1128936/1‑150 (MQ=255)
        gagaTATAGGGGCGATAGTGGTTGTATTTTGAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTg                                                                                                                                               <  1:62706/150‑1 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCCGTTCTGCAGAttt                                                                                                                                               <  2:1237089/150‑2 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCCGTTCTGCAGAttt                                                                                                                                               <  2:629536/150‑2 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTg                                                                                                                                               <  2:961702/150‑1 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTg                                                                                                                                               >  2:754841/1‑150 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCGTTGATTTCATTTGATTCTCCAGTTTTGGTGCTATAGACCTTTTCACCGCTATTGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTg                                                                                                                                               <  2:1061969/150‑1 (MQ=255)
         agaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCCGTTCTGCAGAtttt                                                                                                                                              <  1:507036/150‑3 (MQ=255)
         agaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGc                                                                                                                                              <  1:768943/150‑1 (MQ=255)
          gaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCa                                                                                                                                             <  2:1517376/150‑1 (MQ=255)
          gaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCa                                                                                                                                             <  2:854406/150‑1 (MQ=255)
            tAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATTGTCTGGCACATATTCCCATGAGTTATAATCTGCCGTTCTGCAGAttttaac                                                                                                                                           <  2:212927/150‑6 (MQ=255)
            tAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATa                                                                                                                                           >  1:794286/1‑150 (MQ=255)
             aaaGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATaa                                                                                                                                          >  2:912686/1‑150 (MQ=255)
                ggggCGATAGTTGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGATTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCa                                                                                                                                       <  1:1395949/150‑1 (MQ=255)
                 gggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAg                                                                                                                                      >  1:207223/1‑150 (MQ=255)
                  ggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGc                                                                                                                                     >  1:1168577/1‑150 (MQ=255)
                     gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                                                                                                                                  <  1:1225131/150‑1 (MQ=255)
                     gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                                                                                                                                  >  1:757673/1‑150 (MQ=255)
                        tgTGGTTGTTTTTTCAGGGTAATCACCCGGTTCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTTTAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtg                                                                                                                               <  1:748176/149‑1 (MQ=255)
                        aGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATGCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtg                                                                                                                               <  2:1207936/150‑1 (MQ=255)
                         gTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                                                                                                                              >  2:920610/1‑150 (MQ=255)
                          tGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAATGCATAACCAGCCTTAtgtg                                                                                                                             <  1:1160850/150‑1 (MQ=255)
                            gTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTa                                                                                                                           <  2:457112/150‑1 (MQ=255)
                              tGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAcc                                                                                                                         <  2:599299/150‑1 (MQ=255)
                                tATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATTTTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑g                                                                                                                      >  2:184020/1‑150 (MQ=255)
                                  ttttCTGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTTCTATAGACCGTTTCACCGAGATTGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑ggg                                                                                                                    <  1:1210231/150‑1 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑ggg                                                                                                                    <  2:315664/150‑1 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑ggg                                                                                                                    <  2:371022/150‑1 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑ggg                                                                                                                    >  1:315664/1‑150 (MQ=255)
                                   tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGTTTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑gggg                                                                                                                   <  1:148371/150‑1 (MQ=255)
                                   tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑gggg                                                                                                                   <  1:1028461/150‑1 (MQ=255)
                                    ttgaGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGc                                                                                                                  <  1:1149145/147‑1 (MQ=255)
                                     tCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTAACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCa                                                                                                                 <  2:741669/150‑1 (MQ=255)
                                            aaTCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGTTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATTAGTTAAAATCTTCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAAc                                                                                                          <  1:1322836/150‑1 (MQ=255)
                                             aTCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACa                                                                                                         <  1:495877/150‑1 (MQ=255)
                                               cACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAgg                                                                                                       <  2:1115071/150‑1 (MQ=255)
                                                  ccGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAAc                                                                                                    >  2:1024634/1‑150 (MQ=255)
                                                    ggAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACAt                                                                                                  >  1:1694169/1‑150 (MQ=255)
                                                       gCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGcc                                                                                               >  1:472078/1‑150 (MQ=255)
                                                       gCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGcc                                                                                               >  2:1042334/1‑150 (MQ=255)
                                                        cTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCg                                                                                              >  1:808418/1‑150 (MQ=255)
                                                        cTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCg                                                                                              <  1:1366245/150‑1 (MQ=255)
                                                         tttGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCgg                                                                                             <  1:54848/150‑1 (MQ=255)
                                                            gATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAAt                                                                                          >  1:1560109/1‑150 (MQ=255)
                                                                tCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCg                                                                                      >  1:1116511/1‑150 (MQ=255)
                                                                tCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCg                                                                                      <  2:1364089/150‑1 (MQ=255)
                                                                 attttGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGa                                                                                     <  2:1411924/149‑1 (MQ=255)
                                                                  ttttGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGac                                                                                    >  1:1493198/1‑150 (MQ=255)
                                                                  ttttGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCACCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGaa                                                                                    >  1:1034191/1‑149 (MQ=255)
                                                                     tgtttCTCCATTTTTGGTTCTATTGACCGTTTCCCCGCGTTGGTCTGGCCCATATTCCCATTAGTTAAAATTTGCAGAAAGGCAGATTGCATTACCAGGCTTATGTGTAACA‑GGGGCATATAAACAGGAACATGCCGGAATGCCGACAcc                                                                                 <  1:1102019/147‑1 (MQ=255)
                                                                         tctcCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACg                                                                             >  1:1287952/1‑150 (MQ=255)
                                                                          ctcCAGTTTTGGAGCTATAGACCTTTTCAACGCGATTGTCTGGAACATATTCCCATTAGTTAAAATCTTCATAACGCCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACgg                                                                            <  1:62885/150‑1 (MQ=255)
                                                                            cctgtattGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCAGATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAATGCATAACCACCCTTATGTGTACCA‑GGGGCAGCTAAACAGGCACATGCCGGAATGCCGACACCACCGGCa                                                                          >  2:828096/7‑150 (MQ=255)
                                                                            ccAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCa                                                                          >  2:1620477/1‑150 (MQ=255)
                                                                                 tttGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGata                                                                     <  2:1425510/150‑1 (MQ=255)
                                                                                      tGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACAACAACGGCAAGATATTCAt                                                                >  2:1030513/1‑150 (MQ=255)
                                                                                         tataGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTg                                                             >  2:704164/1‑150 (MQ=255)
                                                                                            aGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                                          >  2:430175/1‑150 (MQ=255)
                                                                                            aGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                                          <  2:1554833/150‑1 (MQ=255)
                                                                                             gACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                                          >  2:1333855/1‑150 (MQ=255)
                                                                                             gACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGt                                                         <  2:794286/150‑1 (MQ=255)
                                                                                             gACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGt                                                         >  2:785933/1‑150 (MQ=255)
                                                                                             gACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGt                                                         <  2:287940/150‑1 (MQ=255)
                                                                                                 gTTTCACCGCGATGGTCTGGCACATATTCCCAAGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGaa                                                     >  2:180212/1‑150 (MQ=255)
                                                                                                   ttCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAat                                                   <  2:750478/150‑1 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGCCACCACCGGCAAGATATTCATTTGAAGTGGAAata                                                  >  1:179505/1‑150 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAata                                                  >  1:1622072/1‑150 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAata                                                  <  1:1435779/150‑1 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAata                                                  <  2:1352112/150‑1 (MQ=255)
                                                                                                     cACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAatat                                                 >  2:578777/1‑150 (MQ=255)
                                                                                                       ccGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAatatat                                               >  2:1178790/1‑150 (MQ=255)
                                                                                                          cgATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTcc                                            <  1:362138/150‑1 (MQ=255)
                                                                                                            aTGGTCTGGCACATATACCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCg                                          >  2:1368896/1‑150 (MQ=255)
                                                                                                             tGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCgt                                         <  2:757673/150‑1 (MQ=255)
                                                                                                               gTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCgtgt                                       >  1:1256781/1‑150 (MQ=255)
                                                                                                                tCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACTGCAGATTGCATAACCATCATTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGtt                                      <  1:1521697/150‑1 (MQ=255)
                                                                                                                   ggCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATACCGACACCAACGGCAAGATATTCATTTGAATTGGAAATATATTCCCGTGTTTCa                                   >  1:878079/1‑150 (MQ=255)
                                                                                                                     cacaTATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATACCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCAcc                                 >  1:994284/1‑150 (MQ=255)
                                                                                                                      acaTATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCa                                >  2:1114365/1‑150 (MQ=255)
                                                                                                                       caTATTACCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGCCACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCAt                               >  2:549204/1‑150 (MQ=255)
                                                                                                                          atTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGt                            >  1:452202/1‑150 (MQ=255)
                                                                                                                            tCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAg                          <  2:1694169/150‑1 (MQ=255)
                                                                                                                               cATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTta                       >  2:860806/1‑150 (MQ=255)
                                                                                                                               cATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTta                       >  1:134647/1‑150 (MQ=255)
                                                                                                                               cATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTta                       <  1:1128936/150‑1 (MQ=255)
                                                                                                                                     atgcaaTCTGCCGTTCTGCAGATTTCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACg                 >  2:347948/5‑150 (MQ=255)
                                                                                                                                     atgcaaTCTGCCGTTCTGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATTTATTCCCGTGTTTCACCATCGTAGTTATAAACg                 >  2:1323720/5‑150 (MQ=255)
                                                                                                                                       aaaaTCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGt               <  1:1896/150‑1 (MQ=255)
                                                                                                                                         aaTCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGACCATGCCGGAATGCCGACACGAACGGCAAGATATTCATTTGAAGTGGAAAGATATTCCCGTGTTTCACCCTCGTAGTTTTAAACGGTaa             >  2:572734/1‑150 (MQ=255)
                                                                                                                                           tCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCTCGTGTTTCACCATCGTAGTTATAAACGGTAACa           <  1:1042334/150‑1 (MQ=255)
                                                                                                                                             tGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCCCCATCGTAGTTATAAACGGTAACATc         >  1:1442661/1‑150 (MQ=255)
                                                                                                                                                aGAACGGCAGATTTCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCATAGTTATAAACGGTAACATcccc      <  2:556115/150‑1 (MQ=255)
                                                                                                                                                aGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATcccc      >  1:348721/1‑150 (MQ=255)
                                                                                                                                                 gAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGCCACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCt     >  1:1617533/1‑150 (MQ=255)
                                                                                                                                                    cGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCTGcc  <  1:536723/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCTGCC  >  NC_000913/1208781‑1209078

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