Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,208,927 A→T 5.2% S65T (TCT→ACT)  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,9270AT 5.2% 345.7 / 4.8 114S65T (TCT→ACT) tfaEe14 prophage; putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base A (57/51);  new base T (2/4);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 114 total aligned reads displayed.

ATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCT  >  NC_000913/1208781‑1209075
                                                                                                                                                  |                                                                                                                                                     
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)
    tATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAg                                                                                                                                                <  1:1153040/150‑1 (MQ=255)
     aTGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGa                                                                                                                                               >  1:1307429/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)
        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:649697/1‑150 (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)
                        tgTGGTTGTTTTTTCAGGGTAATCACCCGGTTCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTTTAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtg                                                                                                                            <  1:748176/149‑1 (MQ=255)
                         gTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                                                                                                                           >  2:920610/1‑150 (MQ=255)
                             ttGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAc                                                                                                                       <  2:1436547/150‑1 (MQ=255)
                              tGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAcc                                                                                                                      <  2:599299/150‑1 (MQ=255)
                                  ttttCTGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTTCTATAGACCGTTTCACCGAGATTGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑ggg                                                                                                                 <  1:1210231/150‑1 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑ggg                                                                                                                 <  2:371022/150‑1 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑ggg                                                                                                                 >  1:315664/1‑150 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑ggg                                                                                                                 <  2:315664/150‑1 (MQ=255)
                                   tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGTTTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑gggg                                                                                                                <  1:148371/150‑1 (MQ=255)
                                   tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATTTGTACCA‑gggg                                                                                                                >  1:177135/1‑150 (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)
                                             aTCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACa                                                                                                      <  1:495877/150‑1 (MQ=255)
                                              tCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAg                                                                                                     <  1:1447814/150‑1 (MQ=255)
                                               cACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAgg                                                                                                    <  2:1307429/150‑1 (MQ=255)
                                               cACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAgg                                                                                                    <  2:1115071/150‑1 (MQ=255)
                                                    ggAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACAt                                                                                               >  1:1694169/1‑150 (MQ=255)
                                                       gCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGcc                                                                                            >  2:1042334/1‑150 (MQ=255)
                                                        cTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCg                                                                                           <  1:1366245/150‑1 (MQ=255)
                                                        cTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCg                                                                                           >  1:808418/1‑150 (MQ=255)
                                                         tttGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCgg                                                                                          <  1:54848/150‑1 (MQ=255)
                                                            gATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAAt                                                                                       >  1:1560109/1‑150 (MQ=255)
                                                                tCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCg                                                                                   <  2:1364089/150‑1 (MQ=255)
                                                                tCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCg                                                                                   >  1:1116511/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)
                                                                          ttctaGTTTTGGTGCTATATACCGTTTCACCGAGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGACTTATGTGTACCA‑GGTGAATCTAAACAGGAACATGCCGGAATGCCGACACCAACgg                                                                         <  1:84423/146‑1 (MQ=255)
                                                                          ctcCAGTTTTGGAGCTATAGACCTTTTCAACGCGATTGTCTGGAACATATTCCCATTAGTTAAAATCTTCATAACGCCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACgg                                                                         <  1:62885/150‑1 (MQ=255)
                                                                            cctgtattGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCAGATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAATGCATAACCACCCTTATGTGTACCA‑GGGGCAGCTAAACAGGCACATGCCGGAATGCCGACACCACCGGCa                                                                       >  2:828096/7‑150 (MQ=255)
                                                                            ccAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCa                                                                       >  2:27921/1‑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)
                                                                                       gCTATAGACCGTTTCACCTCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGTCAGATTGCATAACCATCCTTATGTGTACCA‑GGGGCATCTAAACAGTAACATGCCGGAATGCCGACACCAACGGCAAGATATTCAtt                                                            <  2:819829/150‑1 (MQ=255)
                                                                                         tataGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGTGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCAtttt                                                          >  1:56975/1‑149 (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:785933/1‑150 (MQ=255)
                                                                                             gACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGt                                                      <  2:287940/150‑1 (MQ=255)
                                                                                              aCCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTg                                                     >  2:1084375/1‑150 (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                                               <  2:1352112/150‑1 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAata                                               >  1:1622072/1‑150 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAata                                               <  1:1435779/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                    <  1:1128936/150‑1 (MQ=255)
                                                                                                                               cATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTta                    >  1:134647/1‑150 (MQ=255)
                                                                                                                               cATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTta                    >  2:860806/1‑150 (MQ=255)
                                                                                                                                    gTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCACGTGTTTCACCATCCTAGTTATAAAc               >  1:446265/1‑150 (MQ=255)
                                                                                                                                     atgcaaTCTGCCGTTCTGCAGATTTCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACg              >  2:347948/5‑150 (MQ=255)
                                                                                                                                     atgcaaTCTGCCGTTCTGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATTTATTCCCGTGTTTCACCATCGTAGTTATAAACg              >  2:1323720/5‑150 (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)
                                                                                                                                               cAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATccc    >  2:510236/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)
                                                                                                                                                  |                                                                                                                                                     
ATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCA‑GGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCT  >  NC_000913/1208781‑1209075

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